蜂蜜柚子茶的博客

完整代码如下： #include #include<math.h> using namespace std; double fun(double x) { //return sin(x)+xcos(x);//原函数 return sin(x)/x; } double romberg(double a,double b,double(f)(double),double eps){//上下...

实验原理： 雅克比和高斯的不同之处是，k=1下面那个框框，高斯是a[i][j]*y[j]，雅克比是a[i][j]*x[j]，就是因为高斯用的是更新之后的x，雅克比没有用更新之后的x。 在代码编写过程中，我没有用到动态数组，感觉输入系数矩阵啥的用动态数组更好些。但是，我写了转化对角占优的部分，这部分也是非常重要的。 代码如下： /** 运用雅克比和高斯赛德尔公式求解方程组 a 比较...

#include<stdio.h> #include<math.h> double a[10][10]; double b[10]; int max_i(int i, int n){ int ans = i; int tmp_i; double ans_val=fabs(a[i][i]); for ( tmp_i=i+1;tmp_i&l...