%this file permits to calculate the value for 
%omega1 and mu1 in Mary G finn "Perfect Competition
%and the effects of energy price increases on economic activity"
%(2000) Journal of Money Credit and Banking vol.32
%
%
% equation 24 may be written as:
% omega1*delta+p*mu1*a=(1-theta)*y/k
% 
%we can rewrite equation 21 as:
%
%delta=h^(omega1)/omega1
%
%by the same token, we can rewrite equation 22 as:
%
%a=h^(mu1)/mu1
%
% h is a scalar such that h=s*u
%
%given the different parameter page 407 we have a system of 3 equations:
%0.025*omega1+1*mu1*0.0059-(1-0.7)/7.3217=0
%h^(omega1)/omega1-0.025=0
%h^(mu1)/mu1-0.0059=0
%
%we define x(1)=omega1 ; x(2)=mu1 ; x(3)=h. We can solve this system:


function F = myfun(x)
F=[0.025*x(1)+0.0059*x(2)-(.3/7.3217);
    (x(3)^x(1))/x(1)-0.025;
    (x(3)^x(2))/x(2)-0.0059];
end
    
