js8call/lib/peakup.f90
2018-02-08 21:28:33 -05:00

9 lines
119 B
Fortran

subroutine peakup(ym,y0,yp,dx)
b=(yp-ym)/2.0
c=(yp+ym-2.0*y0)/2.0
dx=-b/(2.0*c)
return
end subroutine peakup