From: Raymond T. <rt...@us...> - 2010-01-15 12:44:56
|
Update of /cvsroot/maxima/maxima/src In directory sfp-cvsdas-1.v30.ch3.sourceforge.com:/tmp/cvs-serv16811 Modified Files: cpoly.lisp Log Message: Fix typo in QUADIT-SL. The loop for computing the estimated error only had one term in it instead of adding up terms with all of the polynomial coefficients. Index: cpoly.lisp =================================================================== RCS file: /cvsroot/maxima/maxima/src/cpoly.lisp,v retrieving revision 1.20 retrieving revision 1.21 diff -u -d -r1.20 -r1.21 --- cpoly.lisp 12 Jan 2010 21:45:45 -0000 1.20 +++ cpoly.lisp 15 Jan 2010 12:44:35 -0000 1.21 @@ -650,7 +650,7 @@ zm (sqrt (abs *v*)) ee (* 2.0 (abs (aref *qpr-sl* 0))) t1 (- (* *szr* *b*))) - (do ((i 1 (1+ *n*))) + (do ((i 1 (1+ i))) ((> i *n*)) (setq ee (+ (* ee zm) (abs (aref *qpr-sl* i))))) (setq ee (+ (* ee zm) (abs (+ *a* t1))) ee (- (* (+ (* 5.0 *mre*) (* 4.0 *are*)) ee) |