function dNp_dt = tankb(t,Np,Q,Cin,V)

dNp_dt = Q.*(Cin-(Np./V));