function [Mm] = oun_p0vonp0m2Mm(gamma,p0vonp0m) %oun_p0vonp0m2Mm % Ma=1; Mb=5; p0vonp0ma=oun_p0vonp0m(gamma,Ma); p0vonp0mb=oun_p0vonp0m(gamma,Mb); Mc=(Ma+Mb)/2; p0vonp0mc=oun_p0vonp0m(gamma,Mc); eps=1e-3; while abs(p0vonp0mc-p0vonp0m)>eps if p0vonp0mc>p0vonp0m Ma=Mc; p0vonp0ma=p0vonp0mc; else Mb=Mc; p0vonp0mb=p0vonp0mc; end Mc=(Ma+Mb)/2; p0vonp0mc=oun_p0vonp0m(gamma,Mc); end Mm=Mc; end