public double f(double t, double dt) { if (t == 0) return ZR(t); else return (ZR(t) * t - ZR(t - dt) * (t - dt)) / dt; }