Расчёт в wxMaxima

а есть ли решение? :)

Должно быть, через любые две точки можно построить окружность, если радиус позволяет.

что-то сам не знаю, тоже не получается... а попробуйте пример что ли попроще =) с цифрами поменьше. на который точно знаете ответ

ага... в общем да.. (я уже что-то плоховато школьную геометрию помню).

но этот радиус единственный, или их бесконечное множество? то есть окружность такую можно одну провести, или их тоже много?..

в вашем примере радиус вы точно задали 97, я что-то думаю, что вполне возможно такого и не существует...

вот такой «детский» пример:

solve([(1-x)^2+(2-y)^2=5^2,(1+x)^2+(2-y)^2=5^2],[x,y]);

В ответ

[[x=0,y=2-2*sqrt(6)],[x=0,y=2*sqrt(6)+2]]

Похоже на правду.

А вот мой упрощённый пример:

solve([(2-x)^2+(2-y)^2=1, (3-x)^2+(1-y)^2=1],[x, y]);

[[x=2,y=1],[x=3,y=2]]

Работает же!

но этот радиус единственный, или их бесконечное множество? то есть окружность такую можно одну провести, или их тоже много?..

С разными радиусами — бесконечно много, главное, чтобы радиус был больше или равен длине отрезка между точками, а то не получится через обе точке построить окружность. Но радиус в моём случае строго указан, поэтому можно построить только две (так, что дуга между точками получится вогнутая и выпуклая)

в вашем примере радиус вы точно задали 97, я что-то думаю, что вполне возможно такого и не существует...

Графически (набросок в тетрадке) его существование подтверждается.

Спасибо за помощь!

Завтра, похоже, придётся самостоятельно подсчитать. Только вторая лаба по разработке УП для ЧПУ, а я застрял на весь вечер с этими дурацкими кругами %)

В wxmaxima можно тоже графики строить... в том числе через gnuplot. Возможно, построить два графика и найдёте решение? :)

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

то есть вы хотите сказать, что уравнение решение _имеет_, но maxima его не взяла? а маткад взял? интересно...

мб повод для баг-реппорта?

Ага, имеет, в этом легко убедиться, построив окружности через эти точки. Подрисовано коряво, но позволяет понять, что вообще за задача передо мной стояла.

Возможно, похоже на баг.

Упрощенный код ошибки:

f1:(x-1)^2+(y)^2=2^2;
f2:(x)^2+(y-1)^2=2^2;
solve([f1,f2],[x,y]);
[f1_1,f1_2]:map(lambda([f],last(f)),solve(f1,y));
[f2_1,f2_2]:map(lambda([f],last(f)),solve(f2,y));
plot2d([f1_1,f1_2,f2_1,f2_2],[x,-2,3]);
build_info();

Во, а я как раз собирался поднять тот багрепорт, но разработчики, разумеется, уже в курсе.

use_grobner=true

Primarily, this option provides a workaround for weakness in the function algsys.

Спасибо за найденное решение.