python3 <<'PY'
import math
a=0.25; b=1.0/12
y1c=0.024697868008998203; y2c=0.23045484489188828
def Aub(y1,y2):
    al=2*a+2*b-y1-y2
    return math.sqrt(al)*(a+b+math.sqrt(max(0,b*(b-y1)))+math.sqrt(max(0,a*(a-y2))))
def BC(y1,y2):
    lo=0.0; hi=a; gr=(math.sqrt(5)-1)/2
    def f(y3):
      beta=a+b-y3; B=math.sqrt(beta)*(b+math.sqrt(max(0,a*(a-y3))))
      gam=y1+y2+y3; C=math.sqrt(gam)*(math.sqrt(max(0,b*y1))+math.sqrt(max(0,a*y2))+math.sqrt(max(0,a*y3)))
      return B+C
    for _ in range(80):
      c=hi-gr*(hi-lo); d=lo+gr*(hi-lo)
      if f(c)<f(d): lo=c
      else: hi=d
    return f((lo+hi)/2)
for d in [-0.002,-0.001,-0.0005,-0.0002,-0.0001,-0.00005,0,0.00005,0.0001,0.0002,0.0005,0.001,0.002]:
  y1=y1c+d; y2=y2c-3*d
  if 0<=y1<=b and 0<=y2<=a:
    ac=Aub(y1,y2); bc=BC(y1,y2)
    print("d",d,"Ac",ac,"BC",bc,"ubtotal",ac+bc)
PY