数值分析实验报告(3)

数值分析实验报告

cout<<\迭代法解第二个方程:\ jacobi(p[1],b,X2,x2,8);

cout<<\迭代法解第三个方程:\ jacobi(p[2],c,X3,x3,10);

cout<<\迭代法解第一个方程:\ GS(p[0],a,X1,x1,10);

cout<<\迭代法解第二个方程:\ GS(p[1],b,X2,x2,8);

cout<<\迭代法解第三个方程:\ GS(p[2],c,X3,x3,10);

cout<<\迭代法解第一个方程:\ SOR(p[0],a,X1,x1,10);

cout<<\迭代法解第二个方程:\ SOR(p[1],b,X2,x2,8);

cout<<\迭代法解第三个方程:\ SOR(p[2],c,X3,x3,10); }

int jacobi(float *p,float b[],float X[],float x[],int n) {

int k,i,j; float m,R,r,e;

cout<<\请输入精度e:\ cin>>e;

for(k=0;k

R=0;

for(i=0;i

m=0;

for(j=0;j

m=m+(*(p+i*n+j))*X[j]; x[i]=X[i]+(b[i]-m)/(*(p+i*n+i)); r=x[i]-X[i]; if(r<0)

r=X[i]-x[i]; if(r>R) R=r; }

if(R<=e&&R>0) {

print(x,n);

cout<<\迭代次数为:\ return k; }

for(j=0;j<10;j++) X[j]=x[j]; }

11

// jacobi 数值分析实验报告

print(x,n);

cout<<\迭代次数为:\

cout<<\方程解发散,无法用Jacobi方法解此方程!\ return 0; }

int GS(float *p,float b[],float X[],float x[],int n) // GS { int i,j,k; float t,R,r,e;

cout<<\请输入精度e:\ cin>>e;

for(k=0;k

for(i=0;i

for(j=0;j

if(ji) t+=(*(p+i*n+j))*X[j]; }

x[i]=(b[i]-t)/(*(p+i*n+i)); }

for(i=0;i<10;i++) {

r=x[i]-X[i]; if(r<0)

r=X[i]-x[i]; if(r>R) R=r; }

if(R<=e&&R>0) {

print(x,n);

cout<<\迭代次数为:\ return k; }

for(j=0;j<8;j++) X[j]=x[j]; } print(x,n);

cout<<\迭代次数为:\

cout<<\方程解发散,无法用Gauss-Seidel方法解此方程!\ return 0; }

12

数值分析实验报告

int SOR(float *p,float b[],float X[],float x[],int n) // SOR {

int i,j,k;

float t,R,r,e,w;

cout<<\请输入松弛因子w(0>w;

cout<<\请输入精度e:\ cin>>e;

for(i=0;i

R=0;

for(i=0;i

t=0;

for(j=0;j