In [1]:
%pylab inline
In [2]:
def transmonFreq(x, b):
return sqrt(8 * b[3] * (b[0] * cos(pi * x + b[2])*sqrt(1 + b[1]**2 * tan(pi * x + b[2])**2)))
In [8]:
x = linspace(-.85, .85, 101) * .6
y = transmonFreq(x, [15, 1/3, 0, .226])
plt.plot(x, y)
Out[8]:
In [ ]: