Прямоугольник и вписанный треугольник Ч. 2

Аналитика оказалась настолько сложной, что оставил её до лучших времён.
Но любимый мой метод Монте Карло опять творит чудеса! Достаточно составить простенькую программу

s1=28:s2=48:s3=80
s38=s1+s2+s3
a0=10.1:b0=10:x0=2:y0=4
z=0.001
s5=10^10
for i=1 to 100000000
a=a0*(1+z*(ran()-.5))
b=b0*(1+z*(ran()-.5))
x=x0*(1+z*(ran()-.5))
y=y0*(1+z*(ran()-.5))
s31=x*b/2:s32=(a-x)*y/2:s33=(b-y)*a/2
z1=abs(s31-s1)+abs(s32-s2)+abs(s33-s3)
s4=a*b-s38
 l=sqrt(x^2+b^2):k=sqrt((a-x)^2+y^2)
t=sqrt((b-y)^2+a^2)
s3a=1/4*sqrt((l+k+t)*(l+k-t)*(l-k+t)*(-l+k+t))
if z1<s5 then
print a,b,x,y,s3a,s38,z1
s5=z1:a0=a:b0=b:x0=x:y0=y
if z1<0.0005 then z=0.000001:fi
fi
next i

Концовка вычислений и комментарии - на рисунке. Исходные данные взял от фонаря. Решение проверил на бумаге. Оказалось все верно!
Но так хочется получить компактную формулу.

12 июля 2022 г.


Рецензии