[GiNaC-list] Crash Ginsh when trying to solve equation involving integral numerically

Joerg Arndt arndt at jjj.de
Sun Aug 22 15:29:08 CEST 2010


* Alexei Sheplyakov <alexei.sheplyakov at gmail.com> [Aug 22. 2010 14:57]:
> 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

> [...]

Sehr gut!  Erfolg!  Jawohl!

% ginsh 
ginsh - GiNaC Interactive Shell (ginac V1.5.8)
  __,  _______  Copyright (C) 1999-2010 Johannes Gutenberg University Mainz,
 (__) *       | Germany.  This is free software with ABSOLUTELY NO WARRANTY.
  ._) i N a C | You are welcome to redistribute it under certain conditions.
<-------------' For details type `warranty;'.

Type ?? for a list of help topics.
> fsolve((1/(sqrt(2*Pi)))*integral(t,0,x,exp(-1/2*t^2))==0.5,x,0,100);
6.4971304904725288433
> 


More information about the GiNaC-list mailing list