plot (r*h–r+d)/h + (h–1)*(m–r)^2/(2*s^2*h^2) at r=0.01, d=0.01, m=0.06, s=0.14, h from 0 to 15