xbvp1.m

Matlab在化学工程中的应用

function xBVP1
% 求解两点边值问题的示例
%
%   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
format long e

a = 0;
b = 1;

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

% 结果显示
DispResult_xBVP

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

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