4楼: Originally posted by dbb627 at 2014-05-09 21:51:36
举个例子吧
y''-3y'+2y=x， y‘(0)=1,y’(1)=2

%解析解
yy=dsolve('D2y-3*Dy+2*y=x','Dy(0)=1','Dy(1)=2','x');
x=0:0.01:1;hold on;plot(x,eval(yy),'bo')
%bvp 数值解
f1=@(x,y);
solinit = b ...

6楼: Originally posted by xcenxcen at 2014-05-12 10:16:45
function  = soodesv(q1,q2,q3,alpha,beta,a,b,dx)
N = (b-a)/dx+1;
y = zeros(N,1);
D = zeros(N);
f = zeros(N,1);
for ir=1:N
    x = a+(ir-1)*dx;
    if ir==1
        D(ir,ir) = -q2(x)-2.0/(dx^2) ...