home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #3 / NN_1993_3.iso / spool / sci / math / numanal / 3974 < prev    next >
Encoding:
Internet Message Format  |  1993-01-28  |  1.9 KB
  1. Xref: sparky sci.math.num-analysis:3974 sci.engr:2538 sci.engr.mech:845
  2. Path: sparky!uunet!portal!lll-winken!uwm.edu!zaphod.mps.ohio-state.edu!howland.reston.ans.net!paladin.american.edu!gatech!rpi!ghost.dsi.unimi.it!univ-lyon1.fr!scsing.switch.ch!bernina!nmr.lpc.ethz.ch!tilo
  3. From: tilo@nmr.lpc.ethz.ch (Tilo Levante)
  4. Newsgroups: sci.math.num-analysis,sci.engr,sci.engr.mech
  5. Subject: Re: Algorithm(s) need for Bessel functions
  6. Message-ID: <1993Jan26.172944.26069@bernina.ethz.ch>
  7. Date: 26 Jan 93 17:29:44 GMT
  8. References: <1993Jan25.201447.2082@athena.mit.edu> <1993Jan26.100935.113272@eratu.rz.uni-konstanz.de>
  9. Sender: news@bernina.ethz.ch (USENET News System)
  10. Organization: NMR-Group, Physical Chemistry Lab., ETH-Zuerich, Switzerland
  11. Lines: 34
  12.  
  13. In article <1993Jan26.100935.113272@eratu.rz.uni-konstanz.de>, vaxinf@V36.CHEMIE.UNI-KONSTANZ.DE () writes:
  14. |> 
  15. |> In article <1993Jan25.201447.2082@athena.mit.edu>, cai@athena.mit.edu (Liang-Wu Cai) writes:
  16. |> |>Dear netter,
  17. |> |>
  18. |> |> I am in need (badly) of an accurate algorithm Bessel functions of various 
  19. |> |>kinds with real argument and integer order. What I need actually is a 
  20. |> |>algorithm which is able to provide the result with at least 7 significant 
  21. |> |>digits for J_n (z) and Y_n (z) with n <=100 and z <= 100.
  22. |> 
  23. |> Look into Abromowitz/Stegun Handbook of Mathamatical Functions:
  24. |> 
  25. |> 1.) set j(100)=0
  26. |> 2.) set j(99)=1
  27. |> 3.) do recursion to j(0)
  28. |> 4.) calculate the true j(0) in order to define a normalisation faktor
  29. |> 5.) Multiply all value with that factor
  30. |> 
  31. |> You get all values J(0)...j(100) in that way.
  32. |> 
  33. |> Eberhard Heuser-Hofmann
  34.  
  35. This algorithm (Miller algorithm)
  36. is based on a recurrence relation for besselfunctions.
  37. You can find a detailed description in 
  38.  
  39.  Computation with recurrence relations
  40.  / Jet Wimp
  41.  Boston a.o. : Pitman, 1984
  42.  XII, 310 p. : tab. ; 24 cm
  43.  (Applicable mathematics series)
  44.  ISBN 0-273-08508-5
  45.  
  46. Tilo
  47.