function dNp_dt = membrane(t,Np)

global B
global Cout
global Rmax
global V
global Kp5

dNp_dt = B.*Cout - ((Rmax.*Np./V)./(Np./V + Kp5));