5: 

6: off exp;

7: 
7: % a programme to calculate the surface energy tensor or the Curzon field in
7: % quasi-cylindrical coords. The surface is a cylinder x(0)=time, x(1)=azimuthal angle,
7: % x(2)=radial, x(3)=axial.
7: 
7: operator q, x, y, z;

8: for all i,j such that i neq j
8:    let {df (x (i), x(j)) => 0,
8:       df (y(i), y(j)) => 0,
8:       df (q(i), q(j)) => 0};

9: 
9: array g(3,3), gp(3,3), h(3,3), cs1(3,3,3), cs2(3,3,3),
9:    n(3), twop(2,2), two(2,2) ;

10: 


11: 
