#LAPLACE
#USEFLOAT=FALSE
#SIMPLIFY=TRUE
#ERRBND=20
#EXPDSP=7
#PREC=2
#F1=),R_a(
#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=
#F12=
#F13=
#F14=
#F15=
#F16=
#F17=
#F18=
#F19=
#F20=
R_VM=15000,R_AM=4000
dI=0.015*50e-6,dU=0.15
R_a(U,I)=U/(I-U/R_VM)
R_b(U,I)=(U-R_AM*I)/I
dR_a(U,I)=abs(diff(R_a,I))*dI+abs(diff(R_a,U))*dU
dR_b(U,I)=abs(diff(R_b,I))*dI+abs(diff(R_b,U))*dU
R_a(6,16.5e-3)
dR_a(6,16.5e-3)
R_b(8,1.2e-6)
dR_b(8,1.2e-6)
R=set(R_a(10,0.095),R_a(9,0.090),R_a(8,0.084),R_a(7,0.078),R_a(6,0.072),R_a(5,0.064),R_a(4,0.057),R_a(3,0.048),R_a(2,0.040),R_a(1,0.028),R_a(0.5,0.021),R_a(0.2,0.014))
dR=set(dR_a(10,0.095),dR_a(9,0.090),dR_a(8,0.084),dR_a(7,0.078),dR_a(6,0.072),dR_a(5,0.064),dR_a(4,0.057),dR_a(3,0.048),dR_a(2,0.040),dR_a(1,0.028),dR_a(0.5,0.021),dR_a(0.2,0.014))

dR_1=set(dR_a(11.5,0.091),dR_a(11,0.087),dR_a(10,0.080),dR_a(9,0.072),dR_a(8,0.063),dR_a(7,0.056),dR_a(6,0.048),dR_a(5,0.040),dR_a(4,0.032),dR_a(3,0.024),dR_a(2,0.016),dR_a(1,0.009))
R_1=set(R_a(11.5,0.091),R_a(11,0.087),R_a(10,0.080),R_a(9,0.072),R_a(8,0.063),R_a(7,0.056),R_a(6,0.048),R_a(5,0.040),R_a(4,0.032),R_a(3,0.024),R_a(2,0.016),R_a(1,0.009))
average(R_1)
eval(error(R_1))

dR_2=set(dR_b(11.5,0.087),dR_b(11,0.084),dR_b(10,0.077),dR_b(9,0.069),dR_b(8,0.061),dR_b(7,0.054),dR_b(6,0.046),dR_b(5,0.038),dR_b(4,0.031),dR_b(3,0.024),dR_b(2,0.016),dR_b(1,0.008))
R_2=set(R_b(11.5,0.087),R_b(11,0.084),R_b(10,0.077),R_b(9,0.069),R_b(8,0.061),R_b(7,0.054),R_b(6,0.046),R_b(5,0.038),R_b(4,0.031),R_b(3,0.024),R_b(2,0.016),R_b(1,0.008))
average(R_2)
eval(error(R_2))

dR_2a=set(dR_a(11.5,0.82e-3),dR_a(11,0.79e-3),dR_a(10,0.72e-3),dR_a(9,0.64e-3),dR_a(8,0.58e-3),dR_a(7,0.50e-3),dR_a(6,0.44e-3),dR_a(5,0.36e-3),dR_a(4,0.30e-3),dR_a(3,0.22e-3),dR_a(2,0.15e-3),dR_a(1,0.08e-3))
R_2a=set(R_a(11.5,0.82e-3),R_a(11,0.79e-3),R_a(10,0.72e-3),R_a(9,0.64e-3),R_a(8,0.58e-3),R_a(7,0.50e-3),R_a(6,0.44e-3),R_a(5,0.36e-3),R_a(4,0.30e-3),R_a(3,0.22e-3),R_a(2,0.15e-3),R_a(1,0.08e-3))
average(R_2a)
eval(error(R_2a))
count(dR_2a)-count(R_2a)
dR_2b=set(dR_b(11,48.5e-6),dR_b(10,44.0e-6),dR_b(9,40.0e-6),dR_b(8,35.5e-6),dR_b(7,31.0e-6),dR_b(6,26.5e-6),dR_b(5,22.5e-6),dR_b(4,18.0e-6),dR_b(3,13.5e-6),dR_b(2,9.0e-6),dR_b(1,5.0e-6))
R_2b=set(R_b(11,48.5e-6),R_b(10,44.0e-6),R_b(9,40.0e-6),R_b(8,35.5e-6),R_b(7,31.0e-6),R_b(6,26.5e-6),R_b(5,22.5e-6),R_b(4,18.0e-6),R_b(3,13.5e-6),R_b(2,9.0e-6),R_b(1,5.0e-6))
average(R_2b)
eval(error(R_2b))