24小时热门版块排行榜    

查看: 238  |  回复: 1

JinkuiZhao

新虫 (初入文坛)

[求助] maple,求大神帮我看看这个程序,为何没结果 已有1人参与

k := 1;
R := 1;
alpha := 1;
H := 1;
w[0] := 1;
eq[m] := A[m]*(HankelH1(m-1, k*R)-HankelH1(m+1, k*R))+sum(A[n](-exp(-(1/2)*i(m+n)*pi)*HankelH1(-m-n, 2*k*H)*BesselJ(1-m, k*R)+exp(-(1/2)*i(m+n)*pi)*HankelH1(-m-n, 2*k*H)*BesselJ(-1-m, k*R)), n = -2 .. 2)
= -i^m*w[0]*(BesselJ(m-1, k*R)-BesselJ(m+1, k*R))*(exp(-i*H*k*sin(alpha)-i*m*alpha)+exp(i*H*k*sin(alpha)+i*m*alpha));
eq[m1] := subs(m = 1, eq[m]);
eq[m2] := subs(m = 2, eq[m]);
eq[m0] := subs(m = 0, eq[m]);
eq[m-1] := subs(m = -1, eq[m]);
eq[m-2] := subs(m = -2, eq[m]);
solve({eq[m-2], eq[m-1], eq[m0], eq[m1], eq[m2]}, {A[-2], A[-1], A[0], A[1], A[2]});
回复此楼

» 猜你喜欢

微纳米力学中的数学问题。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
如果 A[-2],...,A[2]都是一次线性的,可以尝试改写成系数矩阵和向量形式

用 LinearAlgebra:-LinearSolve() 求求看
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
2楼2016-10-29 14:48:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 JinkuiZhao 的主题更新
信息提示
请填处理意见