function cn = anbnToCn(a0, an, bn)
    % Converts Fourier coefficients from an, bn to cn
    N = length(an);
    cn = zeros(2*N+1, 1); % Preallocate cn
    cn(N+1) = a0/2; % c0 at the center
    for n = 1:N
        cn(N+1+n) = 0.5 * (an(n) - 1i * bn(n)); % cn for n > 0
        cn(N+1-n) = conj(cn(N+1+n)); % cn for n < 0
    end
end