function []=logPH(); set(gca,'xscale','log'); set(gca,'yscale','log');