追趕法比較特殊,比前面的演算法實現起來都要簡單,唯一要注意的是gamma的座標從1開始
程式分別用3個陣列存放三條對角線
**如下:
///對於係數矩陣是三對角矩陣的方程組可以用追趕法
#include #include #include using namespace std;
const int maxn = 1000;
double a[maxn];
double d[maxn];
double c[maxn];
double b[maxn];
double a[maxn];
double b[maxn];
double x[maxn];
double y[maxn];
int main()
///輸出
/*for(int i = 0;i
追趕法解方程
演算法 詳情見追趕法公式 這個部落格只是用來儲存。package 計算方法大作業 public class zhuiganfa double b double c double b2 儲存貝塔 double y 儲存y double f 獲取所有的貝塔 for i 1 i 3 i else 獲取所有...
MATLAB實現追趕法
function x chase method a,b chase method 追趕法求三對角矩陣的解 a為三對角矩陣的係數,b為等式右端的常數項,返回值x即為最終的解 注 a盡量為方陣,b一定要為列向量 求追趕法所需l及u t a for i 2 size t,1 t i,i 1 t i,i 1...
數值分析之龍貝格求積法
花了一晚上的時間把數值分析裡面定積分求解看懂了,累得早上頭疼發燒。唉革命尚未成功啊!利用richardson外推算法,得到如下的求積方法,其只產生四個序列 其結束迭代準則為 並認為有一定數值分析基礎,不難寫出如下程式 include include include include using nam...