function odedemo1
%Commented solution of zdot = z
clear
disp(['The time is now:   ' datestr(now)])

%Set up the ODE solver
tspan = [0 1]; %Solve between 0<t<1
zzero = 1;     %initial condition

%The next line solves the ODE
 [t zarray] = ode23(@myrhs,tspan, zzero);

%Output and plot
 plot(t,zarray); axis([tspan 0 max(zarray)])
 xlabel('t'); ylabel('z');
 disp(['  Last z = '  num2str(zarray(end))])  % output
end  % end of function odedemo1




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%The right hand side is defined in this function
function zdot =myrhs(t,z)
 zdot = z;  % This is, essentially, the ODE
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
