function i=phi(u)


c=0.00093231545017;  
k= 13.41566535607806;


i=c*(exp(k*u)-1);

