function F = call_csolve4(F_bar,sp_bar,MU)

%uses the csolve command and the function "frac_default" to solve for sigma_bar
%sigma_bar = csolve(@(sigma_bar) frac_default2(F_bar,sp_bar,MU, sigma_bar),0.5,[],1e-8,200);
sigma_bar = csolve(@(sigma_bar) frac_default2(F_bar,sp_bar,MU, sigma_bar),0.2,[],1e-8,200);
F = sigma_bar;
end