%1A)
%Skapar givna värde
mh=90;
mr=10;
ma=10;
kh=8950000;
ka=10^6;

%Skapar de givna matriser M och K
M=[(mh+mr),0;0,ma];
K=[(kh+ka),-ka;-ka,ka];

%Beräknar egenvinkelfrekvensen för huset samt absorbatorn
wh=sqrt(kh/(mh+mr))
wa=sqrt(ka/ma)

%Beräknar egenvärde lamda för att sedan få egenvinkelfrekvensen
w=eig((inv(M)*K));
wanswer=sqrt(w)

[x,y]=eig((inv(M)*K))

%1B)
mh=90;
mr=10;
ma=10;
kh=8950000;
w=pi*100;
wa=w*[0.3:0.01:1.5];
M=[(mh+mr),0;0,ma];

%Går igenom alla värden i wa för att spara admittansvärdena i
% vektorerna A11 och A12 för att sedan kunna plotta dessa.
for j=1:length(wa)
    ka=wa(j)^2*ma;
    K=[(kh+ka),-ka;-ka,ka];
    A=inv((-w^2*M)+K);
    A11(j)=A(1,1);
    A12(j)=A(1,2);
end

plot(wa,A11,'r')
hold on
plot(wa,A12,'g')
legend('A11', 'A12')
xlabel('vinkelhastighet, rad/s')
ylabel('Admittans')
title('Admittanskomponenterna')

A110=find(A11==0);
wa(A110)