showgan@tg
says to
YSITD
S(0,L)(xr/L)²pi*0.5dx