Matlab - Problemas al programar el Met de Euler..

Hola, saludos..

Estoy programando Euler para: dy/dx=x+y; y(0)=1 y me calcula muchos valores y al final reporta unas variables "inf inf nan".
¿Que pasa en el Codigo, porque no converge el Método?
¿Que estoy haciendo mal?
Si alguien me puede informar al respecto se lo agradecería.

Muchas Gracias..!


clear all
clc
xk=0; yk=1; %Condición Inicial
h=0.1;
yk_1=yk;
%% Ciclo Iterativo
tol=10^(-3);
Er=10^(-2);
while Er > tol
yk_1=yk_1+h*(xk+yk_1);
xk=xk+h;
yr=2*exp(xk)-xk-1;
Er=abs(yk_1-yr);
disp([yr yk_1 Er])
end