SUBROUTINE MAP(N,F,X) ************************************************************************ * * * Watson's Example 4 * * * ************************************************************************ * REAL*8 F(N),X(N),XT * XT=0 DO 10 I=1,N XT=XT+(X(I)-I+2)*(X(I)-I+2) 10 CONTINUE XT=EXP(XT) * DO 20 I=1,N F(I)=(X(I)-I+2)*XT 20 CONTINUE * RETURN END * * * SUBROUTINE NABMAP(N,T,X) ************************************************************************ * * * Jacobian of the mapping * * * ************************************************************************ * * REAL*8 T(N,N),X(N),XT * XT=0 DO 10 I=1,N XT=XT+(X(I)-I+2)*(X(I)-I+2) 10 CONTINUE XT=EXP(XT) * DO 20 I=1,N DO 25 J=1,N T(I,J)=2*XT*(X(I)-I+2)*(X(J)-J+2) 25 CONTINUE T(I,I)=T(I,I)+XT 20 CONTINUE * RETURN END