C      if (je.eq.1) then 



      der(1) = exp( xlog(t-1,1)* b(2,je) + thedlog(t) * b(3,je) 

     £        + xlog(t-1,2) * b(4,je))   

c     £        + (xlog(t-1,1)**2) * b(5,je))

      der(2) = b(1,je)*der(1) * xlog(t-1,1)

      der(3) = b(1,je)*der(1) * thedlog(t)

      der(4) = b(1,je)*der(1) * xlog(t-1,2)

c      der(5) = b(1,je)*der(1) * (xlog(t-1,1)**2)

      

C      elseif (je.eq.2) then



C      der(1) = exp( xlog(t-1,2)* b(2,je) + thedlog(t) * b(3,je)) 

C      der(2) = b(1,je)*der(1) * xlog(t-1,2)

C      der(3) = b(1,je)*der(1) * thedlog(t)

C

C      endif