Eulers metod

Uppgiften handlar om att implementera Eulers metod i MATLAB och lösa ekvationen 

$y\text{'}-e^{-x}\cos \left(y\right)=\sin \left(2x\right) , y\left(0\right)=1$

Jag har en en mall som jag har börjat fylla i men jag vet inte hur jag ska fortsätta med själva lösningen

a=0; b=4; 			% [a,b] intervallet där lösning beräknas

ya=1; 				% begynnelsevärdet, y värdet i vänster ändpunkt

f=@(x,y)exp(-x)*cos(y)+sin(2*x); % funktionen som berskriver differentialekvationen

M=4; 				% antalet indelningspunkter längs x-axeln

x=linspace(a,b,M);		 % indelningspunkterna på x-axeln

y=euler(x,ya,f); 		% anropar lösaren

plot(x,y)		 	% plottar lösningen

function y=euler(x,y0,f)

y=zeros(size(x)); 	% y en vektor med nollor, lika lång som x vektorn.
                  	% i denna vektor ska y-värdena läggas in, ett och ett
              	    	% det blir en loop
y(1)=y0; 		% lägg in begynnelsevärdet på första position.

    % Här ska jag alltså skriva nån kod som löser ekvationen
    
    
   


    
end

Efter lite googlande hittar jag en massa olika sätt att utnyttja Eulers formel och tror jag förstår hur Eulers formel fungerar men inte hur jag ska kunna skriva in det i den här mallen , och jag vill inte gärna kopiera en massa kod som jag inte förstår vad den faktiskt gör