function y=u(t);
% file name: u.m
% create a cosine with amplitude u0
% this function is used by linear.m and reactor.m
global u0;
b=1;
y=zeros(2,1);  % create a column vector
y=u0*cos(b*t);

