def eq1(): r = Symbol("r") e = relativity.Rmn.dd(0,0) e = e.subs(relativity.nu(r), -relativity.lam(r)) print dsolve(e, [relativity.lam(r)])
def eq2(): r = Symbol("r") e = relativity.Rmn.dd(1,1) C = Symbol("CC") e = e.subs(relativity.nu(r), -relativity.lam(r)) print dsolve(e, [relativity.lam(r)])