delip - elliptic integral
returns the value of the elliptic integral with parameter ck :
x real and positive. When called with x a real vector r is evaluated for each entry of x.
ck=0.5; delip([1,2],ck) deff('y=f(t)','y=1/sqrt((1-t^2)*(1-ck^2*t^2))') intg(0,1,f) //OK since real solution!