% bms925, process modeling 1

global Q
global Cin
global V

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

tstart = 0;
tend = 5;

Np_init = 0;

%Np_init = 3;

[t,y] = ode45(@tank2,[tstart tend],Np_init);

%[t,y] = ode15s(@tank2,[tstart tend],Np_init);

plot(t,y)