From: David B. <bil...@us...> - 2006-09-04 12:27:19
|
Update of /cvsroot/maxima/maxima/share/diffequations In directory sc8-pr-cvs7.sourceforge.net:/tmp/cvs-serv20849/share/diffequations Modified Files: ode2.mac Log Message: Use maxima language file ode2.mac in place of translated file ode2.lisp. Remove ode2.lisp and all references to it in the build process. Change Bessel functions from obsolete form %j[n](x) and %y[n](x) to bessel_j(n,x) and bessel_y(n,x) in ode2.mac and in ode2 solutions in testsuite files. Autoload ode2.mac on references to ode2(), ic1(), ic2() and bc2(). Index: ode2.mac =================================================================== RCS file: /cvsroot/maxima/maxima/share/diffequations/ode2.mac,v retrieving revision 1.2 retrieving revision 1.3 diff -u -d -r1.2 -r1.3 --- ode2.mac 9 Jan 2005 07:44:54 -0000 1.2 +++ ode2.mac 4 Sep 2006 12:27:16 -0000 1.3 @@ -282,7 +282,7 @@ dc: -dc, ip: sqrt(dc)/2, return(y = (x-pt)^rp * (%k1*sin(ip*log(x-pt)) + %k2*cos(ip*log(x-pt)))))$ -bessel2(a,b):=block([nu,b1,intp,radexpand:'all,%k1,%y,%k2,%j], +bessel2(a,b):=block([nu,b1,intp,radexpand:'all,%k1,%k2], if not(freeof(x,y,b1: ratsimp((1-b)*(x-pt)^2))) then return(false), if ratsimp(a*(x-pt)) # 1 then return(false), nu: sqrt(b1), method: 'bessel, @@ -291,8 +291,8 @@ loop, if not(intp='y or intp='n) then (intp: readonly("is", nu, "an integer? type y or n."), go(loop)), - if intp = 'y then return(y = %k1*%j[nu](x-pt) + %k2*%y[nu](x-pt)), - return(y = %k1*%j[nu](x-pt) + %k2*%j[-nu](x-pt)))$ + if intp = 'y then return(y = %k1*bessel_j(nu,x-pt) + %k2*bessel_y(nu,x-pt)), + return(y = %k1*bessel_j(nu,x-pt) + %k2*bessel_j(-nu,x-pt)))$ ic1(soln,xc,yc):= block([%c], |