How to vectorize integral function in a code?

16 次查看（过去 30 天）
Priya Khot
Priya Khot 2017-2-21
syms r syms p fun1=@(r)(r.^((1-bta)/alpha)/pi*alpha).*exp(-r.^(1/alpha)).*(r*sin(pi*(1-bta))-z*sin(pi*(1-bta+alpha)))/((r.^2-2*r*z*cos(pi*alpha)+z.^2)); % w=r*(eps^(1/alpha))*sin(r/alpha)+r*(1+(1-bta)/alpha);
fun2=@(p)((eps^(1+(1-bta)/alpha))/(2*pi*alpha))*((exp((eps^(1/alpha))*cos(p/alpha)).*... (cos(p*(eps^(1/alpha))*sin(p/alpha)+p*(1+(1-bta)/alpha))+1i*sin(p*(eps^(1/alpha))*sin(p/alpha)+p*(1+(1-bta)/alpha)))))/(eps*exp(1i*p)-z); if (bta<=1) sm=integral(fun1,0,ro,'ArrayValued',true); else sm=(integral(fun1,1,ro,'ArrayValued',true)+integral(fun2,(-alpha*pi),(alpha*pi),'ArrayValued',true)); end

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by