По точкам.

Пример 1. Уравнение [MATH]f(x;y;z)←x^2+y^2+z^2+sin(4*x)+cos(3*y)+sin(2*z)-4.486[/MATH]

Документ: Blobby2.sm
Как это выглядит в Mathcad'е: Blobby2.pdf

Пояснения: Берётся одна начальная точка на поверхности. Далее алгоритм скользит вдоль поверхности и вычисляет точки вдоль перпендикулярных направлений. Сначала вычисляется одна линия поверхности как результат пересечения поверхности и плоскости, а затем эти найденные точки кривой используются для построения перпендикулярных линий пересечения поверхности и другой плоскости. Т.о. можно получить облако точек, принадлежащих поверхности.

Пример 2. Уравнение Пфаффа "Камешек"


Документ: Камешек.sm

Проблемы с reverse(). Она не работает в связке со stack(). Такая конструкция не работает в программном блоке:

Пример 3. Подушка

Уравнение: [MATH]f(x;y;z)←y^2+(x^2-y^2)^2+z^2-1[/MATH]

Документ: Подушка.sm

Пример 4. Бублик
Уравнение: [MATH]a←1[/MATH] [MATH]c←1.1[/MATH] [MATH]f(x;y;z)←(c-sqrt(x^4+y^4))^2+sin(z)^2-a^2[/MATH]


Будет желание, посчитайте: Бублик.sm

Бублик особенно крут. Классно!

Бублик построен в Mathcad, также как и мой аватар. Если бы вычисления происходили на низком уровне в неуправляемом коде на уровне ассемблера, то можно было бы строить картинки красивее, а так время тратится на интерпретацию и многочисленных посредников.

[YOUTUBE]http://www.youtube.com/watch?v=EGbHbFCEtR8[/YOUTUBE]

Wrote

Пример 4. Бублик
Уравнение: [MATH]a←1[/MATH] [MATH]c←1.1[/MATH] [MATH]f(x;y;z)←(c-sqrt(x^4+y^4))^2+sin(z)^2-a^2[/MATH]
Будет желание, посчитайте: Бублик.sm

Вячеслав Николаевич, чего это Вы? Картинка совсем не такая. Ведь поверхность неограниченная…

Ваша работа вот найду, покажу... по-моему у Вас там такой же получался, я и не стал дальше копать, тем более для Маткада это делать труднее.

Wrote

Ваша работа вот найду, покажу... по-моему у Вас там такой же получался, я и не стал дальше копать, тем более для Маткада это делать труднее.

Так что же, и я когда-то молод был, горяч, спешил… Значит и там, где было, надо бы исправить. Поэтому и говорю о развитии, а не о повторе старого…

Вот документ с выкладками в Mathcad'е: Бублик.pdf

И у меня и у Вас была нарисована правда, просто от стыка видимо отскочил решатель:
[MATH]f(x;y;z)←(1.1-sqrt(x^4+y^4))^2+sin(z^2)-1^2[/MATH] - я поспешил, тут [MATH]z^2[/MATH], а не [MATH]sin(z)^2[/MATH]


Ааа, я начальную точку взял не точно и ориентация секущих плоскостей возможно неудачна - это yOz и xOy. Может надо было покосее взять. Для сфер формулы получаются страшнее, а самое главное, Маткад их считает дольше. Вот и причина таких "казусов".

Картинка у Алексея правильная, я ошибся при переписывании. Всё равно, построенная мной картинка верна правда не даёт полной картинки и вводит в некоторое заблуждение.


Мою реализацию следует использовать немножко по-другому. При использовании вспомогательных сфер меняется только дифурчик.

Если убрать соединительные линии между отрезками, то будет красиво работать:

Документ: implicitplot2d.sm

Тот же самый алгоритм (Лоренсена и Клайна - Марширующие кубы) для пространства 3D:

Документ: implicitplot3d_параболоид.sm
Дополнительные файлы: edgeTable и triTable

П.С.: Если убрать соединительные линии, то будет красиво.
Алгоритм: http://www.graphicon.ru/oldgr/courses/cg03b/assigns/hw6/

Wrote

...

Не хотелось Вас огорчать, но на всякий случай проверьте картинку, Вас здесь представляющую. Если, конечно, это ей соответствует:
1 + (((x-0.61*cos(0.5*Pi*x))^2 +(y-0.61*cos(0.5*Pi*y))^2-0.9)^2 + (2*z)^2)*(((x^8+y^8+z^8)*0.4^8)^6-1) = 0

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

Вот ещё подобная формула (у меня есть несколько вариаций на эту тему):
[YOUTUBE]http://www.youtube.com/watch?v=ejZv9vcBlgI[/YOUTUBE]

Wrote

Пример 4. Бублик
Уравнение: [MATH]a←1[/MATH] [MATH]c←1.1[/MATH] [MATH]f(x;y;z)←(c-sqrt(x^4+y^4))^2+sin(z)^2-a^2[/MATH]


Будет желание, посчитайте: Бублик.sm

Не смог посчитать. Ошибка - "Значение вышло за пределы максимально допустимого положительного числа" (скриншот прилагаю)

---

Нужно подобрать численный решатель для системы ДУ. В примере он выполнен в виде программы и этот вариант устарел уже, хотя должен быть работоспособен. Это решатель с фиксированным шагом. Можно попробовать другие решатели, в т.ч. для жестких систем. Нужно перебрать решатели из плагинов ODESolvers и DotNumerics, они имеют одинаковый интерфейс вызова. Также нужно уменьшить количество точек и шаг сетки. В любом случае, даже если какая-то картика будет, она не очень прояснит структуру поверхности.

Этот пример довольно сложный для SMath Studio.

Уни,Ваш плагин rkfixed с фиксированным шагом прекрасно работает
с Вашими же примерами.Поверхности отображены на плоские графики.

---

Kameshek1.smz (108 KiB) downloaded 108 time(s).