L10 - CN

clc
clear

rel = 0;
while(~rel)
    fprintf('Nr. noduri de interpolare(ni):');
    ni = input(' ');
    rel = 0;
    while(~rel)
        xi = zeros(ni, 1);
        yi = zeros(ni, 1);
        for i = 1: ni
            fprintf('Abscisa x%g = ', i);
            xi(i, 1) = input(' ');
            fprintf('Ordonata y%g = ', i);
            yi(i, 1) = input(' ');
        end
        fprintf('Introduceti 1 pentru continuare sau 0 pentru reintroducere perechile de valori pentru interpolare: ');
        rel = input(' ');
    end
    rel = 1;
    while(~rel)
        fprintf('Pas evaluare (heval): ');
        heval = input(' ');
        xe = xi(1) : heval : xi(ni);
        npev = length(xe);
        if(xe(npev) < xi(ni))
            npev = npev + 1;
            xe = [xe xi(ni, 1)]
        end
        npev
        rel = 0;
        while(~rel)
            fprintf('Grad polinom de interpolare(mg, 1 <= mg <= ni - 1):');
            mg = input(' ');
            ye = zeros(npev, 1);
            for ip = 1: npev
                ye(ip) = INTPL(ni, xi, yi, mg, xe(ip));
                xe(ip)
                ye(ip)
            end
            rel = 0;
            while(~rel)
                plot(xi, yi, 'or', xe, ye, '-b');
                fprintf('Introduceti 1 pentru continuare sau 0 pentru reluare reprezentare grafica:');
                rel = input(' ');
            end
            fprintf('Introduceti 1 pentru continuare sau 0 pentru a relua introducerea altui grad pentru polinomul de interpolare:');
            rel = input(' ');
        end
        fprintf('Introduceti 1 pentru continuare sau 0 pentru reluare cu alt pas de evaluare');
        rel = input(' ');
    end
    fprintf('Introduceti 1 pentru continuare sau 0 pentru reluarea introducerii altei functii:');
    rel = input(' ');
end