How to remove particular harmonic in the time series?

Given a anomaly time series, for which time series need to be computed by removing annual cycle or 1st harmonic or any particular harmonic. Is the method of computing FFT and making ith & (n+1-i)th coefficients making zero and then iFFT, will give the required result or anything need to be taken care?

clc;clear;clf;
t=0:0.2:4*pi+.1;  l = 64 ;
x=3*sin(t)+2*sin(2*t)+sin(6*t);
X = fft(x);
X2=X; X2(3)=0; X2(64-1)=0;
x2 = ifft(X2);
plot(x,'r'); hold on; plot(x2, 'g');
x3=0*sin(t)+2*sin(2*t)+1*sin(6*t);
plot(x3, 'b');

Answer 1

[KSSV]

The following links might be useful, check for a similar discussion:

https://in.mathworks.com/matlabcentral/answers/355037-i-have-a-signal-with-many-frequencies-how-to-remove-a-particular-frequency-and-reconstruct-the-sign

https://in.mathworks.com/matlabcentral/answers/120715-how-to-remove-particular-frequency-component-from-a-time-series

https://in.mathworks.com/help/signal/ug/filtering-data-with-signal-processing-toolbox.html