可以用for迴圈或quadv（對向量一次求積分）
quadv做法
f=@(x,t)sin((x-(0:0.01:t)).^2);
y=quadv(@(x)f(x,10),0,pi);
plot(t,y)
%因為原引數有5個，所以第六個場可以傳引數
因為matlab提供的函式積分區域都為方形，所以需要轉化成這樣的區域再積分，這個時候可以轉換變數，也可以轉換函式的表示式（十分巧妙）
x=-1:0.01:1;
y=-1:0.01:1;
[x y]=meshgrid(x,y);
f=@(x,y)sqrt(max(0,1-x.^2-y.^2));
dblquad(f,-1,1,-1,1)
第二個我還不會寫，再想 MATLAB積分方法
一 相關函式 符號積分 int f,v int f,v,a,b 數值積分 trapz x,y 梯形法沿列方向求函式y關於自變數x的積分 cumtrapz x,y 梯形法沿列方向求函式y關於自變數x的累計積分 quad fun,a,b,tol 採用遞推自適應simpson法計算積分 quad1 fun...
Matlab 積分微分
函式 函式功能 trapz y 梯形法求定積分 cumsum y 矩形法求定積分 y quad fun,a,b,e 求被積函式fun在 a,b 上的定積分，誤差為e，預設為e 6。演算法是變步長演算法 y quadl fun,a,b,e 與quad類似，但是速度，精度更高 q dblquad fun...
Matlab數值積分
function gc,ak,xk guasschebyshev fun,n 高斯 勒讓德數值積分 引數說明 fun 積分表示式，可以是函式控制代碼 inline函式 匿名函式 字串表示式，但是必須可以接受向量輸入 a,b 積分上下限只能 1 1 故忽略 n 積分階數，可以任意正整數，但是不建議設定...