h : iç içe iki küre arasındaki mesafe (içteki küre yarıçapı 1)
a : (0,1) noktası ile dıştaki küreyi kestiği nokta arasındaki vektörün açısı
l : slant path uzunluğu

1. x^2 + y^2 = (1+h)^2
2. y = tan(a)*x + 1

+
——————————-

t : tan(a)

x = (-t + sqrt(t^2 + (1+t^2) h (h+2)) ) / (t^2 + 1)

y = t x + 1

l = sqrt((y – h)^2 + x^2)

 

 

function l = slantpath(a)
    t = tan(a);
    x = (-t + sqrt(t*t + (1+t*t)*h*(h+2))) / (t*t + 1);
    y = t*x + 1
    l = sqrt((y*h)^2 + x^2);
end