Matlab code for the McCallum/Nelson model