double precision wykl /4.0d0/
if (y.lt.ymin) then
rlornmr1prim=-(ymin-y)**(wykl-1)*sigma**wykl*wykl/
- & ((ymin-y)**wykl+sigma**wykl)
+ & ((ymin-y)**wykl+sigma**wykl)**2
else if (y.gt.ymax) then
rlornmr1prim=(y-ymax)**(wykl-1)*sigma**wykl*wykl/
- & ((y-ymax)**wykl+sigma**wykl)
+ & ((y-ymax)**wykl+sigma**wykl)**2
else
rlornmr1prim=0.0d0
endif