Sorry, I've realized that numquad [1] works on Float and not only on machine doubles, not bad!
Am I missing other quadrature packages (tanh-sinh, gaussian weights)? thanks ric [1] https://fricas.github.io/api/NumericalQuadrature.html?highlight=numquad ---------- (1) -> f(x:Float):Float == { m:=10, n:=11; k:= (4/%pi)^n*n*m; return k* atan(x^m)^(n-1)*x^(m-1)/(1+x^(2*m))} Function declaration f : Float -> Float has been added to workspace. Type: Void (2) -> digits(150);aromberg(f, 0::Float,1::Float,1.e-100,1.e-100,5,20,10) Compiling function f with type Float -> Float (2) [ value = 1.0000000000_0000000000_0000000000_0000000000_0000000000_0000000000_00000 00000_0000000000_0000000000_0000000000_0000000039_1966489901_0326774346_6 923646498_990324404 , error = 0.1832202195_2810277526_5224715964_4965694315_3906127288_2189970717_33840 04174_0332421811_9982065024_0030409954_5097041266_6646499458_3349608857_0 722547214_8324045347 E -100 , totalpts = 3842113, success = true] Type: Record(value: Float,error: Float,totalpts: Integer,success: Boolean) (3) -> -- You received this message because you are subscribed to the Google Groups "FriCAS - computer algebra system" group. To unsubscribe from this group and stop receiving emails from it, send an email to fricas-devel+unsubscr...@googlegroups.com. To post to this group, send email to fricas-devel@googlegroups.com. Visit this group at https://groups.google.com/group/fricas-devel. For more options, visit https://groups.google.com/d/optout.