|
odefun=inline('21./x-y./x+x/74','x','y');
[x,y]=ode45(odefun,[1 10],0);[x',y']
ans =
Columns 1 through 7
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
Columns 8 through 14
1.0000 1.0001 1.0001 1.0002 1.0002 1.0003 1.0006
Columns 15 through 21
1.0009 1.0012 1.0015 1.0030 1.0045 1.0060 1.0075
Columns 22 through 28
1.0149 1.0224 1.0299 1.0374 1.0747 1.1121 1.1494
Columns 29 through 35
1.1868 1.3369 1.4871 1.6373 1.7874 2.0124 2.2374
Columns 36 through 42
2.4624 2.6874 2.9124 3.1374 3.3624 3.5874 3.8124
Columns 43 through 49
4.0374 4.2624 4.4874 4.7124 4.9374 5.1624 5.3874
Columns 50 through 56
5.6124 5.8374 6.0624 6.2874 6.5124 6.7374 6.9624
Columns 57 through 63
7.1874 7.4124 7.6374 7.8624 8.0874 8.3124 8.5374
Columns 64 through 70
8.7624 8.9874 9.2124 9.4374 9.6624 9.8874 9.9156
Columns 71 through 77
9.9437 9.9719 10.0000 0 0.0001 0.0001 0.0002
Columns 78 through 84
0.0002 0.0005 0.0007 0.0010 0.0012 0.0025 0.0037
Columns 85 through 91
0.0050 0.0062 0.0125 0.0188 0.0250 0.0313 0.0626
Columns 92 through 98
0.0937 0.1248 0.1558 0.3093 0.4606 0.6097 0.7567
Columns 99 through 105
1.4608 2.1177 2.7319 3.3075 5.2990 6.8860 8.1820
Columns 106 through 112
9.2631 10.5820 11.6351 12.4966 13.2167 13.8263 14.3495
Columns 113 through 119
14.8040 15.2029 15.5560 15.8710 16.1540 16.4099 16.6428
Columns 120 through 126
16.8557 17.0513 17.2319 17.3994 17.5552 17.7008 17.8373
Columns 127 through 133
17.9657 18.0869 18.2015 18.3103 18.4138 18.5125 18.6069
Columns 134 through 140
18.6974 18.7844 18.8680 18.9487 19.0267 19.1023 19.1755
Columns 141 through 146
19.2467 19.3160 19.3245 19.3331 19.3415 19.3500
plot(x,y) |
-
评分
-
1
查看全部评分
-
|