Десять или два года назад стоило написать числовой код, чтобы избежать использования умножения и деления и вместо этого использовать сложение и вычитание. Хорошим примером является использование прямых разностей для оценки полиномиальной кривой...

Я использую C ++ для вычисления различных типов специальных функций (например, функции Ламберта, итерационных методов для оценки инверсий и т. Д.). Во многих случаях очевидно лучший подход к работе с мантиссой и экспонентой напрямую. Я нашел много...

Как преобразовать символическое выражение в функцию Octave из символического пакета? После установки символьного пакета на октаву с pkg install -forge symbolic . Используя символический пакет на октаве, я могу написать это: octave> pkg...

Я пишу программу, которая вычисляет LU-разложение матрицы с частичным поворотом, и я хотел бы, чтобы функция выводила несколько (2 или 3) матриц без запуска программы несколько раз для вывода каждой по отдельности, что является пустой тратой времени....

Я использую Clojure, и мне нужно запустить небольшую симуляцию. У меня есть вектор длины n (n обычно между 10 и 100), который содержит значения. В каждом раунде моделирования (возможно, 1000 раундов вместе) одно из значений в векторе обновляется...

Я хочу решить основной интеграл по заданному интервалу [a, b] для произвольных значений n и m в R (n, m) с интегрированием Ромберга. Я вывел правило Буля из правила трапеций, поэтому знаю, как это сделать на бумаге. Я даже нарисовал блок-схему,...

Я стараюсь сделать это коротко. У меня возникла задача численного моделирования потока в газовой трубе. Я моделирую, используя разные модели для воздушного потока (уравнения переноса), движущегося от одной границы к другой (т.е. справа налево)....

Недавно я начал заниматься численными вычислениями и решать математические задачи численно, программируя на C++ с OpenMP. Но теперь моя проблема слишком велика, и на ее решение даже распараллеленных уходит несколько дней. Итак, я думаю начать...

Проблема Я реализовал алгоритм верлета скорости для вычисления траекторий двух тел, взаимодействующих друг с другом гравитационно (только ньютоновская гравитация). Меньшее тело, вращающееся по орбите, имеет очень маленькую массу, тело,...

Я пытаюсь реализовать разделяй и властвуй SVD верхней двухдиагональной матрицы B, но мой код не работает. Ошибка: "Невозможно выполнить присваивание, так как размер левой стороны 3 на 3, а размер правой стороны 2 на 2. V_bar(1:k,1:k) = V1;"...