www.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