矩阵与数值分析实验报告(2)

大连理工大学矩阵与数值分析课又叫计算机科学计算的上机实验报告,内容详尽,有程序代码和预算结果

cout<<"结果为："<<s2<<endl;

cout<<"N=1000000;S->B:"<<endl;

for(double j=2;j<=n3;j++)

{

s3 += (double)1/(j*j-1);

}

cout<<"结果为："<<s3<<endl;

s3=0;

cout<<"N=1000000;B->S:"<<endl;

for(double j=1000000;j>=2;j--)

{

s3 += (double)1/(j*j-1);

}

cout<<"结果为："<<s3<<endl;

cin>>a;

return 0;

}

计算结果：

结果分析：

二、解线性方程组

1. 分别用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组

2

1

0 0 1 21001 210 0 1 2

1 j 3 x1 x2 x 3 x 4 1 0 = 0 0 迭代法计算的停止条件为：maxx(jk 1) x(jk) 10 6。

Jacobi方法代码：

Jacobi方法计算结果：

Jacobi方法结果分析：

Gauss-Seidel方法代码：