## [Maxima-commits] CVS: maxima/share/numeric dblint_1.dem,1.1.1.1,1.2 kach.dem,1.1.1.1,1.2

 [Maxima-commits] CVS: maxima/share/numeric dblint_1.dem,1.1.1.1,1.2 kach.dem,1.1.1.1,1.2 From: Robert Dodier - 2005-01-26 02:12:21 ```Update of /cvsroot/maxima/maxima/share/numeric In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv28090/numeric Modified Files: dblint_1.dem kach.dem Log Message: Changing uppercase to lowercase in some demo files. Uppercase in comments is unchanged. All demo files were reviewed at this time but most of them are already in lowercase. Index: dblint_1.dem =================================================================== RCS file: /cvsroot/maxima/maxima/share/numeric/dblint_1.dem,v retrieving revision 1.1.1.1 retrieving revision 1.2 diff -u -d -r1.1.1.1 -r1.2 --- dblint_1.dem 8 May 2000 06:09:43 -0000 1.1.1.1 +++ dblint_1.dem 26 Jan 2005 02:09:44 -0000 1.2 @@ -5,65 +5,65 @@ FLOATDEFUNK and QUANC8 */ /* Get the fasl file for DBLINT */ -LOAD('DBLINT); +load('dblint); /* Get the fasl file for QUANC8 */ -LOAD('QQ); +load('qq); /* Use DBLINT to get the double integral of exp(x-y^2) over the region bounded by y=1 and y=2+x^(3/2) and x=0 to x=1 */ /* Define the integrand as a function of two variables */ -F(X,Y):=(MODE_DECLARE([X,Y],FLOAT),EXP(X-Y^2)); +f(x,y):=(mode_declare([x,y],float),exp(x-y^2)); /* Define the lower and upper limits on the inner (y in this case) integral as a function of the outer variable (x in this case) */ -R(X):=(MODE_DECLARE(X,FLOAT),1.0); -S1(X):=(MODE_DECLARE(X,FLOAT),2.0+X^(3/2)); +r(x):=(mode_declare(x,float),1.0); +s1(x):=(mode_declare(x,float),2.0+x^(3/2)); /* Now translate these functions for the sake of efficiency */ -TRANSLATE(F,R,S1); +translate(f,r,s1); /* Call the DBLINT function with quoted arguments for function names, and floating point values for the endpoints of the outer (x) integration */ -DBLINT_ANSWER:DBLINT('F,'R,'S1,0.0,1.0); +dblint_answer:dblint('f,'r,'s1,0.0,1.0); /* Now generate the exact integral over y using RISCH */ -INTY:RISCH(EXP(X-Y^2),Y); +inty:risch(exp(x-y^2),y); /* Now get the integrand of the x integral */ -XINT:EV(INTY,Y:2+X^(3/2))-EV(INTY,Y:1); +xint:ev(inty,y:2+x^(3/2))-ev(inty,y:1); /* Try to do the x integral exactly */ -RISCH(XINT,X); -RADCAN(%); -%,NOUNS; +risch(xint,x); +radcan(%); +%,nouns; /* Still no luck with closed-form */ /* The integral over x can't be done in closed-form, so let's use a one-dimensional integrator to do it numerically. First, we must define the integrand as a floating point function. FLOATDEFUNK does this very conveniently on the expression */ -FLOATDEFUNK(fUN,[X],XINT); +floatdefunk(fun,[x],xint); /* Now translate the function fun */ -TRANSLATE(FUN); +translate(fun); /* Now call QUANC8 on the translated function (see SHARE1;QQ USAGE for details) */ -QUANC8_ANSWER:QUANC8(FUN,0.0,1.0); +quanc8_answer:quanc8(fun,0.0,1.0); /* Compare the answers from DBLINT and QUANC8 */ -(DBLINT_ANSWER-QUANC8_ANSWER)/DBLINT_ANSWER; +(dblint_answer-quanc8_answer)/dblint_answer; /* It appears to be small enough relative error for most uses */ "end of dblint.dm1"\$ Index: kach.dem =================================================================== RCS file: /cvsroot/maxima/maxima/share/numeric/kach.dem,v retrieving revision 1.1.1.1 retrieving revision 1.2 diff -u -d -r1.1.1.1 -r1.2 --- kach.dem 8 May 2000 06:09:43 -0000 1.1.1.1 +++ kach.dem 26 Jan 2005 02:09:44 -0000 1.2 @@ -4,8 +4,8 @@ load("kach")\$ load("haipart")\$ fpprec:16\$ -If Status(Feature,ITS)=True -Then Translate(hach,check,iterate,fx,fq,fh,getl,criterion,accelerate)\$ +if status(feature,ITS)=true +then translate(hach,check,iterate,fx,fq,fh,getl,criterion,accelerate)\$ /* dynamalloc:true\$ */ ```

 [Maxima-commits] CVS: maxima/share/numeric dblint_1.dem,1.1.1.1,1.2 kach.dem,1.1.1.1,1.2 From: Robert Dodier - 2005-01-26 02:12:21 ```Update of /cvsroot/maxima/maxima/share/numeric In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv28090/numeric Modified Files: dblint_1.dem kach.dem Log Message: Changing uppercase to lowercase in some demo files. Uppercase in comments is unchanged. All demo files were reviewed at this time but most of them are already in lowercase. Index: dblint_1.dem =================================================================== RCS file: /cvsroot/maxima/maxima/share/numeric/dblint_1.dem,v retrieving revision 1.1.1.1 retrieving revision 1.2 diff -u -d -r1.1.1.1 -r1.2 --- dblint_1.dem 8 May 2000 06:09:43 -0000 1.1.1.1 +++ dblint_1.dem 26 Jan 2005 02:09:44 -0000 1.2 @@ -5,65 +5,65 @@ FLOATDEFUNK and QUANC8 */ /* Get the fasl file for DBLINT */ -LOAD('DBLINT); +load('dblint); /* Get the fasl file for QUANC8 */ -LOAD('QQ); +load('qq); /* Use DBLINT to get the double integral of exp(x-y^2) over the region bounded by y=1 and y=2+x^(3/2) and x=0 to x=1 */ /* Define the integrand as a function of two variables */ -F(X,Y):=(MODE_DECLARE([X,Y],FLOAT),EXP(X-Y^2)); +f(x,y):=(mode_declare([x,y],float),exp(x-y^2)); /* Define the lower and upper limits on the inner (y in this case) integral as a function of the outer variable (x in this case) */ -R(X):=(MODE_DECLARE(X,FLOAT),1.0); -S1(X):=(MODE_DECLARE(X,FLOAT),2.0+X^(3/2)); +r(x):=(mode_declare(x,float),1.0); +s1(x):=(mode_declare(x,float),2.0+x^(3/2)); /* Now translate these functions for the sake of efficiency */ -TRANSLATE(F,R,S1); +translate(f,r,s1); /* Call the DBLINT function with quoted arguments for function names, and floating point values for the endpoints of the outer (x) integration */ -DBLINT_ANSWER:DBLINT('F,'R,'S1,0.0,1.0); +dblint_answer:dblint('f,'r,'s1,0.0,1.0); /* Now generate the exact integral over y using RISCH */ -INTY:RISCH(EXP(X-Y^2),Y); +inty:risch(exp(x-y^2),y); /* Now get the integrand of the x integral */ -XINT:EV(INTY,Y:2+X^(3/2))-EV(INTY,Y:1); +xint:ev(inty,y:2+x^(3/2))-ev(inty,y:1); /* Try to do the x integral exactly */ -RISCH(XINT,X); -RADCAN(%); -%,NOUNS; +risch(xint,x); +radcan(%); +%,nouns; /* Still no luck with closed-form */ /* The integral over x can't be done in closed-form, so let's use a one-dimensional integrator to do it numerically. First, we must define the integrand as a floating point function. FLOATDEFUNK does this very conveniently on the expression */ -FLOATDEFUNK(fUN,[X],XINT); +floatdefunk(fun,[x],xint); /* Now translate the function fun */ -TRANSLATE(FUN); +translate(fun); /* Now call QUANC8 on the translated function (see SHARE1;QQ USAGE for details) */ -QUANC8_ANSWER:QUANC8(FUN,0.0,1.0); +quanc8_answer:quanc8(fun,0.0,1.0); /* Compare the answers from DBLINT and QUANC8 */ -(DBLINT_ANSWER-QUANC8_ANSWER)/DBLINT_ANSWER; +(dblint_answer-quanc8_answer)/dblint_answer; /* It appears to be small enough relative error for most uses */ "end of dblint.dm1"\$ Index: kach.dem =================================================================== RCS file: /cvsroot/maxima/maxima/share/numeric/kach.dem,v retrieving revision 1.1.1.1 retrieving revision 1.2 diff -u -d -r1.1.1.1 -r1.2 --- kach.dem 8 May 2000 06:09:43 -0000 1.1.1.1 +++ kach.dem 26 Jan 2005 02:09:44 -0000 1.2 @@ -4,8 +4,8 @@ load("kach")\$ load("haipart")\$ fpprec:16\$ -If Status(Feature,ITS)=True -Then Translate(hach,check,iterate,fx,fq,fh,getl,criterion,accelerate)\$ +if status(feature,ITS)=true +then translate(hach,check,iterate,fx,fq,fh,getl,criterion,accelerate)\$ /* dynamalloc:true\$ */ ```