function dNp_dt = tank2(t,Np)

global Q Cin V

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