real function lossfit(x) implicit none real x common /pawpar/ par(4) lossfit = par(1) + par(2)*exp(-x/par(3)) + par(4)*x end