function [pH]=pH(x)
pH=((1-gaussmf(x,[.15 7.41])).*(x>7.41)+(1-gaussmf(x,[.11 7.41])).*(x<7.41));