急求：有2.0以上版本的高手帮忙运行一下程序，1.5版本的跑不了

hongyukelly

小妹的邮箱是691285654@qq.com，不胜感激！

1. Variable x1=[1,], x2=[0,], x3=[,2.981242639847277], t=[0:0.01:1],x4=[,0];
   
2. ODEOptions=[SN=50,A=0,P=2];
   
3. plot x1,x2,x3,x4;
   
4. ODEFunction  x1'= x2+(cos((39^(1/2)*t)/4)-(5*39^(1/2)*sin((39^(1/2)*t)/4))/39)/(2*exp((7*t)/4)) ;
   
5.              x2'=(- (98658077728449830206392979734313631522384319688164975111569408*exp((28769*t)/10000)*(7915705755533312*cos((452220989014229*t)/549755813888000) +
   
6.                   4961373252692701*sin((452220989014229*t)/549755813888000) - 9115229831430144))/(5259268619686907678872730326191907679275300806870751227045469949900953223168*exp((28769*t)/5000) + 434621024932584022745096624004720431321037539799896744197855557638161071865856*exp((28769*t)/10000) - 194282357074795987104433617562072844758930551125810734864269312*cos((452220989014229*t)/549755813888000)^2*exp((28769*t)/5000) + 218232052769054270395520109593353502683087789633521723746811904*cos((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 348902432071189267629930684133713671917909692019215375611963702974917640192000*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 278581280767182523913065996218712995267638133189092027941257216*sin((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 180438718024742897406023448514179223232451285182830782179180544*sin((452220989014229*t)/274877906944000)*exp((28769*t)/5000)
   
7.                   - 218683620606714164366704211765389147722921139007641894058980330532642366685184*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 19680008237368423884372751235849191005687254290127135005647703418985284370432) - 3623/20000)*x1+((1250361226628551173771181388067839950133895921614330182672734825664833583579136*exp((28769*t)/10000) - 823871499857955841258812866845596960936660396359228593615615454976978414731264*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) - 916132880480672965465397890175718976976354179611505219647035943315793321656320*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 113234831396170439448421673169291951727626525353710658535546152997563674918912)/(10518537239373815357745460652383815358550601613741502454090939899801906446336*exp((28769*t)/5000) + 869242049865168045490193248009440862642075079599793488395711115276322143731712*exp((28769*t)/10000)
   
8.                    - 388564714149591974208867235124145689517861102251621469728538624*cos((452220989014229*t)/549755813888000)^2*exp((28769*t)/5000) + 436464105538108540791040219186707005366175579267043447493623808*cos((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 697804864142378535259861368267427343835819384038430751223927405949835280384000*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 557162561534365047826131992437425990535276266378184055882514432*sin((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 360877436049485794812046897028358446464902570365661564358361088*sin((452220989014229*t)/274877906944000)*exp((28769*t)/5000) - 437367241213428328733408423530778295445842278015283788117960661065284733370368*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 39360016474736847768745502471698382011374508580254270011295406837970568740864) + 21381/10000)*x2+x4+(cos((39^(1/2)*t)/4)-(5*39^(1/2)*sin((39^(1/2)*t)/4))/39)/(2*exp((7*t)/4));
   
9.              x3'=((98658077728449830206392979734313631522384319688164975111569408*exp((28769*t)/10000)*(7915705755533312*cos((452220989014229*t)/549755813888000) + 4961373252692701*sin((452220989014229*t)/549755813888000) - 9115229831430144))/(5259268619686907678872730326191907679275300806870751227045469949900953223168*exp((28769*t)/5000) + 434621024932584022745096624004720431321037539799896744197855557638161071865856*exp((28769*t)/10000) - 194282357074795987104433617562072844758930551125810734864269312*cos((452220989014229*t)/549755813888000)^2*exp((28769*t)/5000) + 218232052769054270395520109593353502683087789633521723746811904*cos((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 348902432071189267629930684133713671917909692019215375611963702974917640192000*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 278581280767182523913065996218712995267638133189092027941257216*sin((452220989014229*t)/549755813888000)*exp((28769*t)/5000) -
   
10.              180438718024742897406023448514179223232451285182830782179180544*sin((452220989014229*t)/274877906944000)*exp((28769*t)/5000) - 218683620606714164366704211765389147722921139007641894058980330532642366685184*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 19680008237368423884372751235849191005687254290127135005647703418985284370432) + 3623/20000)*x4+((cos((39^(1/2)*t)/4) - (5*39^(1/2)*sin((39^(1/2)*t)/4))/39)*((98658077728449830206392979734313631522384319688164975111569408*exp((28769*t)/10000)*(7915705755533312*cos((452220989014229*t)/549755813888000) + 4961373252692701*sin((452220989014229*t)/549755813888000) - 9115229831430144))/(5259268619686907678872730326191907679275300806870751227045469949900953223168*exp((28769*t)/5000) + 434621024932584022745096624004720431321037539799896744197855557638161071865856*exp((28769*t)/10000) - 194282357074795987104433617562072844758930551125810734864269312*cos((452220989014229*t)/549755813888000)^2*exp((28769*t)/5000) +
    
11.              218232052769054270395520109593353502683087789633521723746811904*cos((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 348902432071189267629930684133713671917909692019215375611963702974917640192000*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 278581280767182523913065996218712995267638133189092027941257216*sin((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 180438718024742897406023448514179223232451285182830782179180544*sin((452220989014229*t)/274877906944000)*exp((28769*t)/5000) - 218683620606714164366704211765389147722921139007641894058980330532642366685184*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 19680008237368423884372751235849191005687254290127135005647703418985284370432) + 1623/20000))/(2*exp((7*t)/4)) - (5*cos((391^(1/2)*t)/20) - (15*391^(1/2)*sin((391^(1/2)*t)/20))/391)/exp((7*t)/20) - ((cos((39^(1/2)*t)/4) -
    
12.              (5*39^(1/2)*sin((39^(1/2)*t)/4))/39)*((2798683533562686409608157987273346465497443212316367015109278246041893405196288*exp((28769*t)/10000) - 2649348829263438908692229895879980619310565083697151851119680456321100704382976*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) - 765787196925573105246684165624066390216526420665840436015264750081109986902016*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 253453523114048424731818449922972382294833853400082958263898242044695765057536)/(10518537239373815357745460652383815358550601613741502454090939899801906446336*exp((28769*t)/5000) + 869242049865168045490193248009440862642075079599793488395711115276322143731712*exp((28769*t)/10000) - 388564714149591974208867235124145689517861102251621469728538624*cos((452220989014229*t)/549755813888000)^2*exp((28769*t)/5000) + 436464105538108540791040219186707005366175579267043447493623808*cos((452220989014229*t)/549755813888000)*exp((28769*t)/5000) -
    
13.              697804864142378535259861368267427343835819384038430751223927405949835280384000*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 557162561534365047826131992437425990535276266378184055882514432*sin((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 360877436049485794812046897028358446464902570365661564358361088*sin((452220989014229*t)/274877906944000)*exp((28769*t)/5000) - 437367241213428328733408423530778295445842278015283788117960661065284733370368*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 39360016474736847768745502471698382011374508580254270011295406837970568740864) + 3701/20000))/(2*exp((7*t)/4)) ;
    
14.              x4'=-x3+(- (1250361226628551173771181388067839950133895921614330182672734825664833583579136*exp((28769*t)/10000) - 823871499857955841258812866845596960936660396359228593615615454976978414731264*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) - 916132880480672965465397890175718976976354179611505219647035943315793321656320*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 113234831396170439448421673169291951727626525353710658535546152997563674918912)/(10518537239373815357745460652383815358550601613741502454090939899801906446336*exp((28769*t)/5000) + 869242049865168045490193248009440862642075079599793488395711115276322143731712*exp((28769*t)/10000) - 388564714149591974208867235124145689517861102251621469728538624*cos((452220989014229*t)/549755813888000)^2*exp((28769*t)/5000) + 436464105538108540791040219186707005366175579267043447493623808*cos((452220989014229*t)/549755813888000)*exp((28769*t)/5000) -
    
15.              697804864142378535259861368267427343835819384038430751223927405949835280384000*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 557162561534365047826131992437425990535276266378184055882514432*sin((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 360877436049485794812046897028358446464902570365661564358361088*sin((452220989014229*t)/274877906944000)*exp((28769*t)/5000) - 437367241213428328733408423530778295445842278015283788117960661065284733370368*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 39360016474736847768745502471698382011374508580254270011295406837970568740864) - 21381/10000)*x4+               ((cos((39^(1/2)*t)/4) - (5*39^(1/2)*sin((39^(1/2)*t)/4))/39)*((98658077728449830206392979734313631522384319688164975111569408*exp((28769*t)/10000)*(7915705755533312*cos((452220989014229*t)/549755813888000) + 4961373252692701*sin((452220989014229*t)/549755813888000) -
    
16.              9115229831430144))/(5259268619686907678872730326191907679275300806870751227045469949900953223168*exp((28769*t)/5000) + 434621024932584022745096624004720431321037539799896744197855557638161071865856*exp((28769*t)/10000) - 194282357074795987104433617562072844758930551125810734864269312*cos((452220989014229*t)/549755813888000)^2*exp((28769*t)/5000) + 218232052769054270395520109593353502683087789633521723746811904*cos((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 348902432071189267629930684133713671917909692019215375611963702974917640192000*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 278581280767182523913065996218712995267638133189092027941257216*sin((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 180438718024742897406023448514179223232451285182830782179180544*sin((452220989014229*t)/274877906944000)*exp((28769*t)/5000) -
    
17.              218683620606714164366704211765389147722921139007641894058980330532642366685184*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 19680008237368423884372751235849191005687254290127135005647703418985284370432) + 1623/20000))/(2*exp((7*t)/4)) - ((cos((39^(1/2)*t)/4) - (5*39^(1/2)*sin((39^(1/2)*t)/4))/39)*((1250361226628551173771181388067839950133895921614330182672734825664833583579136*exp((28769*t)/10000) - 823871499857955841258812866845596960936660396359228593615615454976978414731264*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) - 916132880480672965465397890175718976976354179611505219647035943315793321656320*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 113234831396170439448421673169291951727626525353710658535546152997563674918912)/(10518537239373815357745460652383815358550601613741502454090939899801906446336*exp((28769*t)/5000) + 869242049865168045490193248009440862642075079599793488395711115276322143731712*exp((28769*t)/10000) -
    
18.              388564714149591974208867235124145689517861102251621469728538624*cos((452220989014229*t)/549755813888000)^2*exp((28769*t)/5000) + 436464105538108540791040219186707005366175579267043447493623808*cos((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 697804864142378535259861368267427343835819384038430751223927405949835280384000*cos((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 557162561534365047826131992437425990535276266378184055882514432*sin((452220989014229*t)/549755813888000)*exp((28769*t)/5000) - 360877436049485794812046897028358446464902570365661564358361088*sin((452220989014229*t)/274877906944000)*exp((28769*t)/5000) - 437367241213428328733408423530778295445842278015283788117960661065284733370368*sin((452220989014229*t)/549755813888000)*exp((28769*t)/10000) + 39360016474736847768745502471698382011374508580254270011295406837970568740864) + 1381/10000))/(2*exp((7*t)/4));
    
19.      

复制代码

dingd

边值估算:
        x3(t=0): 10.3481519515548
        x4(t=0): 40.8299524968696

结果过程:
No.        t           x1(t)                    x2(t)                    x3(t)                    x4(t)                         x1'(t)                   x2'(t)                   x3'(t)                    x4'(t)                   
0        0           1                        0                        10.3481519515548        40.8299524968696              0.5                      42.2182496118202        -42.7657880635508        -199.001575267358       
1        0.1        1.25654238075207        4.33689780737431        7.12117743443723        25.0537921659358              4.69926807880809        46.0736797363238        -23.4330771686729        -122.52886263778         
2        0.2        1.96898706560416        9.52355533550255        5.39498148739676        15.3670400794206              9.7722001412446          59.1917317205102        -12.1492092140406        -74.9786533014292       
3        0.3        3.27045482977999        16.4682038445042        4.52189742324384        9.45642085694031              16.6252217003396        81.0219914898334        -5.98180977190918        -45.608623962181         
4        0.4        5.37910711251874        25.9305985025344        4.09258782183829        5.86585055651611              26.0157946106703        108.812698724304        -3.01465355755542        -27.7035071504243       
5        0.5        8.56983245302351        38.2408090668503        3.85659368791887        3.67603307196848              38.2714576715969        136.927794885473        -1.92957032134668        -17.0367328505606       
6        0.6        13.1212847656019        53.1664796907719        3.67643532774888        2.3096583314855               53.1572735292086        160.812985910308        -1.76135431597678        -10.8637661436189       
7        0.7        19.2754933589434        70.3261638432028        3.49625955219026        1.41578460495884              70.2892977745429        182.586816015169        -1.84629863545884        -7.33667891403959       
8        0.8        27.2579393541992        89.8682365697981        3.31040781665043        0.794869214857879             89.8135818956849        209.862896169023        -1.84315305837268        -5.25040848203988