y_0 = 0.1;
y_list = linspace (0, y_0, 1000);
fy_list = y_list。^ (1.5);
f_inter = @ (y) interp1 (y_list fy_list y);
L = 10;
n = 100;
xmesh = linspace (0, L, n);
猜= @ (x) 1的军医。* (1,1);
nsolve = 100;
为isolve = 1: nsolve
y1end = y_0 * isolve / nsolve;
solinit = bvpinit (xmesh,猜);
索尔= bvp4c (@ (x, y) bvpfcn (x, y, L, f_inter), @(是的,yb) bcfcn(是的,yb, y1end) solinit);
der (isolve) = sol.y (2);
xmesh = sol.x;
猜= @ (x) [interp1 (sol.x, sol.y (1:), x); interp1 (sol.x, sol.y (2:), x));
结束
情节(y_0 * (1: nsolve) / nsolve der)
% % % % %
函数就要= bvpfcn (x, y, L, f_inter)%方程来解决
fy_inter = f_inter (y (1));
就要= 0 (2,1);
它们(1)= y (2);% % (1)= y, y (2) = y ' .....
如果abs (x) < L / 2
它们(2)= fy_inter;
其他的
它们(2)= 0;
结束;
结束
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
函数res = bcfcn(是的,yb, y1end)%边界条件
res =[丫(1);
yb (1) -y1end];% y_0 = 1 e 1
结束