function f=rhslin(t,y)
global a b c d
n=length(y)/2;
y1=y(1:n);       % Vector of n concurrently integrated x values
y2=y(n+1:2*n);   % Vector of n concurrently integrated y values
f=[a*y1+b*y2;c*y1+d*y2];