%right hand side of the equation Y'=-Y+2
function dy=hw4_1_fun(t,Y)
dy=zeros(1,1); %makes output a 1d vector
dy(1)=-Y(1)+2;