bosa na kokosa

function project55

x0 = 1; y0 = 2; a = 3;

sol = dsolve('Dy = (2*y)/x + 4*(x^4)', 'y(1) = 2', 'x');
x = linspace(x0,a);
plot(x, eval(sol), 'k')

grid on
hold on

h = [0.5, 0.4, 0.02];
c = ['b', 'g', 'm'];

function z = f(x,y)
    z = (2*y)/x + 4*(x^4);
end

for j = 1: length(h)
    x = []; y = [];
    x = x0:h(j):a;
    y(1) = y0;
    for k = 1: length(x)-1
        y(k+1) = y(k) + h(j)*f(x(k), y(k));
    end
    plot(x,y,c(j))
end

end