在matlab中求解Simpson求积分积分

Question

我正试图找到一个完整的

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

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

我使用了Simpson正交（quad）内置命令。我得到了正确的结果，但我的问题是，我希望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);