[slgu]

fx=@(,x)(x<=16).*(b(1)*x+b(2)*x.^2+b(3)*x.^3+b(4)*x.^4) +(x>16).*((b(5)+b(6)*x+b(7)*x.^2+b(8)*x.^3+b(9)*x.^4)./(1+b(10)*x+b(11)*x.^2+b(12)*exp(b(13)*x+b(14)*x.^2)));
b=[0.853235684485670 -0.116843442327784 0.00834771049211092 -0.000272675227862312 3.65963514674834 0.0864733971875504 -0.00333133191436954 3.04447085911928e-05 -8.56724687299949e-08 -0.0143860176421765 5.30672783303945e-05 196.513295363632 -0.293982163705569 0.00149387237578180]
R^2=0.99949