| 查看: 221 | 回复: 2 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】请教高手帮忙解决问题所在，重金酬谢！
|
|||
|
E=linspace(4000,200000,1600); S=int('exp(-E/(8.314*T)','T',507,512); S 我本意是想求活化能E在4000-200000范围内，所对应的各个不同的积分值S，但是却出现问题，请问如何解决？非常感谢！ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有70人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12620.7
- 散金: 200
- 帖子: 642
- 在线: 680.6小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kuhailangyu(金币+2,VIP+0):欢迎常来 9-9 11:04
spring180(金币+25,VIP+0): 10-16 09:05
kuhailangyu(金币+2,VIP+0):欢迎常来 9-9 11:04
spring180(金币+25,VIP+0): 10-16 09:05
|
下面的程序计算速度快一些： Ek=linspace(4000,200000,1600); F=inline('exp(-E./(8.314*T))','T','E'); for k=1:length(Ek); Sd(k)=quad(F,507,512,1e-6,[],Ek(k)); k end plot(Ek,Sd) |
3楼2009-09-09 10:40:31
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12620.7
- 散金: 200
- 帖子: 642
- 在线: 680.6小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
2楼2009-09-09 10:31:29
