https://web.snauka.ru/issues/2013/02/20410 Mon, 27 Jan 2025 11:00:56 +0000 hourly 1 http://wordpress.org/?v=3.2.1 https://web.snauka.ru/issues/2013/02/20410#comment-1972 amfirifma Sun, 07 Oct 2018 17:45:59 +0000 https://web.snauka.ru/?p=20410#comment-1972 Огромное спасибо авторам! - Очень интересный материал. Пришлось немного помучаться, чтобы разобраться ( правда, не во всем :). Текст Maple-процедур нужно вводить без ">". При построении графика с синусами Maple 17, почему-то, обрезает первый и последний кружок с числом. Сначала лучше копировать текст процедуры в текстовый редактор (например "Блокнот") и сравнить с тем, что написано на сайте. Почему-то происходит замена или потеря символов. (Например: при попытке построить окружность по 10 чисел, сравнимых по модулю с числом 6) Работающий текст процедуры: cr := proc (t, b) local i, pn, p, c1, w, n, j, ch; n := 0; for i to t do for j to b do n := n+1; pn[n] := [i, evalf(2*j*Pi/b), n-1]; ch[n] := pn[n][3] end do end do; w := {seq(pn[k], k = 1 .. t*b)}; for i to nops(w) do if isprime(ch[i]) then p[i] := plots[textplot]([w[i]], coords = polar, color = red) else p[i] := plots[textplot]([w[i]], coords = polar, color = black) end if end do; c1 := seq(p[k], k = 1 .. t*b); plots[display](c1, axes = none) end proc; cr(10, 6); Огромное спасибо авторам! – Очень интересный материал. Пришлось немного помучаться, чтобы разобраться ( правда, не во всем :).
Текст Maple-процедур нужно вводить без “>”.
При построении графика с синусами Maple 17, почему-то, обрезает первый и последний кружок с числом.
Сначала лучше копировать текст процедуры в текстовый редактор (например “Блокнот”) и сравнить с тем, что написано на сайте. Почему-то происходит замена или потеря символов. (Например: при попытке построить окружность по 10 чисел, сравнимых по модулю с числом 6)
Работающий текст процедуры:
cr := proc (t, b) local i, pn, p, c1, w, n, j, ch; n := 0; for i to t do for j to b do n := n+1; pn[n] := [i, evalf(2*j*Pi/b), n-1]; ch[n] := pn[n][3] end do end do; w := {seq(pn[k], k = 1 .. t*b)}; for i to nops(w) do if isprime(ch[i]) then p[i] := plots[textplot]([w[i]], coords = polar, color = red) else p[i] := plots[textplot]([w[i]], coords = polar, color = black) end if end do; c1 := seq(p[k], k = 1 .. t*b); plots[display](c1, axes = none) end proc;
cr(10, 6);

]]>