% bms925, process modeling 3

global B
global Cout
global Rmax
global V
global Kp5


%close all

B = 1;
Cout = 1;
Rmax = 3;
V = 1;
%V = 0.1;

Kp5 = 0.4;

tstart = 0;
tend = 10;

Np_init = 0;


%Kp5_v = [0.1:0.05:1.0];

%figure(1)
%hold on

%for i = 1:length(Kp5_v)
    
%    Kp5 = Kp5_v(i);

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

    %end

%hold off