clear all;
clc;

options=optimset('display','iter');%,'TolFun',1e-4,'TolX', 1e-4);
options=optimset(options,'MaxFunEvals',10000);

f_steady=@(initial)steady_state(initial);
initial=[71.3721,1,1];
[dopt,fval1,exitflag1]=fsolve(f_steady,initial,options);
exitflag1