Another silly question: is sqrt(3) an integer



> From:  Robert Marik
> 
> Hello is this known problem?
> 
> load('contrib_ode);
> contrib_ode('diff(y,x,2)*x^2+'diff(y,x)*x+(x^2-3)*y=0,y,x);
> 
> asks if sqrt(3) is integer

Thank you for the bug report.  I have had a look and this is actually a 
problem in the ode2() function. 

In share/diffequations/ode/ode2.mac we have:  

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,
   if nu = 1/2 then return(y = (%k1*sin(x-pt) +
%k2*cos(x-pt))/sqrt(x-pt)),
   if featurep(nu,'integer) then intp: 'y else if numberp(nu) then intp:
'n,
 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*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)))$



I thought that the askinteger() function might be part of the solution.
However:

(%i1) askinteger(sqrt(3));
Is sqrt(3) an integer?
no;
(%o1)                                 no

This email is confidential and may also be privileged.  If you are not the intended recipient, please notify us immediately and delete this message from your system without first printing or copying it. Any personal data in this email (including any attachments) must be handled in accordance with the Rio Tinto Group Data Protection Policy and all applicable data protection laws.