public double F2(double T, double dt) { //Relax Pelsser model assumptions double z=f(T, dt) - SIG(T); return Math.Max(0,z); }