6楼: Originally posted by 月只蓝 at 2013-05-15 09:24:36
你的方程少个边界条件吧，是不是少个"dy/dx(x=某数值)= ”的边界条件？ 不然按你目前的边界条件，都是x=0处的初始值，这个微分方程已经不是边值问题了，而是初值问题了。在现有条件下，解该问题用ode可解，程序 ...

8楼: Originally posted by somomo91 at 2013-05-15 18:28:25
不必用1stopt， MATLAB 完全可以做到
这里是程序 ：
===================================
% main program
close all;  clear all;  clc
global a
a = 0.1;
t0    = 0;   tend=5000;
tspan = t0:0.5:tend;
...

11楼: Originally posted by musejianglin at 2013-05-16 21:51:04
谢谢，高手。能再帮我看看么？能不能弄出来？

附件.png
...

13楼: Originally posted by 月只蓝 at 2013-05-17 08:58:08
边界条件还是有问题吧，边界条件顾名思义是将解确定在一定的区间之内，你的边界条件只有右端，没有左端，怎么确定解呢，这样子说不定有无穷多个解。
能不能给出比如"y(r=b)=" 或者"dy/dr(r=b)=&qu ...

14楼: Originally posted by musejianglin at 2013-05-17 10:34:04
哦，这是个球，前面的条件是在r=a处即球的表面处，而在0<r<a范围内都是0.这样可以了么？...

14楼: Originally posted by musejianglin at 2013-05-17 10:34:04
哦，这是个球，前面的条件是在r=a处即球的表面处，而在0<r<a范围内都是0.这样可以了么？...

15楼: Originally posted by 月只蓝 at 2013-05-17 10:50:06
是y(r=0)=0，还是dy/dr(r=0)=0呢？
不可能说0<r<a范围内都是0，如果这样的话，这个区间内原函数的数值全部确定了，都不需要解方程了。...

16楼: Originally posted by 月只蓝 at 2013-05-17 10:51:16
你的原始模型的方程一定要建对，不然解出来的结果都是没有意义的。...