function [Ptopt, Ptotalmin] = findptopt()
global log10pe;
global Pt;
global gamma;
global Ptoptfinal2;

       [Ptopt, Ptotalmin] = fminbnd(@findtotalpower,0.01,100);    