#LAPLACE
#USEFLOAT=FALSE
#SIMPLIFY=TRUE
#ERRBND=15
#EXPDSP=5
#PREC=8
#F1= + abs(diff(c_2,v_2)*dv)
#F2=g(x)=
#F3=h(x)=f(x)/g(x)
#F4=matrix(
#F5=vector(
#F6=diff(f)
#F7=diff(g)
#F8=diff(h)
#F9=taylor(f,0)
#F10=taylor(f,1)
#F11=%(run type s:startup-sequence)
#F12=
#F13=
#F14=
#F15=
#F16=
#F17=
#F18=
#F19=
#F20=
c_0=1e-3
c_1(v_0,v_1)=c_0*(v_1)/(v_0+v_1)
c_2(v_0,v_1,v_2)=c_0*(v_1+v_2)/(v_0+v_1+v_2)
c_3(v_0,v_1,v_2,v_3)=c_0*(v_1+v_2+v_3)/(v_0+v_1+v_2+v_3)
c_4(v_0,v_1,v_2,v_3,v_4)=c_0*(v_1+v_2+v_3+v_4)/(v_0+v_1+v_2+v_3+v_4)
dv=1e-4
v_0=21e-3,v_1=1.4e-3,v_2=1.6e-3,v_3=4e-3,v_4=14e-3
dc_1=abs(diff(c_1,v_0)*dv) + abs(diff(c_1,v_1)*dv)
eval(dc_1)
dc_2=abs(diff(c_2,v_0)*dv) + abs(diff(c_2,v_1)*dv)  + abs(diff(c_2,v_2)*dv)
eval(dc_2)
dc_3=abs(diff(c_3,v_0)*dv) + abs(diff(c_3,v_1)*dv)  + abs(diff(c_3,v_2)*dv) + abs(diff(c_3,v_3)*dv)
eval(dc_3)
 + abs(diff(c_4,v_1)*dv)  + abs(diff(c_4,v_2)*dv) + abs(diff(c_4,v_3)*dv) + abs(diff(c_4,v_4)*dv)
a=abs(diff(c_4,v_0)*dv),a
b=abs(diff(c_4,v_1)*dv),b
c=abs(diff(c_4,v_2)*dv),c
d=abs(diff(c_4,v_3)*dv),d
e=abs(diff(c_4,v_4)*dv),e
a+b+c+d+e