24小时热门版块排行榜    

查看: 684  |  回复: 0

janecool

铁杆木虫 (著名写手)

[交流] 【求助】怎样用MATLAB编写CG系数【已解决】


请教一下高手们,我下面写的这个程序用MATLAB怎么代入各个角动量的值后不能运算啊?
function racah(j1,m1,j2,m2,j,m)
m1+m2==m

term1=sqrt(((2*j+1)*factorial(j1+j2-j)*factorial(j1+m1)*factorial...
(j2+m2)*factorial(j1-m1)*factorial(j2-m2)*factorial(j+m)...
*factorial(j-m)*factorial(j+j2-j1)*factorial(j+j1-j2))...
/factorial(j+j1+j2+1))

term2=0
n=0
while(j1+j2-j-z>0 &j1-m1-z>0 &j2+m2-z>0 &j-j2+m1+z>0 &j-j1+m2+z>0)
n=n+1
end
for z=0:n
   
    term2=term2+(-1)^z/(factorial(z)*factorial(j1+j2-j-z)*...
        factorial(j1-m1-z)*factorial(j2+m2-z)*...
        factorial(j-j2+m1+z)*factorial(j-j1+m2+z))
end
racah=term1*term2

我自己搞定了

[ Last edited by janecool on 2010-9-14 at 19:11 ]
回复此楼

» 猜你喜欢

相信会有属于你的精彩
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 janecool 的主题更新
信息提示
请填处理意见