Exponenta.ru

В коде исправьте sol(x)[2] на rhs(sol(x)[2])

Нет.

1. a это имя списка [[0,0],0] , который мы вводим для дальнейшего поиска позиции (номера строки и столбца) максимального по модулю элемента матрицы А . 2. r это выражение, стоящее в знаменателе, когда ищем нужный угол для поворота (см. формулу в вики). Например, если оказалось, что a=[[2,3],7] , то ...

Проверьте исправленный вариант: Jacobi:=proc(A::Matrix) local n, a, i, j, p, c, s, G, r; uses LinearAlgebra; n:=RowDimension(A); a:=[[0,0], 0]; for i to n do for j from i+1 to n do if abs(A[i,j])>=abs(a[2]) then a:=[[i,j], A[i,j]] fi; od; od; r:=A[a[1,1],a[1,1]]-A[a[1,2],a[1,2]]; p:=`if`(r<>0,2*a[2]...

Вставьте в процедуру строку print(G) после строки G:=...

Вот простая процедура, которая делает одну итерацию в методе вращения Якоби для поиска собственных значений симметрической вещественной матрицы. Сначала в двойном цикле находится максимальный внедиагональный элемент матрицы. Затем использованы формулы из вики для нахождения матрицы G . Применяя эту ...

Ниже примеры как это делается:

plot(x^2, x=0..1, axis=[thickness=0]); # Самые тонкие оси
plot(x^2, x=0..1, axis=[thickness=2]);
plot(x^2, x=0..1, axis=[thickness=4]);

См. справку по запросу plot,axis

Ваш код не работает, т.к. вы забыли разделитель (точка с запятой) когда вводили вектор b . Также я сделал несколько исправлений в коде, связанных в основном с тем, что пакет linalg давно устарел (впрочем в коде он и не нужен). Вместо matrix и vector следует использовать Matrix и Vector . Что означае...

Если я правильно Вас понял, то Вы хотите, чтобы ось y шла не вверх, а вниз. В 2D это сделать невозможно, т.к. направления осей координат всегда стандартное. Но это можно сделать, если выйти в 3D и использовать опцию orientation . Вот пример: plots:-spacecurve([x,x^2,0], x=0..2, color=red, axes=norma...

В Мэйпл это также легко решается:
solve(max(abs(x), abs(x-8))=5);
3, 5

Ну это в 100 раз проще: Primes:=[``, seq(ithprime(k), k=1..100)]: plots:-animate(plots:-textplot,[[0,0, Primes[round(i)]], color=blue, font=[times,roman,30]], i=1..101, frames=500, paraminfo=false, axes=none, size=[100,100]); https://a.radikal.ru/a32/1805/f7/af63e75e9af2.gif

Будет более наглядно (и проще для реализации), чтобы элементы списка появлялись не в одной вертикальной колонке, а построчно, т.е. как обычный текст слева направо. Вот пример такой анимации - список первых 100 простых чисел: restart; L:=[[-50,0,convert([seq(ithprime(k),k=1..100)], string)],color=blu...

Не пробовал. Почему-то не подумал об этом. Команда CurveFitting:-Spline позволяет автоматически получить сплайн нужного порядка (в данном примере - сплайн первого порядка, т.е кусочно линейную функцию).

Вот другое решение, дающее более плавную анимацию. Для этого с помощью процедуры ConvertToPiecewise , каждый список, определяющий ломаные линии, конвертируется в piecewise функции: restart; ConvertToPiecewise:=proc(L::list,v::name:='t') local n; n:=nops(L); piecewise(seq(op([v>=L[i,1] and v<L[i+1,1]...

Вот код для параллельной анимации трёх ломаных: restart; t := [1,2,3,4,5,6,7,8,9,10]: x:= [2,1,-3,-8,-5,-4.-8,-4,6,2,9]: y:= [2,5,3,8,5,4.3,4,6,3,4]: z:= [3,5,3,1,5,4.3,5,4,3,1]: tx:=zip(`[]`, t, x): ty:=zip(`[]`, t, y): tz:=zip(`[]`, t, z): A:=plots:-animate(plot,[tx[1..floor(i)], color=red, thickn...