More probs,
jacobi_p(1,1/2,-2,1.125); --> interval(1.53125,2.6797208657787516*10^-15)
This would be good if everyone else agreed about intervals and knew how
to process them, but they don't, yet.
Right now there appears to be no mention of interval in the documentation.
jacobi_p(1,1/2,-2,1,x);
%,x=1.125 --> 1.53125
which happens to be what I wanted.
where did this come from? consider:
i3(m,k,a,x):=(x^(1+m)*(a^2+x^2)^k*hgfred([(1+m)/2,-k],[1+(1+m)/2],-x^2/a^2))/((1+m)*(1+x^2/a^2)^k);
then z: i3(0,1,2,0,0.5) produces
1.333333333333333*jacobi_p(1.0,0.5,-2.0,1.125).
But it should have computed the value of jacobi_p, no??
After trying many ways to evaluate this, the first that worked was
ratsimp(z),
%,numer.
but that returned an interval expression.
huh?
Also, at the risk of putting too many disparate complaints in the same
message...
? hypergeometric might have some kind of entry in the on-line info.
?? hypergeometric might have a link to hgfred.
Thanks
RJF