function dydt=linkder(t,y)
dydt=zeros(3,1);
A=[1 -2 4; -1 1 -1; 1 1 1]; % a rendszer mátrixa
dydt= A*y;