sophia's profile - activity

that works, thank you!

I have the following code:

var('a')
_.<k> = PolynomialRing(ZZ)
f = k^3+2*k^2+1
g = k^3 + a*k^2 + 1
f.coefficients()
g.coefficients()

the coefficients of f that i get are perfect:

[1, 2, 1]

the coefficients of g should be [1,a,1], but instead i get:

 [[k^3 + 1, 0], [k^2, 1]]

I've tried the solution suggested here:

but it doesn't always produce the coefficients in a logical order.

As always, any help figuring out how to get the correct coefficients of g would be greatly appreciated.

Thank you!

I have a formatting query about the background colors on sagemath cloud. I regularly generate and print pdf of the sage worksheets on a black and white printer. On the printed page, its hard to distinguish between the cell that contains the code and the output of the cell (both have the same background color).

Is there any way to clearly demarcate the cell that contains the code and the output of the cell? Would it be possible to get the output of any cell with a shaded background or with a border around it, so that on a printed paper, the output of the cell appears clearly different from the contents of the cell.

I've tried some of the Editor color scheme options, but I don't see an option where the cell appears with a different background color.

Thank you.

Thank you, I am able to access the latest version of sympy now. My original question was on solving the irrational inequality above. Is that not possible in sagemath? Thanks.

Is there any way to solve the following inequality for a in Sage?

(It's pretty easy to do in Maple.)

abs(-1/4*a - 1/4*sqrt(a^2 - 30*a + 17) + 1/4) < 1

Using solve doesn't work. I tried using the QEPCAD package, but on SageMathCloud, i get an error: "unable to start QEPCAD".

I tried using sympy-solve and sympy-solveset, and that didn't work either. I plan to raise questions regarding sympy on stackexchange, but if anyone has helpful guidelines on using sagemath's solve vs sympy-solve, that'd be greatly appreciated. My general experience has been that sympy-solve is a lot more capable than sagemath's solve.

Any help would be greatly appreciated. Thank you.

EDIT

Thank you for replying. I am able to use qepcad now, but the problem I'm trying to solve takes a very long time. Here is the code:

var('a')
dnf = solve(abs(-1/4*a - 1/4*sqrt(a^2 - 30*a + 17) + 1/4) <= 1, a)
qf = apply(qepcad_formula.or_, map(qepcad_formula.and_, dnf)) # reformat the solution
qf
qepcad(qf, vars='(a)') # simplify

Any help solving this would be greatly appreciated. Also, would it be possible to update the version of sympy on sagemathcloud? Thank you.

that works great, thank you!

your example is correct. I'm trying to accomplish something like this:

the rhs of de1 is just y(t). lets call the rhs of de1 F, so F = y(t). now dF/dy = 1. but dF/dx = 0, and that's what i'm trying to calculate - the derivative of the rhs of de1 w.r.t y and w.r.t x.

Thanks a lot for replying. I tried what you suggested with de1 first, and it works when the rhs is only a function of y(t). so to compute dF/dx, i calculated it as (dF/dt)/(dx/dt) as follows:

sage: de1 = diff(x(t),t) == y(t) sage: F = de1.rhs() sage: (diff(F,t)/diff(x(t),t)).subs(de1)

the result is : D0(t)/y(t), whereas it should just be Zero.

any further suggestions would be greatly appreciated. Thank you.

Thanks for replying. Please do let me know if you find a solution to this.

Dr. Stein explains in the link below that FEniCS can be accessed on sagemath cloud. I followed the instructions in his post but i get an error that says: "No module named dolfin".

http://ask.sagemath.org/question/1100...

I'd appreciate any help figuring out how to import dolfin in sagemath.

Thanks,

that helps me solve what i wanted to, thank you!

i had simplified the problem that i wanted to solve considerably. finding the Jacobian will require solving for the derivative of a functional w.r.t a function. Instead of crowding this post, i've asked a new question here:

http://ask.sagemath.org/question/32876/finding-the-derivative-of-a-functional-wrt-a-function/ (http://ask.sagemath.org/question/3287...)

thank you again for your help.

I've defined a system as:

var('a, t'); function('x, y');

de1 = diff(x(t),t) == y(t); de2 = diff(y(t),t) == -ax(t) - (a-4)/ay(t) - y(t)^3;

I'd like to compute the derivatives of the rhs of de1 and de2 w.r.t to x(t) and y(t) (to eventually form a Jacobian matrix), i.e. the derivative of the rhs of de1 w.r.t to x(t) is Zero and w.r.t y(t) is 1. I've tried the following:

diff(de1.rhs(),y); diff(de1.rhs(),y(t)); derivative(de1.rhs(),y(t));

I get errors on all three. I'd appreciate any help. thank you.

thanks for replying. i don't want to use desolve yet. all i'd like to do is to solve the equations algebraically (explicitly) for the derivatives. thank you,

thanks for replying. i'd like to automate the solution process. i tried substituting x1 for diff(x(t),t), but got an error ( i don't want to manually type in the x1 and y1 part). i tried both:

sage: de1.substitute_function(diff(x(t),t),x1) and sage: de1.subs(diff(x(t),t)=x1)

any suggestions as to how i can sub x1 for diff(x(t),t) directly in de1?

thank you,

I'm not trying to solve the differential equations just yet. i only want to solve the equations algebraically (explicitly) for the derivatives, so I can compute the Jacobian matrix of the right hand sides for further bifurcation analysis. the equations in my actual problem have parameters.

(this is pretty easy to do in Maple)

I have the following defined:

sage: var('a, t');

sage: function('x, y, v, u')

sage: de1 = diff(x(t),t) - diff(y(t),t) == x(t)

sage: de2 = diff(y(t),t) == y(t)

I'd like to solve these two equations algebraically for diff(x(t),t) and diff(y(t),t).

sage: sol = solve(de1, diff(x(t),t)); works, and i get an explicit solution for diff(x(t),t), but when i try solving both equations using:

sage: sol = solve([de1, de2],diff(x(t),t),diff(y(t),t));

i get the following error:

Traceback (most recent call last): File "/projects/sage/sage-6.10/local/lib/python2.7/site-packages/smc_sagews/sage_server.py", line 905, in execute exec compile(block+'\n', '', 'single') in namespace, locals File "", line 1, in <module> File "/projects/sage/sage-6.10/local/lib/python2.7/site-packages/sage/symbolic/relation.py", line 822, in solve variables = tuple(args[0]) TypeError: 'sage.symbolic.expression.Expression' object is not iterable

=============================

I'd appreciate any help with this. thank you.