In [1]:
#definimos nuestras variables
x,u,v,w = var("x,u,v,w")
eqn1 = u + cos(x)*v + sin(x)*w == 0
eqn2 = -sin(x)*v + cos(x)*w == 0
eqn3 = -cos(x)*v - sin(x)*w == sec(x)
sol = solve([eqn1, eqn2, eqn3], [u,v,w], solution_dict=True)
print sol
[{u: sec(x), v: -cos(x)*sec(x)/(cos(x)^2 + sin(x)^2), w: -sec(x)*sin(x)/(cos(x)^2 + sin(x)^2)}]
In [2]:
print sol[0][u]
sec(x)
In [3]:
u(x)=sol[0][u]
show(u)
In [4]:
U(x)=u.integrate(x)
show(U)
In [5]:
v(x)=sol[0][v]
show(v)
In [6]:
v(x)=v.simplify_full()
show(v)
In [7]:
V(x)=v.integrate(x)
show(V)
In [8]:
w(x)=sol[0][w]
show(w)
In [9]:
w(x)=w.simplify_full()
show(w)
In [10]:
W(x)=w.integrate(x)
show(W)
In [ ]: