OK, thanks! tirsdag 3. januar 2017 14.07.18 UTC+1 skrev Dima Pasechnik følgende: > > if you remove 'adaptive=True' then it works---so probably you're hitting a > bug in computing the best box to plot, or something like this. > > On Tuesday, January 3, 2017 at 12:39:05 PM UTC, Fjordforsk A/S wrote: >> >> I compared Wolfram Alpha online and SAGE just to see, and it gives a >> strange result: >> >> Wolfram plots the function as given : >> >> >> https://www.wolframalpha.com/input/?i=plot+(1%2Bx*(e%5E(x%2Bt%2B1)%2Be%5E(2*x%2B2*t%2B2)-(x%5E2*(I*x-t)%2F(2x%2Bt%2B2))*e%5E(3*x%2B3*t%2B3)) >> >> >> while SAGE with command: plot3d(lambda >> x,t:(1+x*(e^(x+t+1)+e^(2*x+2*t+2)-(x^2*(I*x-t)/(x+t+2))*e^(3*x+3*t+3))).real_part(),(x,-3*pi,3*pi),(t,-4,4),adaptive=True,color=cmsel) >> >> gives: >> >> ZeroDivisionError: Symbolic division by zero >> >> >> That is a logically strange result. May one of the two be wrong? >> >> >> >> tirsdag 3. januar 2017 13.24.58 UTC+1 skrev Dima Pasechnik følgende: >>> >>> >>> >>> On Tuesday, January 3, 2017 at 12:02:41 PM UTC, Fjordforsk A/S wrote: >>>> >>>> Hello, I get an idle response from SAGE when trying to plot this: >>>> >>>> var('x t'); >>>> cmsel = [colormaps['gnuplot2'](i) for i in sxrange(0,1,0.02)] >>>> plot3d(lambda >>>> x,t:(1+x*(e^(x+t+1)+e^(2*x+2*t+2)-(x^2*(I*x^2-t)/(x^3+t^2+2))*e^(3*x+3*t+3)).imag_part(),(x,-3*pi,3*pi),(t,-2,2),adaptive=True,color=cmsel) >>>> >>> >>> one missing ')' before .imag_part() >>> >>> Then it works for me. >>> >>
-- 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 https://groups.google.com/group/sage-support. For more options, visit https://groups.google.com/d/optout.
