
16: % input exterior field and normal and compute 2nd fundamental forms;
16: la:=-m/(x(2)^2 + x(3)^2)^(1/2) ;
               - m
la := ---------------------
                2       2
       sqrt(x(3)  + x(2) )


17: nu:=-((m*x(2)/(x(2)^2 + x(3)^2))^2)/2;
                 2  2
           - x(2) *m
nu := --------------------
              2       2 2
       2*(x(3)  + x(2) )


18: gp(0,0):=-e^(2*la);
                        -1
gp(0,0) := ----------------------------
                            2       2
             (2*m)/sqrt(x(3)  + x(2) )
            e


19: gp(1,1):=-x(2)^2/gp(0,0);
                           2       2
            (2*m)/sqrt(x(3)  + x(2) )     2
gp(1,1) := e                         *x(2)


20: gp(3,3):=gp(2,2):=e^(2*(nu-la));
                                        2       2
                         (2*m)/sqrt(x(3)  + x(2) )
                        e
gp(3,3) := gp(2,2) := ------------------------------
                             2  2       2       2 2
                        (x(2) *m )/(x(3)  + x(2) )
                       e


21: for i:=0:3 do g(i,i):=gp(i,i) ;

22: for i:=0:3 do let z(i)=x(i) ;

23: off echo ;
24: 24: 24: 24: 24: 24: 24: 
25: 
26: 26: 
27: 27: 
28: 

30: 30: 

31: 
