Задача линейного,и нелинейного уравнения. - C++ - Ответ 2333627

#include <stdio.h>
#include <conio.h>
#include <math.h>
#define MAX 100
main ()
{
float a[MAX] [MAX],b[MAX],x[MAX],h,eps=1E-6,delta_x[MAX];
int i,j,n=3,k;
 
 printf("n = ");
scanf("%i",&n);
for (i=1; i<=n; i++)
{
  printf("x[%i] = ",i);
  scanf("%g",&x[i]);
}
do
{
a[1][1]=cos(x[1]);
a[2][1]= exp(x[1]);
a[3][1]= -sin(x[1]);
a[1][2]= sin(x[2]);
a[2][2]= cos(x[2]);
a[3][2]=-4*(pow(x[2],3));
a[1][3]=4*(pow(x[3],3);
a[2][3]=-1;
a[3][3]= 4;
b[1]= sin(x[1])-5*cos(x[2])+pow(pow(x[3],4)-9;
b[2]= exp(x[1])+sin(x[2])-pow(x[3]-23;
b[3]= cos(x[1])-(pow(x[2],4)+4*pow(x[3]-5;
 
 for(i=1;i<=n-1;i++)
 {
 for(j=i+1;j<=n;j++)
 {
  a[j][i]=-a[j][i]/a[i][i];
    for(k=i+1;k<=n;k++)
    a[j][k]+=a[j][i]*a[i][k];
    b[j]+=a[j][i]*b[i];
 }
 x[n]=b[n]/a[n][n];
 
   for(i=n-1;i>=1;i--)
   {
   h=b[i];
 
 for(j=i+1;j<=n;j++)
 h-=x[j]*a[i][j];
  x[i]=h/a[i][i];
   }
    for(i=n;i<=n;i++)
    {
 h+=fabs(delta_x[i]);
 x[i]=x[i]+delta_x[i];
    }
 }
 while(n>eps);
 
for(i=1;i<=n;i++)
printf("x%i = %f\n",i,x[i]);
getch ();
return 0;
}