Lab MN

function [a, b] = prob1(x, y, m)
        n = length(x);
    y = y';
        a = zeros(m, 1);
        b = zeros(m, 1);
        for j=1:m
            a(j) = (1/m)*y*cos(j*x);
            b(j) = (1/m)*y*sin(j*x);
        endfor
        a0 = sum(y)/n;
        if n == 2*m
            a(m) = a(m)/2;
        endif
endfunction