Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- p1 = ParametricPlot[{E^-k BesselI[0, k], BesselI[1, k]/BesselI[0, k]}
- , {k, 0, 1000}
- , PlotPoints -> 200
- , PlotRange -> {{0, 1}, All}]
- pts = First@Cases[Normal@p1, Line[a_] :> a, Infinity];
- f[d_, t_] := d Interpolation[pts][t]
- Plot3D[f[d, t], {d, 0, 1}, {t, 0.013, 1}, ViewPoint -> {2.5, 1.5, 1.5}]
- E^-k BesselI[0, k] /. k -> 1000.
- (* 0.01261724045589 *)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement