Experimenting with Numerical Python

example1.m

% Matlab code to plot a cavity resonance

F = 10; % Finesse

f = linspace(-0.5, 1.5, 201);
P = 1./(1 + (2/pi) * F^2 * sin(pi*f).^2);

plot(f, P);
xlabel('free spectral ranges');
ylabel('power buildup');

example1.py

#!/usr/bin/python
# Python code to plot a cavity resonance 

import numpy as numpy
F = 10
f = numpy.linspace(-0.5, 1.5, 201)
P = 1 / (1 + (2/numpy.pi) * F**2 * numpy.sin(numpy.pi * f) ** 2)

import matplotlib.pyplot as plt
plt.plot(f, P)
plt.xlabel("free spectral ranges")
plt.ylabel("power buildup")
plt.show()

example1b.py

#!/usr/bin/python

from numpy import *
from matplotlib.pyplot import *
F = 10
f = linspace(-0.5, 1.5, 201)
P = 1/ (1 + (2/pi) * F**2 * sin(pi*f)**2)
plot(f, P)
xlabel("free spectral ranges")
ylabel("power buildup")
show()

python_vs_matlab.png