function [PaO2]=PaO2(x)
PaO2=((1-gaussmf(x,[80 85])).*(x>85)+(1-gaussmf(x,[22 85])).*(x<85));