Решение уравнения перебором

int X_1 = StrToFloat(Edit1->Text);
int Y_1 = StrToFloat(Edit2->Text);
int Z_1 = StrToFloat(Edit3->Text);
int x;
int y;
int z;
 
for (x = 0; x<360; x=x+1)
    for (y = 0; y<360; y=y+1)
        for (z = 0; z<360; z=z+1)
 
 
 
r = sqrt((pow((X_1-x),2)+ pow((Y_1-y),2) + pow((Z_1-z),2)));
Memo1->Lines->Strings[0] = FloatToStr(x);
Memo1->Lines->Strings[1] = FloatToStr(y);
Memo1->Lines->Strings[2] = FloatToStr(z);

#include <vcl.h>
#pragma hdrstop
#include <math.h>
//#include <stdio.h>
//#include <iostream>
#include "Us0.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
double Xmin=-20,Xmax=20,Ymin=-20,Ymax=20,Zmin=-20,Zmax=20;
double step =1.;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button1Click(TObject *Sender)
{
int X1 = 105, X2=10, X3=-3;     //StrToFloat(Edit1->Text);
int Y1 = -5,  Y2=100,  Y3=0;    //StrToFloat(Edit2->Text);
int Z1 = 10,  Z2=2, Z3=-111;    //StrToFloat(Edit3->Text);
double R1=100, R2=105, R3=110,r1,r2,r3, S1,S2,S3;
double x,y,z, a=0,eps,xmin=R3,xmax=-R3,ymin=R3,ymax=-R3,zmin=R3,zmax=-R3;
eps=(R3*8*step);
r1=R1*R1; r2=R2*R2; r3=R3*R3;
Edit1->Text= FloatToStr(eps);
for (x=Xmin; x<Xmax; x=x+step)
  for (y = Ymin; y < Ymax; y=y+step)
     for (z = Zmin; z< Zmax; z=z+step)
       {
        S1=(X1-x)*(X1-x)+ (Y1-y)*(Y1-y) + (Z1-z)*(Z1-z);
        S2=(X2-x)*(X2-x)+ (Y2-y)*(Y2-y) + (Z2-z)*(Z2-z);
        S3=(X3-x)*(X3-x)+ (Y3-y)*(Y3-y) + (Z3-z)*(Z3-z);
        if((fabs(r1-S1)< eps) && (fabs(r2-S2)< eps) && (fabs(r3-S3)< eps))
          {
           if(x<xmin)xmin=x; else if(x>xmax)xmax=x;
           if(y<ymin)ymin=y; else if(y>ymax)ymax=y;
           if(z<zmin)zmin=z; else if(z>zmax)zmax=z; a=1;
          }
       }
Xmin=xmin; Xmax=xmax; Ymin=ymin; Ymax=ymax; Zmin=zmin; Zmax=zmax; step=step/2;
if(a)
 {
  Label1->Caption="xmin ="+FloatToStr(Xmin)+"  ymin ="+FloatToStr(Ymin)+"  zmin ="+FloatToStr(Zmin);
  Label2->Caption="xmax ="+FloatToStr(Xmax)+"  ymax ="+FloatToStr(Ymax)+"  zmax ="+FloatToStr(Zmax);
 }
  else {Label3->Visible=1; Label3->Caption = "xi-xi";}
}