%LOPASS

%Number of frequencies
Nvals = 500;

%Specify the number of nodes
N = 5;

%initialise all ys
y = zeros(N,N-1,Nvals);
fc= 50.e6;
Zo=50;

frequency = linspace(1.e6,1000e6,Nvals);
omega = 2*pi*frequency;

%specify an L and C in SI units

C1=1.262; L2=1.520; C3=2.239; L4=1.680; C5=2.239; L6=1.520; C7=1.262;

C1 = C1 / (2*pi*fc) /Zo
L2 = Zo*L2 / (2*pi*fc)
C3 = C3 / (2*pi*fc) /Zo
L4 = Zo*L4 / (2*pi*fc)
C5 = C5 / (2*pi*fc) /Zo
L6 = Zo*L6 / (2*pi*fc)
C7 = C7 / (2*pi*fc) /Zo

y(1,2,:) = 1/j./omega/L2; y(3,2,:) = 1/j./omega/L4; y(4,3,:) = 1/j./omega/L6;
y(2,1,:) = 1/j./omega/L2; y(2,3,:) = 1/j./omega/L4; y(3,4,:) = 1/j./omega/L6;

y(1,5,:) = j*omega*C1 + 1/Zo; y(2,5,:) = j*omega*C3; y(3,5,:) = j*omega*C5;
y(5,1,:) = j*omega*C1 + 1/Zo; y(5,2,:) = j*omega*C3; y(5,3,:) = j*omega*C5;

y(4,5,:) = j*omega*C7 + 1/Zo;
y(5,4,:) = j*omega*C7 + 1/Zo;






