itload('merg_egiw.it');

figure(1);
subplot(1,3,1);
hold off
plot(G(:,1),'-')
hold on
plot(P1(:,1),'+-')
plot(M(:,1),'--')
plot([1,100],[1.5,1.5],':')
xlabel('time [t]')
title('Expected value of a_1')

subplot(1,3,2);
hold off
plot(G(:,2),'-')
hold on
plot(P2(:,1),'+-')
plot(M(:,3),'--')
plot([1,100],[0.8,0.8],':')
xlabel('time [t]')
title('Expected value of a_2')

subplot(1,3,3);
hold off
plot(G(:,3))
hold on
plot(P1(:,2),'--')
plot(P2(:,2),'+-')
plot(M(:,2),'--')
plot([1,100],[0.01,0.01],':')
xlabel('time [t]')
title('Expected value of \sigma')