1 | initial version |
So this is what I would have wanted to do.. prior to learning about it from here. :)
var('p, alpha, beta, q, A, B, J, K') h=solve([p==alpha+2betaq, A-Bp==Jq, pq-K-alphaq-beta*q^2==0],p,q,J) show(h[1][0].simplify_full()) show(h[1][1].simplify_full()) show(h[1][2].simplify_full()) show(h[1][0].rhs().diff(A).simplify_full()) show(h[1][1].rhs().diff(A).simplify_full()) show(h[1][2].rhs().diff(A).simplify_full())
2 | No.2 Revision |
So this is what I would have wanted to do.. prior to learning about it from here. :)
var('p, alpha, beta, q, A, B, J, K')
h=solve([p==alpha+2betaq, q,
A-Bp==Jq, q,
pq-K-alphaq-beta*q^2==0],p,q,J)
show(h[1][0].simplify_full())
show(h[1][1].simplify_full())
show(h[1][2].simplify_full())
show(h[1][0].rhs().diff(A).simplify_full())
show(h[1][1].rhs().diff(A).simplify_full())
show(h[1][2].rhs().diff(A).simplify_full())
3 | No.3 Revision |
So this is what I would have wanted to do.. prior to learning about it from here. :)
var('p, alpha, beta, q, A, B, J, K')
4 | No.4 Revision |
So this is what I would have wanted to do.. prior to learning about it from here. :)
var('p, alpha, beta, q, A, B, J, K')
h=solve([p==alpha+2*beta*q, A-B*p==J*q, p*q-K-alpha*q-beta*q^2==0],p,q,J)
show(h[1][0].simplify_full())
show(h[1][1].simplify_full())
show(h[1][2].simplify_full())
show(h[1][0].rhs().diff(A).simplify_full())
show(h[1][1].rhs().diff(A).simplify_full())
show(h[1][2].rhs().diff(A).simplify_full())
and this..
h=solve([p==alpha+2*beta*q, A-B*p==J*q, p*q-K-alpha*q-beta*q^2==0],p,q,J,solution_dict=True)
show([[s.diff(A).simplify_full() for s in sol.values()] for sol in h])