СНУ на СИ

[peace on you]

Привет всем, выручайте пожалуйста, такая трабла.
Си выделяет синим следующую строку и выдает ошибку(прикрепил).
Где математическая ошибка может быть? Кажется до того она мелкая, но глаза уже замазолились.
- расчеты

#include <stdio.h>
#include <conio.h>
#include <math.h>
void main()
{
float f1,f2,e=0.001,x0=0.25,y0=0.5,x,y,k,n,pr,
      dx,dy,deltax,deltay,a,b,c,d;
met1: f1=cos((x0-y0)/3)-2*y0;
      f2=sin((x0+y0)/3)-2*x0;
a=-1/3*sin((x0-y0)/3); b=1/3*sin((x0-y0)/3)-2;
c=1/3*cos((x0+y0)/3)-2; d=1/3*cos((x0+y0)/3);
k=(b*f2*a)/d-f1/a;
n=1-(b*c*a)/d;
dx=k/n;
dy=(-f2-c*dx)/d;
x=x0+dx;
y=y0+dy;
deltax=fabs(dx);
deltay=fabs(dy);
if (deltax>e && deltay>e)
{x0=x,y0=y;
goto met1;}
else {printf("x=%f  y=%f\n",x,y);
 pr=sin((x0+y0)/3)-2*x0;
 printf("pr=%f\n",pr);}
getch( );
}

[Mandrivnyk]

d по ходу дела равно 0.

а все потому, что

a = 1/3;
b = 1./3.;

дадут совершенно разные значения.

И это... плюнь в глаз тому, кто научил использовать в С оператор goto...

[peace on you]

не, d не может быть равно 0, там же производные везде берутся.

[Mandrivnyk]

Сообщение от peace on you

не, d не может быть равно 0, там же производные везде берутся.

а ты вставь после строчек с вычислением вывод этих четырех переменных.
очень удивишься.

вот так, например:

#include <stdio.h>
#include <conio.h>
#include <math.h>
int main()
{
float f1,f2,e=0.001,x0=0.25,y0=0.5,x,y,k,n,pr,
      dx,dy,deltax,deltay,a,b,c,d;
    f1=cos((x0-y0)/3)-2*y0;
    f2=sin((x0+y0)/3)-2*x0;
    a=-1/3*sin((x0-y0)/3);
    b=1/3*sin((x0-y0)/3)-2;
    c=1/3*cos((x0+y0)/3)-2;
    d=1/3*cos((x0+y0)/3);
    printf("%f %f %f %f\n", a, b, c, d);
    getch( );
}

но если ты строчку

a=-1/3*sin((x0-y0)/3)

напишешь, например, так:

a=-sin((x0-y0)/3)/3;

то все кардинально изменится.
Понимаешь, что именно я хочу объяснить?

[peace on you]

код выдал что а и d равны нулю?! ничего не понимаю... мне всего то надо что бы проверка была равна нулю. Прорешивал сейчас другой вариант, всё получилось сразу... а мой косячит...

[peace on you]

Сообщение от Mandrivnyk

Понимаешь, что именно я хочу объяснить?

кажется понимаю, нам что то говорили про это.
Но есть одно но, я поставил вот так

a=-1.0/3.0*sin((x0-y0)/3); b=1.0/3.0*sin((x0-y0)/3)-2;
c=1.0/3.0*cos((x0+y0)/3)-2; d=1.0/3.0*cos((x0+y0)/3);

и теперь он синим обвел

else {printf("x=%f  y=%f\n",x,y);

з.ы. я даже поставил 0,33

[peace on you]

up
всё... нет вариантов более исправления этой ошибки?