function dy=simp_solve(t,y)
%


%-----------------
dy = zeros(1,1); 


dy(1)=(-y(1) + 10*sin(100*t))/20;