Mathematica Graph Plotting Problem in Mathematica

Click For Summary
The discussion revolves around troubleshooting graph plotting in Mathematica, specifically aiming to reproduce Figure 11 from Lattimer's paper on binary systems. The user is attempting to plot the neutron star radius against its mass, using the defined equations for the neutron star mass-to-black hole mass ratio (q) and the radius in relation to the semi-major axis. They have encountered difficulties achieving the expected results and seek assistance. A suggestion is made to simplify the problem by first plotting derivatives of simpler functions to ensure the plotting process works correctly before gradually increasing complexity. The user acknowledges this approach and expresses intent to try it again.
Stella.Physics
Messages
62
Reaction score
13
Hello to all. I've been trying to plot some graphs on Mathematica but I have faced some troubles so far.

I am working on this paper by Lattimer on binary systems:

https://arxiv.org/abs/astro-ph/0002203

I am trying to reproduce Figure 11 but something seems to go wrong.

Here's what I have done so far:

we set ##q=\frac{m_{ns}}{M_{bh}}## the ratio of the neutron star mass and the black hole mass.

then we set the radius

##\frac{R_{l}}{a} = 0.46\bigg(\frac{q}{1+q}\bigg)^{1/3}## with a being the semi major axis that follows
## a<2.8 \bigg[ \frac{m_{ns} m_{bh} (m_{ns}+m_{bh})}{〖M^3_ {\odot}〗} \bigg]^{\frac{1}{4}} R_{\odot}##

The first thing to plot is the neutron star radius versus the neutron star mass.

Then we must plot ##\frac{d \ln{R}}{d \ln{m_{ns}}} ## versus ## m_{ns}## which derivative simply leads to ##\frac{m_{ns}}{R} \frac{dR}{dm_{ns}}##

I've been trying to plot these functions for days but what I get is far from the desired results...

Any help would be very appreciated :)
 
Physics news on Phys.org
It's a good idea to simplify a complicated problem, reducing it to it's basic concept and just get the basics working then add on step by step, pieces that would make the basic look more and more like the problem you're trying to solve. So just a quick overview of what you're trying to do looks like plotting a derivative. So if this was mine, I would first get working, plotting the derivative of a simple function like:

Code:
myf[x_] := x^2 + 3 x + 4;
myd[x_] = D[myf[x], x]
Plot[myd[x], {x, 0, 5}]

Now, I might add a more complicated function like myf[x_]:=Sin[x]+Log[x] and get that working and so on until I guide my code towards the problem I really want to solve so that if I just increment it a little bit by bit and something goes wrong I can go back a step and undo the last change I made then analyze why the newest changed caused my code to fail.
 
Thank you for your reply, aheight! Although I have tried this method, I will give it another shot... Thanks again :)
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 48 ·
2
Replies
48
Views
5K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
5K
  • · Replies 44 ·
2
Replies
44
Views
4K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 17 ·
Replies
17
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 39 ·
2
Replies
39
Views
3K