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


13: 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) )


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


15: 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)


16: 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


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

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

19: off echo ;
20: 20: 20: 20: 20: 20: 20: 
21: 
22: 22: 
23: 23: 
24: 

26: 26: 

27: 
