Hi Liam, On Saturday 22 March 2008 20:15:06 Liam Healy wrote: > Yet I try gsl_sf_ellint_E_e > for phi = 0.5pi and k=0.5, I find the value returned is 1.46746 (which > agrees with the result from the complete elliptic integral as it > should), > and for phi= 0.6pi and k=0.5, value is 1.19394. In fact, for phi=pi, > I get essentially 0, when it looks like I should get 2*1.46746 because > sin^2 is symmetric about pi/2.
I tried what you decribed and got results that seem to be correct also for phi=0.6*pi and phi=pi. Could you send a complete (but short) program that shows the behavior you mentioned? Thanks, Frank _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
