function [dydt] = prob_7_6_fn(t,y)
% 
% Written by:
% -- 
% John L. Weatherwax                2006-07-31
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

% 
% 
dydt = -1000*(y - sin(t)) + cos(t);