辛普森積分法求解積分在MATLAB

Question

我試圖找到一個完整的

b = ((8*WL* Requ)/(pie* WC *Eequ))^.5 

z = quad('exp(sqrt(1-((x.^2)./b^2)))',-b,b) 

我使用內置的命令，此辛普森正交（四）。我得到了正確的結果，但我的問題是，我希望quad命令自動從我的代碼中寫入（b）的公式中自動獲取（b）的變化值。但它沒有那樣做。我必須將（b）的每個變化值都放入函數中。

Answer 1

quad旨在數值積分。在您的發佈代碼中，被積函數被寫爲符號表達式（字符串）。

b = ((8*WL* Requ)/(pie* WC *Eequ))^.5; 

F = @(x)exp(sqrt(1-((x.^2)./b^2)); 
z = quad(F, -b, b);