| 查看: 313 | 回复: 0 | ||
xiaohsu2006木虫 (著名写手)
|
[求助]
含泪跪求1stopt高版本跑个代码 已经写好了
|
|
代码1（550）： ConstStr r1=k1*CA, r2=k2*CB*CD, r3=k3*CB*CD, r4=k4*CB, r5=k5*CB, r6=k6*CB, r7=k7*CB, r8=k8*CC*CD, r9=k9*CC*CD, r10=k10*CD*CF, r10r=k10/3.64*CE*CH, r11=k11*CE*CF, r11r=k11/273.7*CD*CG; Variable t,CE,CF,CG,CH; InitialODEValue t=0, CA=0.730635793, CB=0, CC=1.498246463, CD=3.332429483, CE=0, CF=0, CG=0, CH=0; ODEFunction CA’=-r1; CB’=r1-r2-r3-r4-7*r5-4*r6-11*r7; CC’=-r8-r9; CD’=-r2-2*r3-r8-2*r9-r10+r10r+r11-r11r; CE’=r2+2*r3+6*r4+r8+2*r9+r10-r10r-3*r11+3*r11r; CF’=r2+3*r5+r8-r10+r10r-r11+r11r; CG’=3*r6+r11-r11r; CH’=r3+3*r7+r9+r10-r10r; Data; //t/s CE CF CG CH 0 0 0 0 0 30 0.000149796 0.000760865 0.000217173 0.005495421 60 0.002059783 0.002324358 0.002416821 0.011683547 180 0.006129219 0.002600548 0.004684001 0.014443429 300 0.007712601 0.002375765 0.005132747 0.0146369 600 0.009297344 0.002233589 0.005363171 0.018179198 1200 0.012671286 0.002038349 0.006919546 0.020292707 1800 0.014069376 0.002154183 0.007150669 0.021269902 代码2（700）： ConstStr r1=k1*CA, r2=k2*CB*CD, r3=k3*CB*CD, r4=k4*CB, r5=k5*CB, r6=k6*CB, r7=k7*CB, r8=k8*CC*CD, r9=k9*CC*CD, r10=k10*CD*CF, r10r=k10/1.62*CE*CH, r11=k11*CE*CF, r11r=k11/1.65*CD*CG; Variable t,CE,CF,CG,CH; InitialODEValue t=0, CA=0.730635793, CB=0, CC=1.498246463, CD=3.332429483, CE=0, CF=0, CG=0, CH=0; ODEFunction CA’=-r1; CB’=r1-r2-r3-r4-7*r5-4*r6-11*r7; CC’=-r8-r9; CD’=-r2-2*r3-r8-2*r9-r10+r10r+r11-r11r; CE’=r2+2*r3+6*r4+r8+2*r9+r10-r10r-3*r11+3*r11r; CF’=r2+3*r5+r8-r10+r10r-r11+r11r; CG’=3*r6+r11-r11r; CH’=r3+3*r7+r9+r10-r10r; Data; //t/s CE CF CG CH 0 0 0 0 0 30 0.004375416 0.003086408 0.003315778 0.015026091 60 0.019617178 0.003836818 0.007720506 0.019703998 180 0.038487661 0.003891664 0.011588798 0.031720487 300 0.048176609 0.003933202 0.015502672 0.036733943 600 0.062063433 0.004022325 0.016942726 0.044129082 1200 0.075194185 0.004214514 0.023919896 0.055351964 1800 0.076085076 0.004151429 0.020840457 0.053030061 待求参数范围均在[0,1]内，这个不知道怎么写到代码里，求大神帮忙，跪谢！用matlab和mandona搞了好久得不到好的结果，不知道1stopt咋样 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有103人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有76人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
30