Ask Your Question
1

integral() failing with "segmentation fault"

asked 2014-01-27 03:37:16 -0500

marv gravatar image

updated 2014-01-27 03:57:58 -0500

tmonteil gravatar image

Hi sage community

My question again arises from exercises in technical chemistry done in sage. I wanted to calculate an definite integral (number of transfer units).

var('x')
ys=0.06+2*(x-0.0275)
xs=ys/1.516

#my function
i=1/(xs/(1+xs)-x/(1+x))
#print i
#show(plot(i,x,0,0.03))

#I, the indefinite integral of i
I=integral(i,x)
#print I
#show(plot(I,x,0,0.03))

#three methods too calculate the integral
print n(I(x=0.0275)-I(x=0))
print numerical_integral(i,0,0.0275)
print integral(i,x,0,0.275)

The first two methods work fine. The last method gives me the infinite of the following messages.

;;;
;;; Detected access to protected memory, also kwown as 'bus or segmentation fault'.
;;; Jumping to the outermost toplevel prompt
;;;

This is a bug, isn't it? However, i have no idea how integrate() works.

Greetings, marv

edit retag flag offensive close merge delete

2 answers

Sort by » oldest newest most voted
2

answered 2014-01-27 03:57:37 -0500

tmonteil gravatar image

updated 2014-01-27 04:08:05 -0500

I can reproduce your behaviour, i confirm this is a bug :(

It is reported as trac ticket 15747.

edit flag offensive delete link more

Comments

thanks for making a track ticket

marv gravatar imagemarv ( 2014-01-27 07:36:46 -0500 )edit
1

answered 2014-01-27 14:46:28 -0500

kcrisman gravatar image

I can confirm that this problem also occurs with the integral to 11/40.

That said,

sage: i
1/((1.31926121372032*x + 0.00329815303430079)/(1.31926121372032*x + 1.00329815303430) - x/(x + 1))

and usually Maxima does not particularly like integrating floating point things exactly, because it sort of doesn't make sense.

sage: numerical_integral(i,0,0.275)
(12.325862793647284, 4.191806193893343e-11)
edit flag offensive delete link more

Comments

you're right, the problem doesn't occur when defining all decimals as fractions

marv gravatar imagemarv ( 2014-01-28 21:31:34 -0500 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 2014-01-27 03:37:16 -0500

Seen: 298 times

Last updated: Jan 27 '14