| 查看: 606 | 回复: 2 | |||
karry2karry新虫 (初入文坛)
|
[求助]
请问这个问题，用MATLAB代码怎么实现？ 已有1人参与
|
请问这个问题用MATLAB怎么实现 @月只蓝 @beefly 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有247人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
oldoldstone
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2113
- 散金: 466
- 红花: 4
- 帖子: 166
- 在线: 297.2小时
- 虫号: 414842
- 注册: 2007-06-28
- 专业: 生物、医学光子学
【答案】应助回帖
感谢参与，应助指数 +1
|
我去，matlab数值法最大只能算出170的阶乘，符号法最大也只能算出999的阶乘，应该需要进一步近似的。。。 写了一个符号法的，代码放在同一个文件pp.m %pp.m function pp N=0:999; A=zeros(length(N),1); for n=N A(n+1)=fun(n); end plot(N,A) function result=fun(N) syms n m lambda c n=sym(num2str(N)); lambda=sym('3000'); f=lambda^m/factorial(m); %求和函数 S=lambda^n/factorial(n)/symsum(f,m,0,n); %分子的3000我也是用lambda代替了 result=vpa(S,8);%转化为数值 貌似运行需要几分钟，出来几乎是一条直线。 |
2楼2016-12-27 15:17:53
3楼2016-12-27 18:25:06
20