% This file calls the new version of the ODE solver, ode45 (4th 
% and 5th Runge-Kutta), to solve the dimensionless basic tank 
% model in section 3.2.
%
[T, Z] = ode45n('dfun32',[ 0 2.0], 1.0);

%
% F is the dimensionless inflow rate.
%
F = 0.1*sin(10*T)+1;

plot(T,Z,'-', T,F,'--');
title('Dimensionless Tank Response (tank32)');
xlabel('Dimensionless time');
text(0.25, 1.07, 'Dimensionless inflow rate');
text(0.5, 0.99, 'Dimensionless tank height');
