clear

D = 4*.0254;
h = 3*.0254;

V0 = [pi*(D/2)^2*h];

t0 = 0;
tf = 700;
N = 700;

[t,V] = rk4fixed('onecan',[t0 tf],V0,N);

plot(t,V); title('Volume vs. Time'); xlabel('Time (sec)'); ylabel('Volume (m^3)');