function f=func(t,y)
f=1-y;

% Example for system of ODEs
%function dy=func(t,y)
% omega=1.0;
% dy(1,1) = y(2);
% dy(2,1) = -omega^2*y(1);