Texnikes_Ergasia1_Thema1_Zitoumeno3

clear all
clc

% Πίνακας με τα διάφορα l
l = [0.005, 0.01, 0.02];
epsilon = 0.001;
styles = ['o', '*', 'x'];
for index = 1:length(l)
    plotAkrwn(l(index), epsilon, styles(index))
end

function result = plotAkrwn(l, epsilon, style)
    syms x
    f1 = (x-2)^2 - sin(x+3)
    akro1 = 2;  akro2 = 5;
    % Πάντα πρέπει ε < l
    epsilon = 0.001;    counter = 0;
    akra1 = [];    akra2 = [];
    while (akro2 - akro1 > l)
        counter = counter + 1;
        akra1(counter) = akro1;
        akra2(counter) = akro2;
        % Βήμα 1
        x1k = (akro1 + akro2) / 2 - epsilon;
        x2k = (akro1 + akro2) / 2 + epsilon;
        % Βήμα 2
        if subs(f1, x1k) < subs(f1, x2k)
            % ακ+1 = ακ (δεν αλλάζει) και βκ+1 = χ2κ
            akro2 = x2k;
        else
            % ακ+1 = χ1κ και βκ+1 = βκ (δεν αλλάζει)
            akro1 = x1k;
        end
    end

    % Εδώ, έχουμε βρει το τελικό διάστημα οπού υπάρχει το x* = θέση ελαχίστου
    % Πάντα μιλάμε για συγκεκριμένα l,ε που επιλέγω εγώ απαρχής.
    counter = counter + 1;
    akra1(counter) = akro1;
    akra2(counter) = akro2;
    akra1
    akra2
    display(' '); display(' '); display(' '); display(' ');
    result = [];
    result(1) = akra1(length(akra1));
    result(2) = akra2(length(akra2));
    display(strcat('x* between ', num2str(result(1)), ' and ', num2str(result(2))));
    plot(akra1, style)
    xlabel('Steps k')
    ylabel('The two edges: ak and bk')
    title('f1 = (x-2)^2-sin(x+3),  l = {5,10,20}/1000 and epsilon = 1/1000')
    hold on
    plot(akra2, style)
    hold on
    display('*******************************************************')
end