MathGroup Archive 2004

[Date Index] [Thread Index] [Author Index]

Search the Archive

LogLog plot of NDSolve solution

  • To: mathgroup at smc.vnet.net
  • Subject: [mg52557] LogLog plot of NDSolve solution
  • From: Adam Getchell <agetchell at physics.ucdavis.edu>
  • Date: Thu, 2 Dec 2004 02:21:06 -0500 (EST)
  • Sender: owner-wri-mathgroup at wolfram.com

Hello all,

I'm solving an ODE via NDSolve (neutron abundance for Big Bang 
Nucleosynthesis -- from "Modern Cosmology", S. Dodelson). The result 
comes out as an InterpolatingFunction, which plots fine using Plot. 
However, if I wish to do a Log-Log plot, neither LogLogPlot nor 
LogLogListPlot work. What am I missing?

\!\(Dodelson := \(\[CapitalChi]'\)[
        x] - \(x\/H\) \(255\/\(\(\[Tau]\_n\) x\^5\)\) \((12 + 6  x +
            x\^2)\) \((\[ExponentialE]\^\(-x\) - \[CapitalChi][
                x] \((1 + \[ExponentialE]\^\(-x\))\))\)\)

In[329]:=
\!\(\(\(\[IndentingNewLine]\)\(ans =
    NDSolve[{Dodelson \[Equal]
            0 /. {H \[Rule] 1.13, \[Tau]\_n \[Rule] 886.7}, \[CapitalChi][
            1.293] \[Equal]  .25}, \[CapitalChi], {x,  .1, 1.3},
      MaxSteps \[Rule] 10\^6]\)\)\)

Using the Evaluate trick suggested by 
http://www.nhn.ou.edu/~morrison/Mathematica/TipSheets/2Dgraphics.pdf 
doesn't work.

This example doesn't cover LogLogPlots: 
http://documents.wolfram.com/v4/MainBook/3.8.2.html

Other search results on Google haven't turned up anything so far.

Thanks,

--Adam Getchell



  • Prev by Date: Invoking a function using full package name fail, but works without full package name
  • Next by Date: JLink to remote kernel not possible?
  • Previous by thread: Re: Invoking a function using full package name fail, but works without full package name
  • Next by thread: Re: LogLog plot of NDSolve solution