function xBVP2
% 求解兩點邊值問題的示例
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/04/30 $

clear all
clc
a = 0;
b = 1;

% solution is obtained using an initial guess of y1(x)=0, y2(x)=0 
solinit = bvpinit(linspace(a,b,10),[0 0]);
sol = bvp4c(@ODEfun,@BCfun,solinit);
x = [0:0.1:1];
y = deval(sol,x);

% 結果顯示
DispResult_xBVP

% ------------------------------------------------------------------
function dydx = ODEfun(x,y)
dydx = [y(2); (1+x^2)*y(1)+1];

% ------------------------------------------------------------------
function bc = BCfun(ya,yb)
bc = [ya(1)-1;yb(1)-3];