Table NDSolve Plot
- To: mathgroup at smc.vnet.net
- Subject: [mg121885] Table NDSolve Plot
- From: Howie <hcohl001 at gmail.com>
- Date: Thu, 6 Oct 2011 04:19:20 -0400 (EDT)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
I know this is a basic question, but I am having a hard time figuring
this out on my own.
I would like to numerically solve an ordinary differential equation
with a parameter which varies over a range of numbers and plot all the
solutions with different values of the parameter on the same plot.
Also, I would like to be able to then over-plot another solution on
the same plot.
For instance:
h=0
n=0
k=0.25
s = NDSolve[{(h - k^2 n (1 + n) z^2) w[z] -
z (1 + k^2 - 2 k^2 z^2) Derivative[1][w][
z] + (1 - z^2) (1 - k^2 z^2) (w^\[Prime]\[Prime])[z] == 0,
w[2] == 1, w'[2] == -1}, w, {z, 1.1, 3.0}]
Plot[Evaluate[w[z] /. s], {z, 1.1, 3.0}, PlotStyle -> Automatic]
Plot[1/2 Log[(z + 1)/(z - 1)], {z, 1.1, 3.0}, PlotStyle -> Automatic]
I would like to do this for instance 10 increments of k from 0 to 1,
and also the Log plot, all on the same plot.
Thanks!