      SUBROUTINE MAP(N,F,X)
************************************************************************
*                                                                      *
*     Mathiesen's 4 variables Walrasian equilibrium model              *
*                                                                      *
************************************************************************
*
      REAL*8 F(N),X(N),A,B2,B3
*
      A=0.9
      B2=5
      B3=3
*
*      A=0.75
*      B2=1
*      B3=0.5
*     
      F(1)=-X(2)+X(3)+X(4)
      F(2)=X(1)-A*(B2*X(3)+B3*X(4))/X(2)
      F(3)=A*B2-X(1)-(1-A)*B3*X(4)/X(3)
      F(4)=B3-X(1)
*
      RETURN
      END
*
*
*
      SUBROUTINE NABMAP(N,T,X)
************************************************************************
*                                                                      *
*     Jacobien                                                         *
*                                                                      *
************************************************************************
*
*
      REAL*8 T(N,N),X(N),A,B2,B3
*
      A=0.9
      B2=5
      B3=3
*
*      A=0.75
*      B2=1
*      B3=0.5
*
      T(  1,  1)=0
      T(  1,  2)=-1
      T(  1,  3)=1
      T(  1,  4)=1
      T(  2,  1)=1
      T(  2,  2)=A*(B2*X(3)+B3*X(4))/X(2)/X(2)
      T(  2,  3)=-A*B2/X(2)
      T(  2,  4)=-A*B3/X(2)
      T(  3,  1)=-1
      T(  3,  2)=0
      T(  3,  3)=(1-A)*B3*X(4)/X(3)/X(3)
      T(  3,  4)=-(1-A)*B3/X(3)
      T(  4,  1)=-1
      T(  4,  2)=0
      T(  4,  3)=0
      T(  4,  4)=0
*
      RETURN
      END
