Services & ResourcesWolfram ForumsMathGroup Archive

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

Re: Schaums outline Mathematica 2nd ed, prob. 11.3 help

• To: mathgroup at smc.vnet.net
• Subject: [mg125737] Re: Schaums outline Mathematica 2nd ed, prob. 11.3 help
• From: Bob Hanlon <hanlonr357 at gmail.com>
• Date: Fri, 30 Mar 2012 04:38:06 -0500 (EST)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• References: <201203290759.CAA21786@smc.vnet.net>

First, in the definition of s1 replace X with x.

Do does not provide an output. Either wrap its Plot with Print
(behavior of Plot changed since book was written) to cause an output
or replace the Do with Table.

Do[
 Print[g[k] = Plot[s1[[k, 1, 1, 2]],
    {x, 0, 1},
    PlotStyle -> Thickness[0.007],
    PlotRange -> All]],
 {k, 1, 5}]

or

Table[
 g[k] = Plot[s1[[k, 1, 1, 2]],
   {x, 0, 1},
   PlotStyle -> Thickness[0.007],
   PlotRange -> All],
 {k, 1, 5}]


Bob Hanlon

On Thu, Mar 29, 2012 at 3:59 AM, pennsylvaniajake at gmail.com
<pennsylvaniajake at gmail.com> wrote:
> I cannot get this code to show the solutions to the equation, only the vector plot.  Whats wrong?  It came from Schum's Outline, Mathematica 2n ed, solved problem 11.3 on page 272-3.  The only changes I made was using v1 for vf and s1 for solutions.
>
> v1 = VectorPlot[{1, x^2 + y}, {x, 0, 1}, {y, 0, 12},
>   Axes -> Automatic, VectorScale -> Tiny];
>
> s1 = Table[
>   DSolve[{y'[x] == X^2 + y[x], y[0] == k}, y[x], x], {k, 0, 4}];
>
> Do[g[k] =
>  Plot[s1[[k, 1, 1, 2]], {x, 0, 1}, PlotStyle -> Thickness[0.007],
>   PlotRange -> All], {k, 1, 5}]
>
> Show[g[1], g[2], g[3], g[4], g[5], v1, Frame -> True,
>  AspectRatio -> 1]
>

• References:
  • Schaums outline Mathematica 2nd ed, prob. 11.3 help
    • From: "pennsylvaniajake@gmail.com" <pennsylvaniajake@gmail.com>