function [PaCO2]=PaCO2(x)
    
PaCO2=((1-gaussmf(x,[10 40])).*(x>40)+(1-gaussmf(x,[10 40])).*(x<40));