Question

用matlab 用数值的方法求积分 题目如下

clear;syms r R z y theta ;theta=0:0.01*pi:2*pi;q=-r*cos(theta)*(4*pi*0.0000001./(2*pi))./sqrt( (sqrt((r.^(2)+y.^(2)-2*r*y*sin(theta)))+R).^(2)+z.^(2))*(ellipticK((4*R*(r.^2 - 2*sin(theta)*r*y + y.^2).^(1/2))/((R + (r.^2 - 2*sin(theta)*r*y + y.^2).^(1/2)).^2 + z.^2))+(ellipticE((4*R*(r.^2 - 2*sin(theta)*r*y + y.^2).^(1/2))/((R + (r.^2 - 2*sin(theta)*r*y + y.^2).^(1/2)).^2 + z.^2)))*((r.^2+y.^2-2*r*y*sin(theta)-R.^2+z.^2)./((sqrt(r.^2+y.^2-2*r*y*sin(theta))-R).^2+z.^2)))
;trapz(theta,q)
结果matlab报错了
Undefined function or method 'ellipticK' for input arguments of type
'sym'.
貌似matlab中是没有ellipticK ellipticE 函数的 得调用maple中大 但我具体该怎么做 求大神赐教 很急 麻烦大家了

Answer 1

用matlab（R2010a）中有特殊函数mfun()可以求。

MFUN Numeric evaluation of a special function.

例如，mfun('EllipticE',k)可以得到第二类完全椭圆积分数值解。
k=sin60°=0.8660254
mfun('EllipticE',0.8660254)
ans =
1.2111

追问

谢谢 你不过我的EllipticE 里是表达式 不是一个具体的数值 那该怎么处理呢

追答

那只能用逼近法去处理喽。