[Maxima-commits] CVS: maxima/tests rtest16.mac, 1.68, 1.69 rtesthyp.mac, 1.48, 1.49

[Raymond T. <rt...@us...>]

Update of /cvsroot/maxima/maxima/tests
In directory 23jxhf1.ch3.sourceforge.com:/tmp/cvs-serv3529/tests

Modified Files:
	rtest16.mac rtesthyp.mac 
Log Message:
Bug [ 2687962 ] hgfred([-3/2,1],[-1/2],-t) division by zero

src/hyp.lisp:
o Add implementation of A&S 15.2.7 and 15.2.8.
o For completeness add A&S 15.2.2 and 15.2.9.  Not currently used.
o Fix bug in HYP-ATANH by using 15.2.7 and 15.2.8 so we adjust both a
  and c by the same amount.

tests/rtest16.mac:
o Add test case for bug.

tests/rtesthyp.mac:
o If rtesthyp is run by itself, there are several failures.  Fix these
  failures by setting triginverses to all for the tests that were
  failing.  (Some earlier test must be setting triginverses:all before
  rtesthyp is run.)


Index: rtest16.mac
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/rtest16.mac,v
retrieving revision 1.68
retrieving revision 1.69
diff -u -d -r1.68 -r1.69
--- rtest16.mac	19 Mar 2009 16:26:45 -0000	1.68
+++ rtest16.mac	20 Mar 2009 03:40:21 -0000	1.69
@@ -645,5 +645,13 @@
 float((2^1000-1)/2^1000)-1;
 0.0;
 
+/*
+ * Bug [ 2687962 ] hgfred([-3/2,1],[-1/2],-t) division by zero
+ *
+ * Solution from functions.wolfram.com
+ */
+ratsimp(hgfred([-3/2,1],[-1/2], t));
+1+3*t-3*t^(3/2)*atanh(sqrt(t));
+
 end;
 end$

Index: rtesthyp.mac
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/rtesthyp.mac,v
retrieving revision 1.48
retrieving revision 1.49
diff -u -d -r1.48 -r1.49
--- rtesthyp.mac	23 Oct 2008 11:13:52 -0000	1.48
+++ rtesthyp.mac	20 Mar 2009 03:40:21 -0000	1.49
@@ -131,29 +131,29 @@
 2^(2*a-1)/sqrt(1-z)*(1+sqrt(1-z))^(1-2*a)$
 
 /* A&S 15.1.15 */
-hgfred([a,1-a],[3/2],sin(z)^2);
+hgfred([a,1-a],[3/2],sin(z)^2),triginverses:all;
 sin((1-2*a)*z)/sin(z)/(1-2*a)$
 
 /* A&S 15.1.16 */
-hgfred([a,2-a],[3/2],sin(z)^2);
+hgfred([a,2-a],[3/2],sin(z)^2),triginverses:all;
 sin((2-2*a)*z)/(2-2*a)/(sin(z)*cos(z))$
 
 /* A&S 15.1.17 */
 assume(not(equal(sin(z),0)));
 [not equal(sin(z),0)]$
 
-hgfred([-a,a],[1/2],sin(z)^2);
+hgfred([-a,a],[1/2],sin(z)^2),triginverses:all;
 cos(2*a*z)$
 
 /* A&S 15.1.18 */
-hgfred([a,1-a],[1/2],sin(z)^2);
+hgfred([a,1-a],[1/2],sin(z)^2),triginverses:all;
 cos((2*a-1)*z)/cos(z)$
 
 /* A&S 15.1.19 */
 assume(not(equal(tan(z),0)));
 [not equal(tan(z), 0)]$
 
-hgfred([a,a+1/2],[1/2],-tan(z)^2);
+hgfred([a,a+1/2],[1/2],-tan(z)^2),triginverses:all;
 cos(z)^(2*a)*cos(2*a*z)$
 
 /* A&S 15.4.3 */