Problema1 - Lab10

function [c, a, b] = interpolare_trigonometrica(m, y)
    for j=1:2*m
        x(j) = -pi + ((j-1)/m)*pi;
    end
    n = length(x);
    x = x';
    w = 2*pi/n;
    a = zeros(m, 1);
    b = 1;
    for j=1:m
        a(j) = y*cos(j*x);
        b(j) = y*sin(j*x);
    end
    a = 2*a/n;
    b = 2*b/n;
    a0 = sum(y)/n;
    if n == 2*m
        a(m) = a(m)/2;
    end
    c(1) = a0 + a(1)*cos(x(1));
    c(2) = b(j)*cos(x(1))
    for j=3:m
        if(mod(j, 2) == 1)
            c(j) = c((j+1)/2-1) + a((j+1)/2)*cos((j+1)/2*x);
        else
            c(j) = c(j/2-1) + b(j/2)*sin(j/2*x);
        end
    end
end