func8p2.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 18 行

solinit = bvpinit(linspace(0,1,5),[0 1]);
sol = bvp4c(@twoode,@twobc,solinit);
x = linspace(0,1);
y = deval(sol,x);
z=(0:0.1:1)
h=0.1
for i=1:10
    z=-exp(1)/(1+exp(1))*exp(-i*h)-1/(1+exp(1))*exp(h*i)+1
end
plot(x,y(1,:));
%调用matlab求微分方程边值函数bvp4c,
%计算值：0.0379  0.0708  0.0941  0.1082  0.1132  0.1091  0.0959  0.0735  0.0417 
%0.0046  0
%真实值：0.1000  0.0413  0.0730  0.0954  0.1087  0.1132  0.1087  0.0954  0.0730
%0.0413  0   ，在中间部分误差较小
%figure,twode.m,twobc.m enclosed