function Sfm=My_FM_moulation1(fc,k,t,Sm);
% t=[0:0.001:1];
% Sm=cos((2*pi*fc*t)-msin(2*pi*fm*t)); %Massage Signal%
subplot(3,1,1);
plot(t,Sm);
xlabel('Time');
ylabel('Amplitude');
title('Massage Signal');
grid on ;

Sc=cos(2*pi*fc*t); %Carrier Signal%
subplot(3,1,2);
plot(t,Sc);
xlabel('Time');
ylabel('Amplitude');
title('Carrier Signal');
grid on ;

Sfm=cos(1+k*Sm).*Sc;  %FM Modoulation%
subplot(3,1,3);
plot(t,Sfm);
xlabel('Time');
ylabel('Amplitude');
title('FM Modoulation');
grid on ;