function yp = lake2(t,y)

% Model for Lake Pollution
% Parameters for the model 
f = 0.001+0.0006*sin(0.0172*t);
yp = -f*(y - 8*exp(-0.002*t));
end

% [t,y] = ode23(@lake2,[0,500],0);
