source: MML/trunk/machine/SOLEIL/common/toolbox/chebfun_v2_0501/chebfun/chebtests/bvptest.m @ 17

function pass = bvptest

% Rodrigo Platte Jan 2009
% This routine tests BVP solvers bvp4c and bvp5c

[d,x] = domain(0,4);
y0 = [ x.^0, 0 ];

solinit = bvpinit([0 1 2 3 4],[1 0]); 

% Test bvp4c using default tolerance (RelTol = 1e-3)
y = bvp4c(@twoode,@twobc,y0);         % Chebfun solution
sol = bvp4c(@twoode,@twobc,solinit);  % Matlab's solution
pass = max(max(abs(sol.y' - feval(y,sol.x')))) < 2e-2;

%% Test bvp5c using new tolerance 

% Only test if Matlab ver > 7.4
matlabver = ver('matlab');
if str2double(matlabver.Version) > 7.45
    
    % Set new tolerance:
    opts = odeset('RelTol', 1e-6);
    
    % Test bvp5c
    y = bvp5c(@twoode,@twobc,y0,opts);         % Chebfun solution
    sol = bvp5c(@twoode,@twobc,solinit,opts);  % Matlab's solution
    pass2 = max(max(abs(sol.y' - feval(y,sol.x')))) < 1e-4;
    pass = pass && pass2;
    
end