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')