function [out1] = LINfunction(t,y,flag,P)
% ----------------------------------------------------------------------
%   y(1) ... x
%   y(2) ... y
out1(1)= P.A(1,1)*y(1) + P.A(1,2)*y(2); 
out1(2)= P.A(2,1)*y(1) + P.A(2,2)*y(2); 
out1= out1';    % wants output as a column vector