|
|
[求助]
如何计算这个二重积分(有图有程序)
![如何计算这个二重积分(有图有程序)]()
tz.JPG
下面是我对这个公式中的H1写的一段程序,但是怎么计算也算不出来
a1=1/1000;
a2=11/1000;
b1=5/1000;
b2=30/1000;
h1=0.5/1000;
h2=0.7/1000;
z=0;
x=(1:30)/1000;
%y=(1:30)/1000;
y=1/1000;
syms a b h
R11=sqrt(((b2-b1)*a/(a2-a1)+y).^2+(z-h).^2+(x-a).^2);
R12=sqrt(((b2-b1)*a/(a2-a1)-y).^2+(z-h).^2+(x-a).^2);
HB11=(((b2-b1)*a/(a2-a1)-y)./R11+((b2-b1)*a/(a2-a1)+y)./R12).*(z-h)./((z-h).^2+(x-a).^2);
F=vpa(int(HB11,h,h1,h2));
H=vpa(int(F,a,a1,a2));
然后尝试了下面用泰勒展开
k=taylor(HB11,h);
kk=taylor(k,a);
r=int(kk,a,a1,a2);
rr=int(r,h,h1,h2);
计算出来的结果感觉不对
请大家帮我看看应该怎么改才好呢?谢谢了 |
» 猜你喜欢
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有18人回复
售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
2027博士申请-临床检验诊断学、病原生物学、微生物方向
已经有4人回复
售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急
已经有17人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有22人回复
越来越担心硕士毕业后没有竞争力了
已经有14人回复
论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章
已经有3人回复
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
闲聊
已经有6人回复
|