Graph Plotting Problem in Mathematica

Click For Summary
SUMMARY

The discussion focuses on graph plotting in Mathematica, specifically reproducing Figure 11 from Lattimer's paper on binary systems. The user is attempting to plot the neutron star radius versus neutron star mass and the derivative of the radius with respect to mass. They have encountered difficulties achieving the desired results and are advised to simplify their approach by first plotting basic functions before incrementally adding complexity. This method allows for easier debugging and understanding of the plotting process.

PREREQUISITES
  • Familiarity with Mathematica 12.0 syntax and functions
  • Understanding of derivatives in calculus
  • Knowledge of binary star systems and related astrophysical concepts
  • Ability to interpret and implement mathematical equations in code
NEXT STEPS
  • Learn how to use Mathematica's Plot function effectively
  • Study the implementation of derivatives in Mathematica using the D function
  • Explore examples of plotting complex functions in Mathematica
  • Review Lattimer's paper for deeper insights into the equations used
USEFUL FOR

Astronomy researchers, astrophysicists, and students working with Mathematica for plotting complex mathematical functions related to binary systems.

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
2K
  • · Replies 48 ·
2
Replies
48
Views
6K
  • · 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
4K