Changeset d795ba
 Timestamp:
 04/12/12 03:15:47 (3 years ago)
 Branches:
 master
 Children:
 6ab522
 Parents:
 c0e12d
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

qdbessel.lisp
rc0e12d rd795ba 376 376 (format t " term = ~S~%" term) 377 377 (format t " sum = ~S~%" sum)))))) 378 378 379 ;; TODO: 380 ;; o For z <= 1 use the series. 381 ;; o Currently accuracy is not good for large z and halfinteger 382 ;; order. 383 ;; o For real v and z, return a real number instead of complex. 384 ;; o Handle the case of Re(z) < 0. (The formulas are for Re(z) > 0: 385 ;; bessel_j(v,z*exp(m*%pi*%i)) = exp(m*v*%pi*%i)*bessel_j(v, z) 386 ;; o The paper suggests using 387 ;; bessel_i(v,z) = exp(v*%pi*%i/2)*bessel_j(v, %i*z) 388 ;; when Im(z) >> Re(z) 389 ;; 379 390 (defun besselj (v z) 380 391 (let ((vv (ftruncate v))) … … 383 394 (integerbesseljexparc v z)) 384 395 (t 396 ;; Need to finetune the value of bign. 385 397 (let ((bign 100) 386 398 (vpi (* v (floatpi (realpart z)))))
Note: See TracChangeset
for help on using the changeset viewer.