Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- (*First ODE gives y as output*)
- s = NDSolve[{y'[x] == 0, y[0] == 1}, y, {x, 0, 1}]
- yy[x_] := Evaluate[y[x] /. s]
- Plot[yy[x], {x, 0, 2}, PlotRange -> All]
- (*Second ODE must use y as input*)
- s2 = NDSolve[{z'[x] == yy[x], z[0] == 0}, z, {x, 0, 1}]
- zz[x_] := Evaluate[z[x] /. s2]
- Plot[zz[x], {x, 0, 2}, PlotRange -> All]
- s = NDSolve[{y'[x] == 0, y[0] == 1}, y, {x, 0, 1}]
- yy[x_] := Evaluate[y[x] /. s] // First
- Plot[yy[x], {x, 0, 2}, PlotRange -> All]
- s2 = NDSolve[{z'[x] == yy[x], z[0] == 0}, z, {x, 0, 1}]
- zz[x_] := Evaluate[z[x] /. s2]
- Plot[zz[x], {x, 0, 2}, PlotRange -> All]
Add Comment
Please, Sign In to add comment