Написать программу в с++ решения линейных уравнений методом неполной релаксации ! )

#include< iostream.h> #include< math> int maximal(int n, double r0[]){ int i,f; f=0.0; for(i=0; i< n-1; i++){ if(r0[i+1]> r0[i]) f=i+1; } return f; } void main(){ int i,j,n,f,k,iter; double s,det; cout< < "vvedite razmer kvadratnoy matrici= "; cin> > n; double *x=new double [n]; double **b=new double *[n]; for(i=0; i< n; i++)   b[i]=new double[n+1]; double **a=new double *[n]; for(i=0; i< n; i++)   a[i]=new double[n+1]; cout< < "vvedite kolichestvo iteraciy: "; cin> > iter; cout< < "vvedite matritcu"; for(i=0; i< n; i++){ for(j=0; j< =n; j++) cin> > b[i][j]; } cout< < "podgotovka k "; for(i=0; i< n; i++){ for(j=0; j< n; j++) a[i][j]=-b[i][j]/b[i][i]; a[i][n]=b[i][n]/b[i][i]; } for(i=0; i< n; i++){ for(j=0; j< n+1; j++) cout< < " "< < a[i][j]< < " || "; cout< < "\n"; } double *x0=new double [n]; for(i=0; i< n; i++) x[i]=0.0; double *r0=new double [n]; cout< < "vvedite znachenie nachal`nih priblizheniy: \n"; for(i=0; i< n; i++) cin> > x0[i]; s=0.0; for(i=0; i< n; i++){ for(j=0; j< n; j++) s=s+a[i][j]*x0[i]; } for(i=0; i< n; i++){ r0[i]=a[i][n]-x0[i]+s; cout< < "r("< < i< < ")="< < r0[i]< < " | "; } f=maximal(n,r0); det=r0[f]; for(k=0; k< iter; k++){ cout< < "det{"< < k< < "}="< < det< < "\n"; for(i=0; i< n; i++){ if(i! =f) r0[i]=r0[i]+a[i][f]*det; else r0[i]=r0[i]-det; } for(i=0; i< n; i++) cout< < "r["< < i+1< < "]="< < r0[i]< < " "; x[f]=x[f]+det; f=maximal(n,r0); det=r0[f]; } cout< < "\n"; for(i=0; i< n; i++) cout< < "x{"< < i+1< < "}="< < x[i]< < "\n"; delete []x; delete []r0; delete []x0; delete []a; cin.get(); cin.get(); }

40 гигабайт это  40960 мб  40960/600 = 68  ответ 68 дисков