Co=0.1
kc=0.1
V_NaOH=0
x=V_NaOH/20

L0=L1=(1/1000*kc)*(349.8*Co+76.3*Co)

print(paste("VNaOH=0 L=",L0,"L0dil=",L0*20/(V_NaOH+20)))
## [1] "VNaOH=0 L= 0.004261 L0dil= 0.004261"
V_NaOH=10
x=V_NaOH/20
L10=(1/1000*kc)*(349.8*Co*(1-x)+76.3*Co+50.1*x*Co)
print(paste("VNaOH=10 L=",L10,"L10dil=",L10*20/(V_NaOH+20)))
## [1] "VNaOH=10 L= 0.0027625 L10dil= 0.00184166666666667"
V_NaOH=20
x=V_NaOH/20
L20=(1/1000*kc)*(50.1*x*Co+76.3*Co+199.0*(x-1)*Co)
print(paste("VNaOH=20 L=",L20,"L10dil=",L20*20/(V_NaOH+20)))
## [1] "VNaOH=20 L= 0.001264 L10dil= 0.000632"
V_NaOH=30
x=V_NaOH/20
L30=(1/1000*kc)*(50.1*x*Co+76.3*Co+199.0*(x-1)*Co)
print(paste("VNaOH=30 L=",L30,"L30dil=",L30*20/(V_NaOH+20)))
## [1] "VNaOH=30 L= 0.0025095 L30dil= 0.0010038"