clear;clc; 

n=12 ; 
vi=150e3;  

C = 0.3e-6 ;
c_1 = 0.02e-6 ;
c_2= 1000e-12 ;   

alpha1_inv=68.2e-6;
alpha1=1/alpha1_inv;

alpha2_inv=0.405e-6;
alpha2=1/alpha2_inv;

inv=(alpha1_inv+alpha2_inv);
inv_s=inv^2;

Cn = c_2+c_1;

R_front=(0.5/c_2)*( inv - sqrt(inv_s - (4*Cn/(alpha2*alpha1*c_1))))

R_tail=(0.5/Cn)*( inv + sqrt(inv_s - (4*Cn/(alpha2*alpha1*c_1))))

tmax = log(alpha2/alpha1)/((alpha2)-(alpha1))
k = R_front*c_2;
vo = 125e3;
Vmax = (n*vo/k)*(1/(alpha2-alpha1))*(exp(-alpha1*tmax)-exp(-alpha2*tmax))

