clc, clear
plot_analytic=1;    %1: plot
%blue, green, red, cyan, magenta, yellow, black (k)

a=0;
b=4;
ya=1;
n=8;

%f = @(x,y) 4*exp(0.8*x)-0.5*y;
f = @(x,y) -2*x^3+12*x^2-20*x+8.5;

[x,y]=heun(f,a,b,ya,n);
if n<20
    E=[x' y']
end
plot(x,y,'gs-','LineWidth',3);
grid on
hold on %Ard arda çizdirilen grafiklerin üst üste konumlandırılması için
% Analitik çözümle karşılaştırma için 
if plot_analytic==1
    x_a = linspace(a,b,40);
    y_a = -0.5*x_a.^4 + 4*x_a.^3 - 10*x_a.^2 + 8.5*x_a + 1;
    plot(x_a,y_a,'r-','LineWidth',3);
    hold on
end