brenogil's profile - activity

I was trying to convert this into sage. Basically print function, first derivative, second derivative all sharing the x axis. Instead of scipy.misc.derivative I'd like to use sage derivative(). I have tries for hours an got several different errors. Any ideas?

d={}
with open(source) as f:
    for line in f:
       (key, val) = line.split()
       d[str(key)] = val
       print key
       f = lambda t : eval(val)
       x = np.arange(1,20,0.01)
       first = derivative(f,x,dx=1,n=1)
       second = derivative(f,x,dx=1,n=2)
       fig,ax = plt.subplots(3,1,sharex=True)
       #plt.suptitle("\n".join(wrap(str(key + " - \n" + val), 50)))
       ax[0].plot(x,f(x))
       ax[0].set_ylabel(r'$f(x)$')
       ax[1].plot(x,first)
       ax[1].set_ylabel(r'$f\/\prime(x)$')
       ax[2].plot(x,second)
       ax[2].set_ylabel(r'$f\/\prime\prime(x)$')
       ax[2].set_xlabel(r'$x$')
       plt.savefig()
       plt.close()

The results should be something like this (i would upload my own image but I realized I cannot): http://matplotlib.org/_images/subplot...

Deleted that because it was not really an answer. I just meant that even facing similar results, sage does seem to behave differently when trying to solve our problems.

I guess your case is a bit different because even p0 and p2 being sage.symbolic.expression.Expression

bool(p0+p2-1<0)still yields a True

This is quite close to what I was expecting. Thank you!

Thank you very much!

The following positions the labels at the ends of the axes.

plot(E,t,(x,0,15),axes_labels=['$t$ axis','$E$ axis'])

You could use the text command to position a label.

t = var('t')
E(t)=(10.40*(1-exp(-0.25*t)))^2.81    
p=plot(E,t,(x,0,15),axes_labels=['','$E$ axis'])
p+=text('$t$ axis',(8,-30),color='black')
show(p)

You may also be able to use matplotlib to get what you want.

Or this may be what you want:

t = var('t')
E(t)=(10.40*(1-exp(-0.25*t)))^2.81    
plot(E,t,(x,0,15),frame=True,axes_labels=['$t$ axis','$E$ axis'],axes=False)

Hello,

How do I add labels (names) to my axes when doing a simple 2d plot such as:

t = var('t')
E(t)=(10.40*(1-exp(-0.25*t)))^2.81    
plot(E,t,(x,0,15))

I would like to do so the labels are aligned directly under the x axis and vertically over the y axis. Is there a way?

Thanks for the hand! I'll pay you a beer whenever I am in Canada!

Thank you! I do no know if I am missing some package or something but I cannot make the graph 'behave' like the picture above.

Hello,

I have this marix:

M=matrix(ZZ,6,9,[[1,1,1,0,0,0,0,0,0],[0,0,0,1,1,1,0,0,0],[0,0,0,0,0,0,1,1,1],[-1,0,0,-1,0,0,-1,0,0],[0,-1,0,0,-1,0,0,-1,0],[0,0,-1,0,0,-1,0,0,-1]])

and I am looking to get something like this:

How can I do it?

How can this be?

sage: x = var('x')
sage: Q = function('Q', x)
sage: assume(x>0)
sage: assume(Q>0)
sage: bool(Q>0)
False

And even

sage: bool(abs(Q)>0)
False

Am I forgetting or neglecting something? Thank you

I also have this doubt

I also have this doubt