function [Na]=Na(x)
Na=(1-gaussmf(x,[18 140])).*(x>140)+(1-gaussmf(x,[15 140])).*(x<140);