source: MML/trunk/machine/SOLEIL/common/toolbox/chebfun_v2_0501/chebfun/chebtests/ivp_ty_test.m @ 4

function pass = ivp_ty_test

% Rodrigo Platte Jan 2009, modified by Asgeir Birkisson
% This routine tests IVP solver: ode45, ode113, ode15s

% Test ode15s Using default tolerances (RelTol = 1e-3)
[t,y] = ode15s(@vdp1,domain(0,20),[2;0]); % chebfun solution
[tm,ym] = ode15s(@vdp1,[0,20],[2;0]); % Matlab's solution

pass1 = max(max(abs(ym - feval(y,tm)))) < 2e-2;

% Test ode45 Using default tolerances (RelTol = 1e-3)
[t,y] = ode45(@vdp1,domain(0,20),[2;0]); % chebfun solution
[tm,ym] = ode45(@vdp1,[0,20],[2;0]); % Matlab's solution

pass2 = max(max(abs(ym - feval(y,tm)))) < 1e-2;

% Test with different tolerance
%opts = odeset('RelTol', 1e-6);

% Test ode113
%[t,y] = ode113(@vdp1,domain(0,20),[2;0],opts); % chebfun solution
%[tm,ym] = ode113(@vdp1,[0,20],[2;0],opts); % Matlab's solution
%
%pass3 = max(max(abs(ym - feval(y,tm)))) < 1e-5;
pass3 = 1; % commented out to save time

pass = pass1 && pass2 && pass3;