SHARE
TWEET

Untitled

a guest Jun 18th, 2019 50 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. rmax = 100, g = 3, n = 10^10, B = 160^4, T = 170, r0 = 1.477, L=n/(4/3*pi*r^3)
  2. a = 1/(20*r0^3), mu = ((6*L*pi^2)/g)^1/3,
  3. ErĀ® = 3*pr(r)+4B-16/27*mu*T(mu^2-2*mu*T+T^2)
  4.  
  5. s =
  6.   NDSolve[
  7.     {m'[r] == a*r^2 Er[r],
  8.      pr'[r] == -(r0/r)*((pr[r] +Er[r]) (m[r] + a*r^3* pr[r]))/(r - 2 r0*m[r]),
  9.      m[1] == 1, pr[1] == 1},
  10.     {m, pr}, {r, 1, rmax}]
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top