function y=perakkain(t,x)
global K1 K2
y=[-1 0 0 0;1 -K1 0 0;0 K1 -K2 0;0 0 K2 0]*x;
end