function F_c = F_c(t,y)
diff = (y(1)-eta_dot(t));
if (diff < 0)
   F_c = 0;
else
   F_c = -100*diff; 
end
end

function eta = eta(t)
eta_0 =2.0;  %1 meter wave amplitude
T = 10;  %10 second period 
eta = eta_0*cos(2*pi*t/T);
end

function eta_dot = eta_dot(t)
eta_0 =2.0;  %1 meter wave amplitude
T = 10;  %10 second period 
eta_dot = -eta_0*2*pi*sin(2*pi*t/T)/T;
end