p=[2.79,0.9439] & eresl = 0.15 & y = gauss_smooth(p[0]*rmd_heaviside(x-!pi/2*p[1])*rmd_heaviside(!pi*p[1]*sqrt(2)/2-x)/sqrt(abs(x^2-0.25*(!pi*p[1])^2)),eresl/(x[1]-x[0])/sqrt(8.*alog(2)))
