Eigen_Values.m

N= [5,6,10]-1;

%Given constant values
E=76e9;
I=4e-9;
L=5;

%For loop that creates:
    %1 - the lambda values 
    %2 - number of segments
    %3 - number of eignvalues
i=1;
for m=N
    Dx=L/m;
    p=10/(E*I);
    a=Dx^2*p^2-2;
    m=diag(ones(1,m-1),-1)+diag(a*ones(1,m))+ diag(ones(1,m-1),1)
    e=eig(m);
    i;
    maxLam(i)= max(e);
    minLam(i)=min(e);
    i=i+1;
end

maxLam'
minLam'
N'
nodes = N+1;
nodes'