Algoritmo MATLAB para Euler inverso

Según la fórmula de Euler hacia atrás, podemos obtener:

y(k 1)=y(k)-30*h*y(k 1)

El resultado de la deformación es:

y(k 1)=y(k)/(30*h 1)

Por lo tanto, el programa MATLAB es:

h=0.05;

x=[0:h:1];

y(1)=1;

para k=1: longitud(x)-1

y(k 1)=y(k)/ (30*h 1);

end

plot(x, y, 'r.-');

title('Euler hacia atrás');

cuadrícula encendida