%
clear
close all
%
x = -3:3;
y1 = bessel(1,x+3);
y2 = besselj(2,x+3);
y3 = besselj(3,x);
y4 = bessel(2,x+2);
%
figure(1)
subplot(2,2,1)
%
plot(x,y1,'--+','LineWidth',2,'MarkerSize',12)
%
axis([-4 4 -.5 1])
set(gca,'XTick',[-3 -2 -1 0 1 2 3])
grid on
%
xlabel('\it t')
title('Time series of {\it x_1}')
%
subplot(2,2,2)
%
plot(x,y2,'--+','LineWidth',2,'MarkerSize',12)
%
axis([-4 4 -.5 1])
set(gca,'XTick',[-3 -2 -1 0 1 2 3])
grid on
%
xlabel('\it t')
title('Time series of {\it x_2}')
%
subplot(2,2,3)
%
plot(x,y3,'--+','LineWidth',2,'MarkerSize',12)
%
axis([-4 4 -.5 1])
set(gca,'XTick',[-3 -2 -1 0 1 2 3])
grid on
%
xlabel('\it t')
title('Time series of {\it x_3}')
%
subplot(2,2,4)
%
plot(x,y4,'--+','LineWidth',2,'MarkerSize',12)
%
axis([-4 4 -.5 1])
set(gca,'XTick',[-3 -2 -1 0 1 2 3])
grid on
%
xlabel('\it t')
title('Time series of {\it x_4}')
%
%print -deps timeseries.eps
