Changeset e10e40
 Timestamp:
 03/23/12 04:13:07 (3 years ago)
 Branches:
 master
 Children:
 fe8cff
 Parents:
 1d9ec0
 Files:

 2 edited
Legend:
 Unmodified
 Added
 Removed

qdgamma.lisp
r1d9ec0 re10e40 502 502 (gamma v)) 503 503 ( (psi v) (log z))) 504 (loop for k from 0 below n 505 for term = 1 then (* term (/ z k)) 506 for sum = (/ ( 1 v)) then (+ sum (/ term (+ k 1 v))) 507 when (< (abs term) (* (abs sum) eps)) 508 return sum) 509 (loop for k from n 510 for term = 1 then (* term (/ z k)) 511 for sum = 0 then (+ sum (/ term (+ k 1 v))) 512 when (< (abs term) (* (abs sum) eps)) 513 return sum))) 504 (loop for k from 0 505 for term = 1 then (* term (/ z k)) 506 for sum = (/ ( 1 v)) then (+ sum (let ((denom (+ k 1 ( n)))) 507 (if (zerop denom) 508 0 509 (/ term (+ k 1 v))))) 510 when (< (abs term) (* (abs sum) eps)) 511 return sum))) 514 512 (loop for k from 0 515 513 for term = 1 then (* term (/ z k)) 
rttests.lisp
r1d9ec0 re10e40 1470 1470 (checkaccuracy 210.8 e true)) 1471 1471 nil) 1472 1473 (rt:deftest expintegrale.5d 1474 (let* ((x .5d0) 1475 (e (expintegrale 2d0 x)) 1476 (true #q0.3266438623245530177304015653336378358284946903290101)) 1477 (checkaccuracy 51.2 e true)) 1478 nil) 1479 1480 (rt:deftest expintegrale.5q 1481 (let* ((x #q.5) 1482 (e (expintegrale #q2 x)) 1483 (true #q0.326643862324553017730401565333637835828494690329010198058745549181386569998611289568)) 1484 (checkaccuracy 208.4 e true)) 1485 nil)
Note: See TracChangeset
for help on using the changeset viewer.