Advanced Search

Question: tensors for GR in maple

Question:tensors for GR in maple

Andrew Carmichael 216 Maple
array
June 20 2008
0
Dear all, I'm trying to solve a difficult metric in maple. I did the following: restart; with(tensor): coords:=[t, r, theta,phi]: g:=array(symmetric,sparse,1..4,1..4): g[1,1]:=1+mu*a(r); g[2,2]:=-(1+mu*b(r)): g[3,3]:=-r^2*(1+mu*b(r)): g[4,4]:=-r^2*(sin(theta))^2*(1+mu*b(r)): g[4,1]:=mu*omega*f(r)*r^2*sin(theta)^2*(1+mu*b(r))/c: metric:=create([-1,-1], eval(g)); tensorsGR(coords,metric,contra_metric,det_met, C1, C2, Rm, Rc, R, G, C): display_allGR (coords,metric,contra_metric, det_met, C1, C2, Rm, Rc, R, G, C); which produces a lengthy but correct result. However, I want to tell maple to neglect all powers of omega greater than 2. Does anybody know how to do this? Drew
Please Wait...