function[y] = cos0approx(x, n)

y=0;

for j = 0:1:n
    
    y = y+((-1)^j)*((x^(2*j))/(factorial(2*j)));
end
