all 300
cal 1940 1 4

source(noecho) gdpq.rtd

source(noecho) bpfilter.src

set lyy = log(yy)

@bpfilter lyy / lyf

set lyt = lyy - lyf
set yt = exp(lyt)

print / yy yt