Добавил:
Maestro
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:LAB3-4
.H#include <stdio.h>
#include <math.h>
#include <conio.h>
void main()
{
int n=0;
float al=-0.1, bt=-0.4, gm=-0.1, dt=0.6, eps=0.001;
float x0=3, y0=2;
float x,x1,z1;
float y,y1,z2;
x=x0;
y=y0;
clrscr();
do
{
x1=x+(al*(x*x+y*y-4))+(bt*(x-y-2)); //MPI
//x1=x-((-x*x+y*y-2*x*y+4*y+4)/(-2*x-2*y));
z1=x;
y1=y+(gm*(x*x+y*y-4))+(dt*(x-y-2));
//y1=y-((x*x-y*y-2*x*y-4*x+4)/(-2*x-2*y));
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);
getch();
}