% bms925, process modeling 1

Q = 0.01;
Cin = 1;
V = 1;

tstart = 0;
tend = 500;

Np_init = 0;

[t,y] = ode45(@tankb,[tstart tend],Np_init,[],Q,Cin,V);

plot(t,y)