function wsp3(x) COMMON/PAWPAR/PAR(4) wsp3=par(1)/(1.0+exp((x-par(2))/par(3)))+par(4)*x end