> restart:
> grtw();
> qload(statica);
> grcalc(G(dn,dn));
> dsolve(grcomponent(G(dn,dn),[t,t])+Lambda*grcomponent(g(dn,dn),[t,t])=0,a(r));
The standard notation for the constant is -2*m
> a(r):=r/(r-r^3*Lambda/3-2*m);
> dsolve(grcomponent(G(dn,dn),[r,r])+Lambda*grcomponent(g(dn,dn),[r,r])=0,b(r));
Scale t so that the constant is 1
> b(r):=-(-r+r^3*Lambda/3+2*m)/r;
> grdisplay(g(dn,dn));
>