> restart:

Standard Worksheet SSSPF

D-F

Physical quatities are defined in the initialization file

> grts();

Created definition for G(up,dn) 

Created definition for rho

Created definition for Iso

Created definition for p

Created definition for gthetatheta

Created definition for R(dn,dn,up,up) 

Created definition for mass

> qload(df);

Physical parameters

> grcalc(Iso,rho,p,mass):

> grdisplay(_);

Perssure plots (gamma=1)

> pj:=subs(grcomponent(p,[])):

> p1:=subs(R=1,pj):

> plot(p1,xi=0..1,title="Pressure D-F");

Energy density plots

> rhoj:=subs(grcomponent(rho,[])):

> rho1:=subs(R=1,rhoj):

> plot(rho1,xi=0..1,title="Energy Density D-F");

Mass plots

> mj:=subs(grcomponent(mass,[])):

> mj1:=subs(R=1,mj):

> plot(mj1,xi=0..1,title="Mass D-F");

Trapping

Potential impac parameter

> B:=radsimp(sqrt(grcomponent(g(dn,dn),[theta,theta]))/sqrt(-grcomponent(g(dn,dn),[t,t]))):

> Bj:=subs(R=1,B):

> plot(Bj,xi=0..1,title="Trapping D-F");

w - modes

Potential

> V:=1/((Bj^2))*(6+4*Pi*xi^2*R^2*(grcomponent(rho,[])-grcomponent(p,[]))-6*grcomponent(mass,[])/(xi*R)):

> Vj:=subs(R=1,V):

> plot(Vj,xi=0.1..1.0,title="w - modes D-F");

V:=sqrt(dp/dr/drho/dr)

> vs:=sqrt(diff(grcomponent(p,[]),xi)/(diff(grcomponent(rho,[]),xi))):

> vsj:=subs(R=1,vs):

> plot(vsj,xi=0..1,title="V D-F");

>