#include<stdio.h>

#include<math.h>

main()

{ int n=0;

float al=-5/18, bt=0, gm=5/54, dt=-1/3, eps=0.00001;

float x0=1.8,y0=1;

float x,x1,z1;

float y,y1,z2;

x=x0;

y=y0;

printf("n++ x x1 fabs(x1-x) y y1 fabs(y1-y)\n");

do

{

x1=x+(al*(x*x-4) + bt*(x+3*y-8)); // MPI должно вывести х=2, у=2

//x1=x+((-1/(2*x))*(x*x-4)+ (0*(x+3*y-8))); //Oбщий метод должно вывести х=2, у=2

//x1=x-((x*x-4)/3);//Метод Ньютона работает корректно

z1=x;

y1=y+(gm*(x*x-4) + dt*(x+3*y-8)); // MPI должно вывести х=2, у=2

//y1=y+((1/(6*x))*(x*x-4)-(1/3)*(x+3*y-8)); //Oбщий метод должно вывести х=2, у=2

//y1=y-((x*x+3*x+9*y-28)/(6*x));//Метод Ньютона работает корректно

z2=y;

printf("%d %.4f %.4f %.4f %.4f %.4f %.4f\n",n++,x,x1,fabs(x1-x),y,y1,fabs(y1-y));

x=x1;

y=y1;

}

while(fabs(z1-x)>eps || fabs(z2-y)>eps);

}

edited Mar 12 by Rail4icK