poleVault.m~

%Program used to evaluate question 6

%Define 
N= [5,6,10]-1;
P=10;
E=76e9;
I=4e-9;
L=5;

i=1;
for n=N
    Dx=L/n
    p=P/(E*I);
    a=Dx^2*p^2-2;
    m=diag(ones(1,n-1),-1)+diag(a*ones(1,n))+ diag(ones(1,n-1),1)
    e=eig(m);
    i;
    max_lambda(i)= max(e);
    min_lambda(i)=min(e);
    i=i+1
end

max_lambda'
min_lambda'
N'
nodes=N+1;
nodes'