辛普森积分法求解积分在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);