I want to get implus from sin

what is modification on this code to get dirac delta from pure sin wave?

 fc=10,fs=2*fc=20
t=1:1/fs:1
 a=sin(2*pi*fc*t)
 b=abs(fft(a).^2);