clear,clc
fx1=@(b,x)b(1)+b(2)*x+sqrt(b(3)+b(4)*x+b(5)*x.^2);
fx2=@(b,x)b(1)+b(2)*x-sqrt(b(3)+b(4)*x+b(5)*x.^2);
figure(1),clf
b=[304.896155 -0.0468025122 -25124.8078 330.246104 -1.05243991];
X=xlsread('tynh','a1:b156');
x=X(:,1);y=X(:,2);
rg=range(x);
x1=linspace(min(x)-rg/100,max(x)+rg/100,200);
plot(x,y,'o','markersize',10,'markerfacecolor','k')
hold on
y1=fx1(b,x1);y2=fx2(b,x1);
plot(x1,y1,'r-',x1,y2,'r-','linewidth',2.5)
axis tight
figure(2),clf
b=[372.4589157 -0.05152819875 -25265.89178 321.9513439 -0.9973063478];
X=xlsread('tynh','d1:e269');
x=X(:,1);y=X(:,2);
rg=range(x);
x1=linsp