function [ ] = prob_2_17
% 
% Written by:
% -- 
% John L. Weatherwax                2005-05-07
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

nTs = 10;

ts = linspace( 0, 1, nTs );

[ts,ys] = ode45(@f1,ts,0);
figure; plot( ts, ys, '-og' );

[ts,ys] = ode45(@f2,ts,0);
figure; plot( ts, ys, '-og' );

[ts,ys] = ode45(@f3,ts,0);
figure; plot( ts, ys, '-og' );

return;



function [dydt] = f1(t,y)

dydt = y^2 + 1;

return;



function [dydt] = f2(t,y)

dydt = y^2 + t;

return;



function [dydt] = f3(t,y)

dydt = y^2 + t^2;

return;