/* ------------------------------------------
    Catalogue of Spacetimes
    
    Morris-Thorne metric 
        
    Thomas Mueller, Frank Grave
    Visualization Research Center (VISUS)
    University of Stuttgart  
    
    Use with wxmaxima:
       batchload("morristhorne.wxm")
   ------------------------------------------ */

/* load ctensor package */
load(ctensor);

/* define coordinates to use */
ct_coords:[t,l,theta,phi];

/* start with the identity metric */
lg:ident(4);

lg[1,1]:c^2;
lg[2,2]:-1;
lg[3,3]:-(b0^2+l^2);
lg[4,4]:-(b0^2+l^2)*sin(theta)^2;

cmetric();

/* calculate the christoffel symbols of the second kind */
christof(mcs);

/* calculate the riemann tensor */
lriemann(mcs);

/* calculate the ricci tensor */
ricci(mcs);

/* calculate the ricci scalar */
scurvature();

/* calculate the Kretschman scalar */
uriemann(mcs);
rinvariant();
ratsimp(%);