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