clear all
close all

clc

syms y pie ir real;    % definition of variables
syms kappa sigma mu  nu real; % definition of parameters
syms e u real;  % definition of shocks

%equations: expected variables = 0
PC = -pie + kappa * y + e;              % phillips curve
IS = - y + (-1/sigma) * ir + u;         % IS curve
MR = - ir + mu * pie + nu * y;  % monetary rule

%sigma = 0;
%mu = kappa;

sol = solve(PC,IS,MR,'y,pie,ir');   % symbolic solution

% solution for interest rate, inflation and output
disp('Reseni pro ir:')
pretty(sol.ir)           

disp('Reseni pro pie:')
pretty(sol.pie)

disp('Reseni pro y:')
pretty(sol.y)

shock2y =  [diff(sol.y,e), diff(sol.y,u)];
shock2pie = [diff(sol.pie,e), diff(sol.pie,u)];
shock2ir = [diff(sol.ir,e), diff(sol.pie,u)];