function [HR]=HR(x)
HR=((1-gaussmf(x,[35 90])).*(x>90)+(1-gaussmf(x,[25 90])).*(x<90));