I am in the process of learning Sage, coming from Maxima (and Mathematica,
which I do not like much...).
Cut'n'pastes from a notebook running on sagenb.org
version()
version()
==>
'Sage Version 5.4, Release Date: 2012-11-09'
var('t,a,b,d')
## beta density
dbeta(t,a,b)=t^(a-1)*(1-t)^(b-1)/beta(a,b)
## density of the *difference* of two independent beta(1,1)-distributed RVs
## (yes, I mean uniforms(0 1)...)
integrate(dbeta(t,1,1)*dbeta(t+d,1,1),t,max(0, -d), min(1,1-d))
==>
1
Huh ?? This is seriously whacky :
dbeta(t,a,b):=t^(a-1)*(1-t)^(b-1)/beta(a,b)
integrate(dbeta(t,1,1)*dbeta(t+d,1,1),t,max(0,-d),min(1,1-d))
==>
dbeta(t,a,b):=t^(a-1)*(1-t)^(b-1)/beta(a,b)
min(1,1-d)-max(0,-d)
The problem appears to lie with integrate() :
integrate(1,t,max(0,-d),min(1,1-d))
==>
1
integrate(t,t,max(0,-d),min(1,1-d))
==>
1/2
But another one : i tried to give dbeta a proper definition, i. e. with a
domain of definition, therefore allowing convolutions :
db2(t,a,b)=Piecewise([[(0,1),t^(a-1)*(1-t)^(b-1)/beta(a,b)]])
db2
==>
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "_sage_input_7.py", line 10, in <module>
exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n"
+
_support_.preparse_worksheet_cell(base64.b64decode("ZGIyKHQsYSxiKT1QaWVjZXdpc2UoW1soMCwxKSx0XihhLTEpKigxLXQpXihiLTEpL2JldGEoYSxiKV1dKQpkYjI="),globals())+"\\n");
execfile(os.path.abspath("___code___.py"))
File "", line 1, in <module>
File "/tmp/tmpzOS0Jf/___code___.py", line 3, in <module>
__tmp__=var("t,a,b"); db2 = symbolic_expression(Piecewise([[(_sage_const_0
,_sage_const_1 ),t**(a-_sage_const_1 )*(_sage_const_1 -t)**(b-_sage_const_1
)/beta(a,b)]])).function(t,a,b)
File
"/sagenb/sage_install/sage-5.4-sage.math.washington.edu-x86_64-Linux/local/lib/python2.7/site-packages/sage/calculus/all.py",
line 96, in symbolic_expression
return SR(x)
File "parent.pyx", line 804, in sage.structure.parent.Parent.__call__
(sage/structure/parent.c:7228)
File "coerce_maps.pyx", line 82, in
sage.structure.coerce_maps.DefaultConvertMap_unique._call_
(sage/structure/coerce_maps.c:3547)
File "coerce_maps.pyx", line 77, in
sage.structure.coerce_maps.DefaultConvertMap_unique._call_
(sage/structure/coerce_maps.c:3449)
File "ring.pyx", line 284, in
sage.symbolic.ring.SymbolicRing._element_constructor_
(sage/symbolic/ring.cpp:4846)
TypeError
Threfore, whackiness might lie with a more general problem.
Where am I goofing ?
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
Visit this group at http://groups.google.com/group/sage-support?hl=en.
