このときピストンの位置xは
https://i.imgur.com/Sj9kJIp.jpg

であり、伝統的にはこれを三角関数の足し合わせに展開したあと2次(cos2θ)の項までを用い、
https://i.imgur.com/Imj1bPm.jpg

と近似して計算します

両者をMaxima用の数式として表すと

ややこしい方の式：
r*(1-cos(%theta)+%lambda*(1-(1-sin(%theta)^2/%lambda^2)^(1/2)))

近似式：
r*(1-cos(%theta)+1/(4*%lambda)*(1-cos(2*%theta)))

です