function [cx,vat,vatx] = deval_steady_state_helper(alpha,lambda,deltax,delta,sigma,eta,alphax,lambdax,n,epsilon) % Defining the function to be used in the mod-file
options=optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
guess = [0.8, 0.2, 0.2];
[cx,vat,vatx] = fsolve('m1', guess, options);
end