Hallo Gunter,
Erst nach seiner Korrektur loeste sich alles so schoen auf!
Richtig.
Oder anders gesagt, das Skript, das ich nachzuvollziehen versuchte,
enthielt doch Fehler, vielleicht nur Darstellungsfehler.
Die Gleichungen Ihres Skriptes bestehen die einfachste Sichtprobe nicht
In[7]:= With[{q = 4},
Plot[{(BesselJ[-1 + q, u] - BesselJ[1 + q, u])/(
2 u BesselJ[q, u]) - (q/u^2 + BesselJ[1 + q, u]/(u BesselJ[q, u])),
(BesselJ[-1 + q, u] - BesselJ[1 + q, u])/(2 u BesselJ[q, u]) + q/
u^2 + BesselJ[-1 + q, u]/(u BesselJ[q, u])} , {u, 1, 8}]
]
die von Markus natürlich schon
In[11]:= With[{q = 4},
Plot[{(BesselJ[-1 + q, u] - BesselJ[1 + q, u])/(2 u BesselJ[q, u]) -
q/u^2 + BesselJ[1 + q, u]/(u BesselJ[q, u]),
(BesselJ[-1 + q, u] - BesselJ[1 + q, u])/(2 u BesselJ[q, u]) + q/
u^2 - BesselJ[-1 + q, u]/(u BesselJ[q, u])} , {u, 1, 8}]
]
alle ausgedruckten Werte sind kleiner als $MachineEpsilon und den
Gleichungen von Markus sieht man es -wie gesagt- sofort an, dass sie
richtig sind.
Gruss
Udo.
DMUG-Archiv, http://www.mathematica.ch/archiv.html