I had the code based on Bill Page and Waldek Hebisch suggestions that was doing what I wanted but this last hint by Themos simplified it a lot. I will put it here as Marduk seems to be doing a similar exercise in the "Simple Gaussian integral fails thread".
The goal is to see if we can make FriCAS compute the integral that gives the Black-Scholes formula for the price of a call option. So, we define C:=exp(-r*T)/(s*sqrt(2*%pi*T)) f(x)==(S*exp(x)-K) g(x) == exp(-(x-T*m)^2/(2*T*s^2)) m := r-s^2/2 The integral integrate(f(x)*g(x),x=log(K/S)..%plusInfinity) fails, but we can do I := C*integrate(f(x)*g(x),x=log(K/S)..y) Now we need to take the limit as y goes to infinity. For that to work we first have to tell FriCAS that the expression sqrt(1/(2*T))/s is positive. We do it by replacing this expression by c^2. withPos:=subst(I,[sqrt(1/(2*T))::Kernel(Expression(Integer))],[c^2*s]) Nowe we can take the limit: lim:=limit(withPos,y=%plusInfinity) and replace c^2 back by sqrt(1/(2*T))/s BS := (rule ('c^2==sqrt(1/(2*T))/s)) (lim::Expression(Integer)) -- 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.