I found problem integrating many problems using sagemath calling Fricas to 
do the integration when using exponential integral function. These are 
problems from Rubi test files.

It works OK with other CAS systems supported by sagemath (Maxima and GIAC) 
but sagemath 10.3 does not seem to correctly translate the call to Fricas.

I am having hard time finding what the exponential integral function is 
called before I ask at sagemath forum. I looked at the Fricas book and do 
not see anything,. I tried Ei but this does not work. (i.e. does not give 
same answer as other cas systems).

First here is a link to the special function I am taking about

https://reference.wolfram.com/language/ref/ExpIntegralE.html

https://www.maplesoft.com/support/help/maple/view.aspx?path=Ei

Here is a test to what values it should give for some random input. In 
Mathematica it gives 

ExpIntegralE[3, 5.0]
.000877801

In Maple

Ei(3,5.0)
0.0008778008928

IN sagemath 10.3

sage: exp_integral_e(3,5.0)
0.000877800892770638

But I tried Ei(3,5.0)   in Fricas and it gives error.

Here is an example, using sagemath trying to integrate. It works OK with 
maxima and giac but gives error with Fricas. I am using 

>fricas --version
FriCAS 1.3.10
based on sbcl 2.3.11

with sagemath 
>sage --version
SageMath version 10.3, Release Date: 2024-03-19

Starting sagemath and typing:

sage: var('x a b')
(x, a, b)

sage: integrate(exp_integral_e(1,b*x),x,algorithm="giac")
integrate(exp_integral_e(1, b*x), x)

sage: integrate(exp_integral_e(1,b*x),x,algorithm="maxima")
-exp_integral_e(2, b*x)/b

sage: integrate(exp_integral_e(1,b*x),x,algorithm="fricas")
RuntimeError                              Traceback (most recent call last)
TypeError: An error occurred when FriCAS evaluated 
'exp_integral_e(((1)::EXPR INT),(b)*(x))':
   There are no library operations named exp_integral_e

So clearly sagemath did not translate the exp_integral_e to Fricas 
correctly.

What should the translation look like?  

Thanks
--Nasser 



-- 
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 view this discussion on the web visit 
https://groups.google.com/d/msgid/fricas-devel/2db2efb2-2638-4a84-8739-f14f3d3dfe0an%40googlegroups.com.
  • [fricas-devel] wha... 'Nasser M. Abbasi' via FriCAS - computer algebra system

Reply via email to