Вопросы по теме 'maxima'
Создайте %union в списке в максимумах
Функция to_poly_solve в maxima возвращает %union , но я хочу работать со списком.
Как можно превратить %union в списки, если listify не работает?
413 просмотров
schedule
25.04.2023
График решений нелинейных дифференциальных уравнений
У меня есть система дифференциальных уравнений в Maxima. И я пытаюсь найти решения.
diff_eq1: 'diff(p(t),t) = (5/2 + (3^(1/2))/24 - (5/8)*p(t) - ((3^(1/2))/24)*q(t)) * p(t);
diff_eq2: 'diff(q(t),t) = (7/8 + (3*(3^(1/2))/2) - (3*(3^(1/2))/8)*p(t)...
1519 просмотров
schedule
07.11.2022
Sage: Catch MACSYMA-QUIT не определен
Моя ситуация
Когда я использую модуль sage.all только в обычной оболочке python , sage.all выдает ошибку при запуске методов, таких как solve или assume :
>>> import sage.all as sg;a=sg.var('a');sg.assume(a==0);
Traceback...
639 просмотров
schedule
17.03.2023
Найти ВСЕ локальные максимумы функции
Привет всем, я пытаюсь найти ВСЕ максимумы функции и построить точки, где находятся максимумы.
Функция:
U[x_,y_,a_]:=-((1-a)/Sqrt[(x-a)^2+y^2])-a/Sqrt[(x+1-a)^2+y^2]- 0.5*(x^2+y^2);
Конечно, это 3D-функция, но я также доволен, если...
1213 просмотров
schedule
18.10.2023
wxMaxima сравнивает логические выражения
Я новичок в wxMaxima и не могу понять, как сравнивать два логических выражения. Например:
test1:(A and B) or C;
test2:(A or C) and (B or C);
is(test1=test2);
должно быть
true
но результат
false
Что я здесь делаю не так?!
614 просмотров
schedule
13.02.2023
Более компактное решение в Maxima
У меня есть следующий код:
for n:1 thru 11 do for j:1 thru 21 do v[n,j]:1/sqrt(dp)*
(sum(eigenfunctionsort[n,j]*exp(%i*2*%pi*m*x/dp),m,-10,10));
Там, где собственная функция сортировки определена ранее, x — это переменная, по которой я буду...
111 просмотров
schedule
08.06.2023
Матрицы произвольного размера в максимумах
Я хочу сделать некоторые вычисления с матрицами произвольного размера. Простой пример - возьмите две матрицы NxM и MxK с произвольными элементами и просмотрите элемент произведения как сумму.
Но я не могу найти способ выполнять такие символьные...
315 просмотров
schedule
28.04.2022
Диагонализация симметричной матрицы в Maxima
Вот моя проблема: у меня есть реальная симметричная матрица M , зависящая от двух параметров a,b (которые считаются положительными), и я хочу найти ортогональную матрицу P такой, что PMP^{-1} — диагональная матрица. Вот пример того, что я...
1355 просмотров
schedule
07.03.2023
кодирование на языке maxima vs lisp
Я просто подумал написать какую-нибудь функцию, похожую на функции разделения Mathematica с передачей опции в максимумах в виде,
listpartitionpad(l,n,k,d):= block([temp:[],gap,newl,ntemp:[]],...
449 просмотров
schedule
14.07.2023
скрипт vim для вставки формул ascii-art
Я стараюсь использовать Vim как можно чаще, чтобы изучить его. Сегодня мне пришлось написать множество математических функций. Я написал их как f(x)=e^(x/sqrt(x)) (например, бессмысленная функция), но это очень трудно читать.
Моя идея...
606 просмотров
schedule
06.05.2023
Как я должен оценить функцию в точке в Maxima?
Я хочу взять производную функции омега по направлению в точке p в направлении V:
omega:x*y*z*z;
p:[2,3,-1];
V:[1,2,3];
p2:p+t*V;
sp:[x=p[1],y=p[2],z=p[3]];
sp2:[x=p2[1],y=p2[2],z=p2[3]];
deltaomega:subst(sp2,omega)-subst(sp,omega);...
422 просмотров
schedule
20.08.2022
Как передать массив в качестве аргумента функции maxima?
Я работаю над созданием функций максимума, чтобы упростить оператор del для векторов. Как я могу передать список/вектор функции в максимуме? Это работает:
(%i7) dot(a,b) := a[1]*b[1]+a[2]*b[2]+a[3]*b[3];
(%o7) dot(a, b) := a b...
1126 просмотров
schedule
16.08.2022
Как избежать конфликта имен в Maxima
Я учусь, читая руководство Maxima.
Maxima имеет множество встроенных функций и глобальных переменных. Меня беспокоит возможность перезаписи функции или переменной с тем же именем, которая уже существует.
Есть ли удобный способ избежать...
98 просмотров
schedule
21.03.2023
Maxima: как избежать отрицательного символа перед выражением
Вот мой тестовый код:
kill(all)$
assume(beta>0 and gamma>0)$
simpsum:true$
assume(beta<1 and gamma<=1)$
f2:1-sum((1-beta)^l*(1-gamma)^(l-1)*gamma,l,1,inf),ratsimp;
Результат:
На самом деле перед уравнением стоит знак...
487 просмотров
schedule
27.09.2022
Максима: как преобразовать выражение в ожидаемую форму
У меня есть выражение:
(b+2*ab+a+1)/c
Я хочу использовать Maxima для факторизации уравнения, рассматривающего (b+1) как фактор.
то есть я хочу выражение в следующем виде:
[(b+1)(1+a)+ab]/c
Любая помощь будет оценена по достоинству.
831 просмотров
schedule
30.03.2024
График Максима в цикле, необходимо закрыть один график, чтобы увидеть следующий.
График Максима в цикле должен закрыть один график, чтобы увидеть следующий. Я работаю в среде Windows, но в Linux функция будет отображаться в одном представлении.
Например, эта функция:
for d:0.1 thru 1 step 0.1 do...
1348 просмотров
schedule
19.04.2024
Поля направлений, дифференциальные уравнения и решение на одном графике
Предположим, что у меня есть следующее дифференциальное уравнение:
\dot{y} = a*(y-0.5) + b*(y-0.5)^3
Мне любопытно посмотреть, можно ли изобразить на одной диаграмме фактическое дифференциальное уравнение (как указано выше) с точками для y,...
328 просмотров
schedule
18.04.2022
acos() Максимы не работает должным образом
Maxima5.27.0 дает мне результат с аргументом sqrt(3)/2
acos(sqrt(3)/2);
%pi/6
Это хорошо. Однако с аргументом 3/(sqrt(2)*sqrt(6)) он отказывается от вычислений, даже если эти два аргумента эквивалентны.
Почему и как сделать так, чтобы...
131 просмотров
schedule
13.08.2023
wxMaxima: Как постоянно добавлять пробел к знаку минус
Как сообщается в комментарии к ответу на Maxima: Как избежать отрицательного символа перед выражением , добавление строки (setf (get 'mminus 'wxxmlsym) '("<v>- </v>")) в файл maxima-init.lisp должно выполняться при запуске каждого...
108 просмотров
schedule
25.01.2023
Принудительная переоценка переменной maxima
Это может быть очень просто, но я не могу найти способ заставить Maxima сказать мне, что c теперь равно 8, а не 3? Кто-нибудь может помочь?
(%i1) a:1;
(%o1) 1
(%i2) b:2;
(%o2) 2...
205 просмотров
schedule
05.10.2023