function [Ptopt, Ptotalmin] = findptopt2()

global Shlt;

[Ptopt, Ptotalmin] = fminbnd(@findtotalpower2,Shlt,100,optimset('TolX',1e-20));    