kill(all);

load(ctensor);

ct_coords:[r,theta,phi,t];

depends([%nu,%lambda],[r]);

lg:matrix([-exp(%lambda),0,0,0],[0,-r^2,0,0],[0,0,-r^2*(sin(theta))^2,0],[0,0,0,exp(%nu)]);

cmetric(true);

/* last index is up in Gamma */
christof(mcs);

ricci(true);

leinstein(true);

scurvature();

mu;

%alpha;

%lambda;

%mu;