clear
clc

% Amplitude
Amplitude=3;

% n
n=1;

% c
c=2; % m/s

% L
L=4;

% Time
Time=0:1/1000:2;

for n=1:3
   
	% Frequency
	BasicFrequency=(n*c)/(2*L);

	% Angular Frequency
	AngularFrequency=2*pi*BasicFrequency;

	Wave=Amplitude.*sin(AngularFrequency.*Time);
   
   figure(1);
	subplot(3,1,n);
   plot(Time,Wave);
   title(['Normal Mode for n=',num2str(n)]);
   ylabel('Amplitude');
   grid on
   
end

xlabel('String Length');
