That looks correct. Sure wish I could do that :D
Chris
On Tue, Jun 24, 2014 at 9:50 PM, P Purkayastha <ppu...@gmail.com> wrote:
>
>
> 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.
>
--
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.
