[sage-support] Additionnal constraints are never enough

2011-03-07 Thread Maxim
Hello, I'd like to know if it's possible to get Sage to process the
following integrals:

reset()
var('x y z')
assume(z-16)
assume(z0)
integrate(integrate(1/16,y,-2,(z+2*x+4)/(x+2)),x,(-z-8)/4,2)

Sage always return: Is  z+16  positive, negative, or zero?

Maxim

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


Re: [sage-support] Additionnal constraints are never enough

2011-03-07 Thread D. S. McNeil
 integrate(integrate(1/16,y,-2,(z+2*x+4)/(x+2)),x,(-z-8)/4,2)

Seems to work for me:

--
| Sage Version 4.6.1, Release Date: 2011-01-11   |
| Type notebook() for the GUI, and license() for information.|
--
sage: var(x y z)
(x, y, z)
sage: integrate(integrate(1/16,y,-2,(z+2*x+4)/(x+2)),x,(-z-8)/4,2)
[...]
TypeError: Computation failed since Maxima requested additional
constraints (try the command 'assume(z+160)' before integral or limit
evaluation, for example):
Is  z+16  positive, negative, or zero?
sage: assume(z+160)
sage: integrate(integrate(1/16,y,-2,(z+2*x+4)/(x+2)),x,(-z-8)/4,2)
[...]
TypeError: Computation failed since Maxima requested additional
constraints (try the command 'assume(z0)' before integral or limit
evaluation, for example):
Is  z  positive, negative, or zero?
sage: assume(z0)
sage: integrate(integrate(1/16,y,-2,(z+2*x+4)/(x+2)),x,(-z-8)/4,2)
-1/16*z*log(-1/4*z) + 1/16*z*log(4) + 1/16*z + 1

What version of Sage are you running?  (Note that 4.6.2 has just been released.)


Doug

--
Department of Earth Sciences
University of Hong Kong

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org