如何解二维常微方程组？

19# 李宏海
应该是可以的，这是我的Matlab的版本信息

1. 
   
2. >> version
   
3. ans =
   
4. 7.11.0.584 (R2010b)

复制代码

21# qibbxxt
貌似楼上的大侠是2010版，好像那个占用好大的！
>> version
ans =
7.0.0.19920 (R14)
The element type "name" must be terminated by the matching end-tag "</name>".
Could not parse the file: d:\matlab7\toolbox\ccslink\ccslink\info.xml
>>
这是我的版本，貌似落伍好多！呵呵

用1stOpt更简单：

1. Variable t=[0,30], x=0.3327, y=0.0497;
   
2. Plot t[x],x,y;
   
3. ODEFunction x'=(-0.0002t+0.0037)/(0.0891685780439393*2*x+0.938619232671706*y);
   
4.             y'=(-0.0004t+0.0089)/(0.0650571169521575*2*y+0.938619232671706*x);

复制代码

23# dingd
一开始我以为1st0pt只用来优化拟合，原来它也可以解方程啊，真是百花齐放，哈哈！

23# dingd
恩，很想用用1stOpt，可惜没有机会用

代码稍加修改如下:

1. LoopConstant a=[0:0.1:3];
   
2. Variable t=[0,30], x=0.3327, y=0.0497;
   
3. Plot x[x],y,x',y';
   
4. ODEFunction x'=(-0.0002t+0.0037*a)/(0.0891685780439393*2*x+0.938619232671706*y);
   
5.             y'=(-0.0004t+0.0089)/(0.0650571169521575*2*y+0.938619232671706*x);

复制代码

26# dingd
楼上的大侠，我想将t的取值范围改为[0,120]是不是直接将其改范围即可，其余代码不变？因为我得到的x[t],y[t]都是关于年龄的生存函数，年龄范围取值在较符合实际的0至120岁区间，

1. Variable t=[0,120], x=0.3327, y=0.0497;
   
2. Plot x[x],y;
   
3. ODEFunction x'=(-0.0002t+0.0037)/(0.0891685780439393*2*x+0.938619232671706*y);
   
4.             y'=(-0.0004t+0.0089)/(0.0650571169521575*2*y+0.938619232671706*x);

复制代码

各位大侠，为什么我用mathematica8.0画不出来t属于[0,120]区间的图形呢！
In[2]=：Plot[Evaluate[{x[t],y[t]}/.%],{t,0,120},PlotRange -> All]
命令对吗，为什么画不出来呢？

16# qibbxxt
大侠，请教：我希望用MATLAB直接产生一组100个x1和x2分别对应的[0，1]上的服从标准正态分布N(0,1)的点，相当于产生一个2×100的矩阵，其中的每一个元素都是服从标准正态分布N(0,1)的点，并且第一行的第i个元素和第二行的第i个元素对应，i=1,2,...,100
求算法！