Freqeuncy modulation & demodulation 1MHz(certain carrier signal) ECG singal

From [http://eleceng.dit.ie/dorran/matlab/ecg.txt] I got ecg signal.

With this I want make a modulation & demodulation MATLAB code and find graph both real time and freqeuncy-domain.

I can not match time and other function index such as Modulation Index other things...

what I have to next?

clc clear all ; close all ; data_2=load('ecg_data.txt'); ecg_data=(data_2(:,1));

% fm=input('Message Frequency=');

fc=input('Carrier Frequency=');

% mi=input('Modulation Index=');

t=0.0001:0.0001:0.6;

m=(data_2)';

subplot(4,1,1);

plot(t,m);

xlabel('millisecond');

ylabel('Amplitude');

title('ECG Signal');

grid on;

t=0:0.0001:0.6; c=sin(2*pi*fc*t)'; subplot(4,1,2); plot(t,c); xlabel('millisecond'); ylabel('Amplitude'); title('Carrier Signal'); grid on;

t=0:0.0001:0.6;

y=c.*(sin(2*pi*t'));%Frequency changing w.r.t Message

subplot(4,1,3);

plot(t,y);

xlabel('millisecond');

ylabel('Amplitude');

title('Signal_Frequency_Modulation');

grid on;

t=0:0.0001:0.6; y=c.*(sin(2*pi*t')).*(sin(2*pi*t'));%Frequency changing w.r.t Message subplot(4,1,4); plot(t,y); xlabel('millisecond'); ylabel('Amplitude'); title('Signal_Frequency_Demodulation');

grid on;