решение интеграла в квадратурном методе Симпсона в 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);