gusucode.com > Diabetic retinopathy > S_function.m

    

function y=S_function(x,m,n,p1,p2)


 if (x<m)
     y= n*(1-(1-x/m)^p1);
 else
     y=n+(1-n)*((x-m)/(1-m))^p2;
 end
 return