On Wednesday, June 25, 2014 12:39:30 PM UTC+8, P Purkayastha wrote:
>
> On Tuesday, June 24, 2014 8:51:18 PM UTC+8, kcrisman wrote:
>>
>>
>>
>> You, saw a plot? I didn't see a plot.
>>>
>>>
>>
>> Hmm. Can you say more about *exactly* what version of Sage, setup, and
>> so forth you are using? The plot appears for me with Sage 5.12 but not
>> in Sage 6.3.beta3 (though that one may have had #13246 applied). ppurka,
>> do you have a sense as to whether the stuff at #13246 would be responsible
>> for this?
>>
>> Thanks for reporting this - we'll have to track it down a bit more.
>>
>
> It "works" with 6.3.beta4. And it fails with #13246 applied to 6.3.beta4.
>
> But there is a reason why it fails with #13246 applied and "works" without
> it. The plot has only three data points. So, #13246 is doing its job just
> as it was intended to I suppose :-)
>
> sage: var('n,a')
> sage: hbar, m = 1,1
> sage: psi(x,t,n)=sqrt(2/a)*sin(n*pi*x/a)*e^(-i*n^2*pi^2*hbar*t/(2*m*a^2));
> sage: print psi
> sage: Psi(x,t)=1/sqrt(2)*psi(x,t,1)+1/sqrt(2)*psi(x,t,2);
> sage: print Psi
> sage: P(x,t,a) = Psi.conjugate()*Psi
> sage: print P.expand()
> sage: p = plot(P(x,1,1),x,0,1)
> sage: po = p._objects[0]; po.xdata
> [0.0,0.5,1.0]
>
>
Here is a way to fix this problem. Since you know that the output of P() is
always real, simply use the following:
plot(P(x,1,1).real(), x, 0, 1)
And you will get a plot like the one attached.
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To post to this group, send email to sage-support@googlegroups.com.
Visit this group at http://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/d/optout.
