Решение уравнений

//ïðîãà ðåøàåò íåëèíåéíîå óðàâíåíèå ìåòîäîì ñåêóùèõ
#include<stdlib.h>
#include<iostream>
#include<conio.h>
#include<math.h>
 
using namespace std;
 
int main()
{
 //clrscr();
 double x,x1,y,y1,X[5],Y[5],s;
 int n,i;
 x=-40;
 y=0.4+ atan(sqrt(x))-x;
  n=1;
  cout.setf(ios::fixed);
   cout.setf(ios::showpoint);
   cout.precision(10);
 for(x1=0;x1<=40;x1++)
{
  y1=0.4+ atan(sqrt(x1))-x1;
  if((y*y1)<0)
  {
   X[n]=x1;
   Y[n]=y1;
   n++;
  }
   y=y1;
   x=x1;
}
 X[0]=(X[1]-1);
 Y[0]=0.4+ atan(sqrt(X[0]))-X[0];
n=1;
 double X1[]={X[0], X[1], X[1], X[2]};
 double Y1[]={Y[0], Y[1], Y[1], Y[2]};
 do     
{
  x=X1[0] - Y1[0]*((X1[1]-X1[0])/(Y1[1]-Y1[0]));
  y=0.4+ atan(sqrt(x))-x;
 
  if((Y1[0]*y)>1)
{
  X1[0]=x;
  Y1[0]=y;
}
 else
{
  X1[1] = x;
  Y1[1] = y;
}
 n++;
s=y;
}while(fabs(s)>0.0001);
 cout << "x1="  << x << "   " << "y1=" << y << "   " << "n1=" << n <<endl;
 s=n=1;
 /*do     
{
  x=X1[2] - Y1[2]*((X1[3]-X1[2])/(Y1[3]-Y1[2]));
  y=0.4+ atan(sqrt(x))-x;
 
  if((Y1[2]*y)>1)
{
  X1[2]=x;
  Y1[2]=y;
}
 else
{
  X1[3] = x;
  Y1[3] = y;
}
 n++;
s=y;
}while(fabs(s)>0.001);//   здесь задается точность
 cout << "x2="  << x << "   " << "y2=" << y << "   " << "n2=" << n <<endl;*/
 s=n=1;
getch();
return 0;
}