function [HCO3]=HCO3(x)
    
HCO3=((1-gaussmf(x,[13 25])).*(x>25)+(1-gaussmf(x,[10 25])).*(x<25));