Решение системы нелинейных уравнений. - C (СИ) - Ответ 2183049

#include <cstdlib>
#include <iostream>
#include <math.h>
 
 
#if 0  // код в этом блоке компилятор не воспринимает, так как препроцессор его сьедает, можно здесь писать всякую чепуху не по делу и ничего не бояться
xy-tg(x-y)=0
0,5x^2+2y^2-1=0
 
x=0,975247
y=0,512078     
#endif
 
double fun1(double x,double y)
{
    double f;
    f=x*y-tan(x-y); //здесь твоя функция №1
    return f;
}
 
double fun2(double x,double y)
{
    double f;
    f=0.5*x*x+2*y*y-1; //здесь твоя функция №2
    return f;
}
 
double dfun1_dx(double x,double y)
{
    double f;
    f=(y-1)*x-(tan(x-y)*tan(x-y)+1); //здесь твоя производная от функции №1 по х
    return f;
}
 
double dfun1_dy(double x,double y)
{
    double f;
    f=(x-1)*y+(tan(x-y)*tan(x-y)+1);  //здесь твоя производная от функции №1 по y
    return f;
}
 
double dfun2_dx(double x, double y)
{
    double f;
    f=x;   //здесь твоя производная от функции №2 по х
    return f;
}
 
double dfun2_dy(double x,double y)
{
    double f;
    f=4*y;   //здесь твоя производная от функции №2 по y
    return f;
}
 
using namespace std;
 
int main()
{
    double x,y,e,x0,y0,a,b,c,d,lick,suck,dx,dy; // переменные, все те же что и у тебя только без массивов ( массивов в с++ нету как таковых, но можно сделать но это надо отдельно обьяснять)
    
    x=0.9; // (приближенное х)
    y=0.5; // (приближенное у)
    e=0.000001; //точность
 
    do
    {
        x0=x;
        y0=y;
        lick=-fun1(x0,y0);
        suck=-fun2(x0,y0);
        a=dfun1_dx(x0,y0);
        b=dfun1_dy(x0,y0);
        c=dfun2_dx(x0,y0);
        d=dfun2_dy(x0,y0);
        dx=(lick*d-b*suck)/(a*d-b*c);
        dy=(a*suck-lick*c)/(a*d-b*c);
        x=x0+dx;
        y=y0+dy;
    } while (((dx<e && dx>(-e)) || (dy<e && dy>(-e)))==false);  // в этом цикле все абсолютно в точности как у тебя, ничего менять не нужно (кстати в твоем коде ошибки, здесь исправлены)
 
    printf ("x = %f\n",x);
    printf ("y = %f\n",y);
    system("PAUSE");
    return 1;
}