[GiNaC-list] Crash Ginsh when trying to solve equation involving integral numerically
Alexei Sheplyakov
alexei.sheplyakov at gmail.com
Sat Aug 21 18:51:15 CEST 2010
Hello,
On Mon, Aug 16, 2010 at 10:09:01PM +0200, Ernst Moritz Hahn wrote:
> I tried to solve an equation numerically, using Ginsh from the git
> repository. When entering
>
> fsolve((1/(sqrt(2*Pi)))*integral(t,0,x,exp(-1/2*t^2))==0.5,x,0,100)
>
> into Ginsh, it crashes, see below. I got the same problem with the
> version shipped with the latest stable Ubuntu Linux.
I've posted patches which should fix this bug to ginac-devel (the combined
patch is attached for convenience). Could you please check if they work for
you?
Best regards,
Alexei
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fsolve_fix_allin1.diff
Type: text/x-diff
Size: 3431 bytes
Desc: not available
URL: <http://www.cebix.net/pipermail/ginac-list/attachments/20100821/bff525c6/attachment.bin>
More information about the GiNaC-list
mailing list