помогите найти ошибку (си выч.мат)

[Dimoff80]

не могу найти ошибку:решить снау методом ньютона
sin(x+y)-1.6x
x^2+y^2=1
подскажите пожалуйста!

#inclue<stdio.h>
# include <math.h>
const e=0.001;
main()
{
int n=0;
float x,y,xn,yn,dW,f1,f2;
float wx1,wx2,wy1,wy2;
printf("\n\t Vvedite x y \n");
scanf("%f%f",&xn,&yn);
do
{
n++;
x=yn;
y=yn;
f1=sin(x+y)-1.6*x;
f2=(x*x+y*y)-1;
wx1=cos(x+y);
wx2=1;
wy1=2*x;
wy2=2*y;
dW=wx1*wy2-wx2*wy1;
if (dW>0){printf("stop\n");}
 xn=x-((f1*wy2+wx2*f2)/dW);
xn=y-((wx1*f2+f1*wy1)/dW);
}while ((fabs(xn)>e && fabs(yn)>e));
printf ( "\n x=%f y=%f n=%d",xn,yn,n);
}

[Svetlag]

Может нужно так:

do
{
n++;
x=xn;
y=yn;
f1=sin(x+y)-1.6*x;
f2=(x*x+y*y)-1;
wx1=cos(x+y);
wx2=1;
wy1=2*x;
wy2=2*y;
dW=wx1*wy2-wx2*wy1;
if (dW>0){printf("stop\n");}
 xn=x-((f1*wy2+wx2*f2)/dW);
yn=y-((wx1*f2+f1*wy1)/dW);
}while ((fabs(xn)>e && fabs(yn)>e));

[Dimoff80]

извини старый текст дал. компилятор выдаёт Invalid floating point operation

[Svetlag]

у меня вот это компилируется в Visual C++. Но выдает некорректные результаты.Напишите формулы,по которым патаетесь написать программу. Тогда можно и говорить.

#include<stdio.h>
#include<math.h>
const double e=0.001;
int main()
{
	int n=0;
	float x,y,xn,yn,dW,f1,f2;
	float wx1,wx2,wy1,wy2;
	printf("\n\t Vvedite x y \n");
	scanf("%f%f",&xn,&yn);
	do
	{
	n++;
	x=xn;
	y=yn;
	f1=sin(x+y)-1.6*x;
	f2=(x*x+y*y)-1;
	wx1=cos(x+y);
	wx2=1;
	wy1=2*x;
	wy2=2*y;
	dW=wx1*wy2-wx2*wy1;
	if (dW>0){printf("stop\n");}
	 xn=x-((f1*wy2+wx2*f2)/dW);
	yn=y-((wx1*f2+f1*wy1)/dW);
	}while ((fabs(xn)>e && fabs(yn)>e));
	printf ( "\n x=%f y=%f n=%d",xn,yn,n);
	return 0;
}