For the first-order differential equation: (Please solve all parts and provide MATLAB code, Thank You)

function euler_for.m is provided here:

function [x,y] = euler_for(x0,y0,x_end,h,fcn) % % function [x,y]=euler_for(x0,y0,x_end,h,fcn) % % Solve the initial value problem % y' = f(x,y), x0 x_end-h % The vector y will contain the estimates of the solution Y % at the node points in x. % n = fix((x_end-x0)/h)+1; x = linspace(x0,x_end,n)'; y = zeros(n,1); y(1) = y0; for i = 2:n y(i)=y(i-1)+h*feval(fcn,x(i-1),y(i-1)); end