function q = fhansen(c,param)
%
% Input function to CSOLVE in SECTION4.M
%

W_ss 	= param(1);
c_ss 	= param(2);
k_ss 	= param(3);
AN	= param(4);
beta  	= param(5);
sigma	= param(6);
alpha	= param(7);
k_t_1	= param(8);
theta_t	= param(9);
mu1  	= param(10);
mu2  	= param(11);
mu3  	= param(12);
mu4  	= param(13);
mu5  	= param(14);
delta	= param(15);


q = mu3 * ( (1/beta) * c^(-sigma) - W_ss ) + mu4 * (   -c +...
    theta_t * k_t_1^(alpha) * (  AN * c^(sigma) /...
    (theta_t * k_t_1^(alpha) * (1-alpha))  )^((alpha-1)/alpha) + ...
    (1 - delta) * k_t_1 - k_ss   ) + mu5 * log(theta_t);

