function F = call_csolve3(F_bar,omega_bar)

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