import numpy as np
import matplotlib.pyplot as plt
from numpy import sin, cos, tan, pi
from scipy.optimize import fsolve
def eq_circuito(s,q):
A = s[0]
B = s[1]
f = np.empty(2)
f[0] = R*cos(q)+C1*cos(B)+C2*sin(B)-C4+S
f[1] = R*sin(q)-C1*sin(B)+C2*cos(B)-C3
return f
# Constantes:
R = 40
C1 = 75
C2 = 8
C3 = 20
C4 = 100
# Chute inicial:
S0 = np.array([pi/4,30])
n = 100
q = np.linspace(0,pi,n)
Si = np.empty([n,2])
#print(Si)
for k in range(n):
Si[k,:] = fsolve(eq_circuito,S0,q[k])
S0 = Si[k,:]
plt.plot(q,Si[:,0])
plt.figure()
plt.plot(q,Si[:,1])
{"html5":"htmlmixed","css":"css","javascript":"javascript","php":"php","python":"python","ruby":"ruby","lua":"text\/x-lua","bash":"text\/x-sh","go":"go","c":"text\/x-csrc","cpp":"text\/x-c++src","diff":"diff","latex":"stex","sql":"sql","xml":"xml","apl":"apl","asterisk":"asterisk","c_loadrunner":"text\/x-csrc","c_mac":"text\/x-csrc","coffeescript":"text\/x-coffeescript","csharp":"text\/x-csharp","d":"d","ecmascript":"javascript","erlang":"erlang","groovy":"text\/x-groovy","haskell":"text\/x-haskell","haxe":"text\/x-haxe","html4strict":"htmlmixed","java":"text\/x-java","java5":"text\/x-java","jquery":"javascript","mirc":"mirc","mysql":"sql","ocaml":"text\/x-ocaml","pascal":"text\/x-pascal","perl":"perl","perl6":"perl","plsql":"sql","properties":"text\/x-properties","q":"text\/x-q","scala":"scala","scheme":"text\/x-scheme","tcl":"text\/x-tcl","vb":"text\/x-vb","verilog":"text\/x-verilog","yaml":"text\/x-yaml","z80":"text\/x-z80"}