      psi(t,1) =   b(1,1) * exp( zlog(t-1,1) * b(2,1) + 

     &             thetalog(t,1)*b(3,1)) 

