//數值計算實驗 數值積分
#include #include #include #include using namespace std;
//被積函式
double f(double x)
//被積函式 給gsl用的
double g(double x, void * params)
//原函式 用於除錯演算法
double f(double x)
//輸出被積函式的精確解
double jinque(const double a, const double b)
//梯形法 求函式在[a,b]上的定積分,積分區間分為n部分
double tixing(const double & a, const double & b, const int & n)
return sum;
}//拋物線法
double paowuxian(const double & a, const double & b, const int & n)
return sum;
}//柯特斯公式
double cotes(const double & a, const double & b, const int & n)
return sum;
}//gsl解法,參考gsl文件
double gslintegration(double & a, double & b)
int main()
{ double a, b;
int n;
cout<<"請輸入積分區間:"<>a;
cout<<"b = ";
cin>>b;
cout<<"請輸入分割被積區間的數量:";
cin>>n;
if (a > b || n <= 1)
{cout<<"輸入錯誤!"<
window數值分析庫gsl安裝
成功的前提條件,選對編譯系統,在win7下可行 1 2 安裝mingw,目的是使用make.exe,見 1 鏈結,如 1 中test.c可編譯,即安裝成功 安裝msys,為了在window下使用configure,make,建議,將這2個程式放在同1目錄 安裝了msys,在其安裝路徑下,找到msys...
Matlab數值積分
function gc,ak,xk guasschebyshev fun,n 高斯 勒讓德數值積分 引數說明 fun 積分表示式,可以是函式控制代碼 inline函式 匿名函式 字串表示式,但是必須可以接受向量輸入 a,b 積分上下限只能 1 1 故忽略 n 積分階數,可以任意正整數,但是不建議設定...
matlab 數值積分!
1.梯形格式 2.其他格式 matlab提供了quad quadl 等函式命令計算一元函式的數值積分,其中quad 採用遍歷的自適應simpson法計算函式的數值積分,quadl採用遍歷的自適應lobatto法計算函式的數值積分。其呼叫格式一般為 q,fcnt quad function,a,b 其...