Advertisement
Guest User

Untitled

a guest
Oct 6th, 2015
82
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.37 KB | None | 0 0
  1. p1 = ParametricPlot[{E^-k BesselI[0, k], BesselI[1, k]/BesselI[0, k]}
  2. , {k, 0, 1000}
  3. , PlotPoints -> 200
  4. , PlotRange -> {{0, 1}, All}]
  5.  
  6. pts = First@Cases[Normal@p1, Line[a_] :> a, Infinity];
  7.  
  8. f[d_, t_] := d Interpolation[pts][t]
  9.  
  10. Plot3D[f[d, t], {d, 0, 1}, {t, 0.013, 1}, ViewPoint -> {2.5, 1.5, 1.5}]
  11.  
  12. E^-k BesselI[0, k] /. k -> 1000.
  13. (* 0.01261724045589 *)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement