function y=robs(theta1, theta2, seed);
y1=robusts(theta1, seed);
y2=robusts(theta2, seed);

theta1
theta2
T=50; % Length of simulation
% Plotting
subplot(3,2,1);
plot(1:T, y1(:,1), 'b-', 1:T, y2(:,1), 'r-');
title('Consumption: Blue for theta1');

subplot(3,2,2);
plot(1:T, y1(:,2), 'b-', 1:T, y2(:,2), 'r-');
title('Labor: Red for theta2');

subplot(3,2,3);
plot(1:T, y1(:,3), 'b-', 1:T, y2(:,3), 'r-');
title('Tax');

subplot(3,2,4);
plot(1:T, y1(:,4), 'b-', 1:T, y2(:,4), 'r-');
title('Gov Expenditure');

subplot(3,2,5);
plot(1:T, y1(:,5), 'b-', 1:T, y2(:,5), 'r-');
title('Interst Rate');

subplot(3,2,6);
plot(1:T, y1(:,6), 'b-', 1:T, y2(:,6), 'r-');
title('Government Debt');