fit_g1-dis_rss.html

> g1:=a*x^b*(1-x)^3*(1+c/q2);

> dya:=diff(g1,a);

> dyb:=diff(g1,b);

> dyc:=diff(g1,c);

> dy:= sqrt((dya*da)^2+(dyb*db)^2+(dyc*dc)^2+2*dya*dyb*dab+2*dya*dyc*dac+2*dyc*dyb*dcb);

dab=0.015190159,dac=-0.16542670,dcb=-0.028085555

dab=0.0,dac=0.,dcb=0.0

> y1:=subs(a=0.39195223,b=0.067564942,c=0.063631241,q2=1.279,g1);
dy1:=subs(a=0.39195223,b=0.067564942,c=0.063631241,
da=sqrt(0.064588362),db=sqrt(0.0070870338),dc=sqrt(0.46407566),
dab=0.015190159,dac=-0.16542670,dcb=-0.028085555,q2=1.279,dy);

> plot({y1,y1+dy1,y1-dy1},x=0..0.3);

> Gamma1_err:=int(dy1,x=0.00..0.316);

> Gamma1_p:=int(y1,x=0..0.316);

> yn:=subs(a=-0.0000734,b=-1.933556,c=2.837267,q2=1.279,g1);
dyn:=subs(a=-0.0000734,b=-1.933556,c=2.837267,
da=sqrt(0.000000434),db=sqrt(2.134173),dc=sqrt(2643.0566),
dab=0.00035987,dac=0.0304052,dcb=55.85988,q2=1.279,dy);

> plot({yn,yn+dyn,yn-dyn},x=0.01..0.3);

> Gamma1_n_err:=int(dyn,x=0.001..0.316);

> Gamma1_n:=int(yn,x=0.001..0.316);

> yd:=subs(a=1.77767124,b=0.73906325,c= -0.78551499,q2=1.279,g1);
dyd:=subs(a=1.77767124,b=0.73906325,c= -0.78551499,
da=sqrt(3.88586297),db=sqrt(0.16549087),dc=sqrt(0.065409055),
dab=0.65374959,dac=-0.35940478,dcb=-0.024749276,q2=1.279,dy);

> plot({yd,yd+dyd,yd-dyd},x=0.001..0.3);

> Gamma1_d_err:=int(dyd,x=0.001..0.316);

> Gamma1_d:=int(yd,x=0.001..0.316);