matlab ode15s,利用ode45和ode15s求解偏微分方程组相关

odefun=@(t,X)[- 375.47*X(1)^50 + 260.71*X(2)^50 + 3.0217*X(3)^50 + 22.624*X(4)^50 + 10.999*X(5)^50 + 8.6579*X(6)^50 + 6.5514*X(7)^50 + 5.2861*X(8)^50 + 4.3805*X(9)^50 + 3.7214*X(10)^50 + 3.2213*X(11)^50 + 2.8316*X(12)^50 + 2.5208*X(13)^50 + 2.2683*X(14)^50 + 2.0599*X(15)^50 + 1.8855*X(16)^50 + 1.7378*X(17)^50 + 1.6116*X(18)^50 + 1.5027*X(19)^50 + 1.4081*X(20)^50 + 1.3253*X(21)^50 + 1.2524*X(22)^50 + 1.1879*X(23)^50 + 1.1305*X(24)^50 + 1.0794*X(25)^50 + 1.0336*X(26)^50 + 0.99241*X(27)^50 + 0.95533*X(28)^50 + 0.92187*X(29)^50 + 0.89161*X(30)^50 + 0.86422*X(31)^50 + 0.83939*X(32)^50 + 0.81689*X(33)^50 + 0.79648*X(34)^50 + 0.77799*X(35)^50 + 0.76125*X(36)^50 + 0.74613*X(37)^50 + 0.7325*X(38)^50 + 0.72026*X(39)^50 + 0.70931*X(40)^50 + 0.69958*X(41)^50 + 0.691*X(42)^50 + 0.6835*X(43)^50 + 0.67705*X(44)^50 + 0.67159*X(45)^50 + 0.66709*X(46)^50 + 0.66353*X(47)^50 + 0.66087*X(48)^50 + 0.65911*X(49)^50 + 0.65823*X(50)^50 - 32004.0*X(51)

212.69*X(1)^50 - 453.96*X(2)^50 + 210.99*X(3)^50 - 10.34*X(4)^50 + 12.235*X(5)^50 + 4.0773*X(6)^50 + 3.3599*X(7)^50 + 2.3579*X(8)^50 + 1.8384*X(9)^50 + 1.4729*X(10)^50 + 1.2184*X(11)^50 + 1.0311*X(12)^50 + 0.8892*X(13)^50 + 0.77871*X(14)^50 + 0.69082*X(15)^50 + 0.6196*X(16)^50 + 0.561*X(17)^50 + 0.51212*X(18)^50 + 0.47089*X(19)^50 + 0.43576*X(20)^50 + 0.40556*X(21)^50 + 0.3794*X(22)^50 + 0.35659*X(23)^50 + 0.33658*X(24)^50 + 0.31893*X(25)^50 + 0.3033*X(26)^50 + 0.2894*X(27)^50 + 0.27699*X(28)^50 + 0.26589*X(29)^50 + 0.25593*X(30)^50 + 0.24697*X(31)^50 + 0.2389*X(32)^50 + 0.23163*X(33)^50 + 0.22508*X(34)^50 + 0.21917*X(35)^50 + 0.21385*X(36)^50 + 0.20906*X(37)^50 + 0.20476*X(38)^50 + 0.2009*X(39)^50 + 0.19747*X(40)^50 + 0.19443*X(41)^50 + 0.19175*X(42)^50 + 0.18942*X(43)^50 + 0.18742*X(44)^50 + 0.18572*X(45)^50 + 0.18433*X(46)^50 + 0.18323*X(47)^50 + 0.18241*X(48)^50 + 0.18186*X(49)^50 + 0.18159*X(50)^50 - 8364.2*X(51)

2.0816*X(1)^50 + 178.17*X(2)^50 - 383.35*X(3)^50 + 178.04*X(4)^50 - 8.7898*X(5)^50 + 10.286*X(6)^50 + 3.409*X(7)^50 + 2.8101*X(8)^50 + 1.9689*X(9)^50 + 1.5339*X(10)^50 + 1.228*X(11)^50 + 1.0152*X(12)^50 + 0.85876*X(13)^50 + 0.74032*X(14)^50 + 0.64818*X(15)^50 + 0.57494*X(16)^50 + 0.51563*X(17)^50 + 0.46686*X(18)^50 + 0.42621*X(19)^50 + 0.39194*X(20)^50 + 0.36276*X(21)^50 + 0.3377*X(22)^50 + 0.316*X(23)^50 + 0.29709*X(24)^50 + 0.28052*X(25)^50 + 0.26592*X(26)^50 + 0.253*X(27)^50 + 0.24151*X(28)^50 + 0.23128*X(29)^50 + 0.22213*X(30)^50 + 0.21393*X(31)^50 + 0.20657*X(32)^50 + 0.19995*X(33)^50 + 0.194*X(34)^50 + 0.18865*X(35)^50 + 0.18384*X(36)^50 + 0.17952*X(37)^50 + 0.17565*X(38)^50 + 0.17218*X(39)^50 + 0.1691*X(40)^50 + 0.16637*X(41)^50 + 0.16398*X(42)^50 + 0.16189*X(43)^50 + 0.1601*X(44)^50 + 0.15858*X(45)^50 + 0.15734*X(46)^50 + 0.15636*X(47)^50 + 0.15562*X(48)^50 + 0.15514*X(49)^50 + 0.1549*X(50)^50 - 6987.1*X(51)

13.487*X(1)^50 - 7.556*X(2)^50 + 154.08*X(3)^50 - 332.56*X(4)^50 + 153.68*X(5)^50 - 7.9155*X(6)^50 + 8.6678*X(7)^50 + 2.7621*X(8)^50 + 2.2763*X(9)^50 + 1.572*X(10)^50 + 1.2134*X(11)^50 + 0.9627*X(12)^50 + 0.78975*X(13)^50 + 0.66343*X(14)^50 + 0.56842*X(15)^50 + 0.49496*X(16)^50 + 0.43689*X(17)^50 + 0.39012*X(18)^50 + 0.35184*X(19)^50 + 0.32009*X(20)^50 + 0.29343*X(21)^50 + 0.27083*X(22)^50 + 0.25149*X(23)^50 + 0.2348*X(24)^50 + 0.22032*X(25)^50 + 0.20766*X(26)^50 + 0.19655*X(27)^50 + 0.18674*X(28)^50 + 0.17806*X(29)^50 + 0.17034*X(30)^50 + 0.16346*X(31)^50 + 0.15731*X(32)^50 + 0.15182*X(33)^50 + 0.14689*X(34)^50 + 0.14248*X(35)^50 + 0.13853*X(36)^50 + 0.13499*X(37)^50 + 0.13182*X(38)^50 + 0.129*X(39)^50 + 0.1265*X(40)^50 + 0.12429*X(41)^50 + 0.12235*X(42)^50 + 0.12066*X(43)^50 + 0.11921*X(44)^50 + 0.11799*X(45)^50 + 0.11699*X(46)^50 + 0.1162*X(47)^50 + 0.11561*X(48)^50 + 0.11522*X(49)^50 + 0.11503*X(50)^50 - 4983.0*X(51)

5.7794*X(1)^50 + 7.8799*X(2)^50 - 6.7041*X(3)^50 + 135.45*X(4)^50 - 293.27*X(5)^50 + 135.32*X(6)^50 - 7.0757*X(7)^50 + 7.5594*X(8)^50 + 2.3683*X(9)^50 + 1.9503*X(10)^50 + 1.3372*X(11)^50 + 1.0272*X(12)^50 + 0.81099*X(13)^50 + 0.66244*X(14)^50 + 0.55432*X(15)^50 + 0.47328*X(16)^50 + 0.41082*X(17)^50 + 0.36159*X(18)^50 + 0.32206*X(19)^50 + 0.2898*X(20)^50 + 0.26311*X(21)^50 + 0.24076*X(22)^50 + 0.22185*X(23)^50 + 0.20571*X(24)^50 + 0.19183*X(25)^50 + 0.17979*X(26)^50 + 0.16931*X(27)^50 + 0.16012*X(28)^50 + 0.15203*X(29)^50 + 0.14488*X(30)^50 + 0.13854*X(31)^50 + 0.13291*X(32)^50 + 0.12789*X(33)^50 + 0.12341*X(34)^50 + 0.11942*X(35)^50 + 0.11584*X(36)^50 + 0.11266*X(37)^50 + 0.10982*X(38)^50 + 0.10729*X(39)^50 + 0.10505*X(40)^50 + 0.10308*X(41)^50 + 0.10135*X(42)^50 + 0.09985*X(43)^50 + 0.098565*X(44)^50 + 0.097484*X(45)^50 + 0.096597*X(46)^50 + 0.095896*X(47)^50 + 0.095375*X(48)^50 + 0.095031*X(49)^50 + 0.094859*X(50)^50 - 3960.7*X(51)

4.0666*X(1)^50 + 2.3475*X(2)^50 + 7.013*X(3)^50 - 6.2365*X(4)^50 + 120.97*X(5)^50 - 262.26*X(6)^50 + 120.9*X(7)^50 - 6.3814*X(8)^50 + 6.7114*X(9)^50 + 2.078*X(10)^50 + 1.7097*X(11)^50 + 1.1658*X(12)^50 + 0.89204*X(13)^50 + 0.70149*X(14)^50 + 0.57096*X(15)^50 + 0.4762*X(16)^50 + 0.40536*X(17)^50 + 0.3509*X(18)^50 + 0.30809*X(19)^50 + 0.27379*X(20)^50 + 0.24586*X(21)^50 + 0.22281*X(22)^50 + 0.20355*X(23)^50 + 0.1873*X(24)^50 + 0.17345*X(25)^50 + 0.16156*X(26)^50 + 0.15127*X(27)^50 + 0.14233*X(28)^50 + 0.13451*X(29)^50 + 0.12765*X(30)^50 + 0.12159*X(31)^50 + 0.11624*X(32)^50 + 0.11149*X(33)^50 + 0.10727*X(34)^50 + 0.10351*X(35)^50 + 0.10017*X(36)^50 + 0.097202*X(37)^50 + 0.094561*X(38)^50 + 0.092218*X(39)^50 + 0.090147*X(40)^50 + 0.088325*X(41)^50 + 0.086732*X(42)^50 + 0.085352*X(43)^50 + 0.084172*X(44)^50 + 0.08318*X(45)^50 + 0.082367*X(46)^50 + 0.081725*X(47)^50 + 0.081249*X(48)^50 + 0.080933*X(49)^50 + 0.080776*X(50)^50 - 3243.0*X(51)

2.7821*X(1)^50 + 1.749*X(2)^50 + 2.1014*X(3)^50 + 6.1743*X(4)^50 - 5.7188*X(5)^50 + 109.3*X(6)^50 - 237.16*X(7)^50 + 109.26*X(8)^50 - 5.8009*X(9)^50 + 6.0412*X(10)^50 + 1.8558*X(11)^50 + 1.5257*X(12)^50 + 1.0363*X(13)^50 + 0.79062*X(14)^50 + 0.61989*X(15)^50 + 0.50317*X(16)^50 + 0.4186*X(17)^50 + 0.3555*X(18)^50 + 0.30708*X(19)^50 + 0.26908*X(20)^50 + 0.2387*X(21)^50 + 0.21401*X(22)^50 + 0.19367*X(23)^50 + 0.1767*X(24)^50 + 0.1624*X(25)^50 + 0.15024*X(26)^50 + 0.13982*X(27)^50 + 0.13083*X(28)^50 + 0.12302*X(29)^50 + 0.11621*X(30)^50 + 0.11024*X(31)^50 + 0.10498*X(32)^50 + 0.10035*X(33)^50 + 0.096246*X(34)^50 + 0.092613*X(35)^50 + 0.089391*X(36)^50 + 0.086534*X(37)^50 + 0.084002*X(38)^50 + 0.081763*X(39)^50 + 0.079788*X(40)^50 + 0.078054*X(41)^50 + 0.076542*X(42)^50 + 0.075234*X(43)^50 + 0.074117*X(44)^50 + 0.073179*X(45)^50 + 0.072411*X(46)^50 + 0.071805*X(47)^50 + 0.071356*X(48)^50 + 0.071059*X(49)^50 + 0.070911*X(50)^50 - 2735.0*X(51)

2.0483*X(1)^50 + 1.12*X(2)^50 + 1.5807*X(3)^50 + 1.7954*X(4)^50 + 5.5751*X(5)^50 - 5.2646*X(6)^50 + 99.704*X(7)^50 - 216.43*X(8)^50 + 99.681*X(9)^50 - 5.3128*X(10)^50 + 5.4958*X(11)^50 + 1.6788*X(12)^50 + 1.3792*X(13)^50 + 0.93404*X(14)^50 + 0.71101*X(15)^50 + 0.55618*X(16)^50 + 0.45048*X(17)^50 + 0.37401*X(18)^50 + 0.31703*X(19)^50 + 0.27338*X(20)^50 + 0.23918*X(21)^50 + 0.21187*X(22)^50 + 0.1897*X(23)^50 + 0.17147*X(24)^50 + 0.15628*X(25)^50 + 0.14351*X(26)^50 + 0.13266*X(27)^50 + 0.12338*X(28)^50 + 0.11539*X(29)^50 + 0.10846*X(30)^50 + 0.10242*X(31)^50 + 0.09714*X(32)^50 + 0.092504*X(33)^50 + 0.088422*X(34)^50 + 0.084822*X(35)^50 + 0.081642*X(36)^50 + 0.078831*X(37)^50 + 0.076349*X(38)^50 + 0.074159*X(39)^50 + 0.072233*X(40)^50 + 0.070546*X(41)^50 + 0.069077*X(42)^50 + 0.067809*X(43)^50 + 0.066728*X(44)^50 + 0.065822*X(45)^50 + 0.065081*X(46)^50 + 0.064496*X(47)^50 + 0.064063*X(48)^50 + 0.063777*X(49)^50 + 0.063634*X(50)^50 - 2354.8*X(51)

1.6975*X(1)^50 + 0.87323*X(2)^50 + 1.1075*X(3)^50 + 1.4796*X(4)^50 + 1.7467*X(5)^50 + 5.5369*X(6)^50 - 5.2934*X(7)^50 + 99.681*X(8)^50 - 216.45*X(9)^50 + 99.665*X(10)^50 - 5.3266*X(11)^50 + 5.4838*X(12)^50 + 1.6681*X(13)^50 + 1.3698*X(14)^50 + 0.92552*X(15)^50 + 0.70328*X(16)^50 + 0.54912*X(17)^50 + 0.444*X(18)^50 + 0.36803*X(19)^50 + 0.3115*X(20)^50 + 0.26823*X(21)^50 + 0.23437*X(22)^50 + 0.20737*X(23)^50 + 0.18549*X(24)^50 + 0.1675*X(25)^50 + 0.15255*X(26)^50 + 0.13998*X(27)^50 + 0.12933*X(28)^50 + 0.12023*X(29)^50 + 0.1124*X(30)^50 + 0.10562*X(31)^50 + 0.099734*X(32)^50 + 0.094591*X(33)^50 + 0.090086*X(34)^50 + 0.08613*X(35)^50 + 0.08265*X(36)^50 + 0.079586*X(37)^50 + 0.076889*X(38)^50 + 0.074517*X(39)^50 + 0.072437*X(40)^50 + 0.070619*X(41)^50 + 0.069039*X(42)^50 + 0.067678*X(43)^50 + 0.06652*X(44)^50 + 0.06555*X(45)^50 + 0.064758*X(46)^50 + 0.064134*X(47)^50 + 0.063672*X(48)^50 + 0.063366*X(49)^50 + 0.063214*X(50)^50 - 2241.5*X(51)

1.4421*X(1)^50 + 0.69962*X(2)^50 + 0.86279*X(3)^50 + 1.0218*X(4)^50 + 1.4383*X(5)^50 + 1.7143*X(6)^50 + 5.5126*X(7)^50 - 5.3128*X(8)^50 + 99.665*X(9)^50 - 216.47*X(10)^50 + 99.654*X(11)^50 - 5.3367*X(12)^50 + 5.4749*X(13)^50 + 1.6602*X(14)^50 + 1.3626*X(15)^50 + 0.91905*X(16)^50 + 0.69739*X(17)^50 + 0.54372*X(18)^50 + 0.43904*X(19)^50 + 0.36344*X(20)^50 + 0.30723*X(21)^50 + 0.26426*X(22)^50 + 0.23067*X(23)^50 + 0.20391*X(24)^50 + 0.18224*X(25)^50 + 0.16446*X(26)^50 + 0.14968*X(27)^50 + 0.13729*X(28)^50 + 0.12679*X(29)^50 + 0.11784*X(30)^50 + 0.11014*X(31)^50 + 0.1035*X(32)^50 + 0.097734*X(33)^50 + 0.09271*X(34)^50 + 0.08832*X(35)^50 + 0.084474*X(36)^50 + 0.081102*X(37)^50 + 0.078144*X(38)^50 + 0.075551*X(39)^50 + 0.073283*X(40)^50 + 0.071307*X(41)^50 + 0.069593*X(42)^50 + 0.06812*X(43)^50 + 0.066869*X(44)^50 + 0.065822*X(45)^50 + 0.064968*X(46)^50 + 0.064296*X(47)^50 + 0.063799*X(48)^50 + 0.063471*X(49)^50 + 0.063307*X(50)^50 - 2148.2*X(51)

1.2483*X(1)^50 + 0.57875*X(2)^50 + 0.69072*X(3)^50 + 0.78873*X(4)^50 +

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值