numerical-analysis

1зной

1ответ

C++ - Назад Замена через Boost - Bad Index Error

У меня есть следующий код для алгоритма обратной подстановки: #include <boost\numeric\ublas\matrix.hpp> typedef boost::numeric::ublas::matrix<double> dM; dM bSub(dM A, dM b) { unsigned int n =

0зной

1ответ

Как передать твердое кодированное дифференциальное уравнение через Рунге-Кутту 4

Я пытаюсь реализовать Runge-Kutta, например, проблемы dy/dt = y - t^2 + 1 и dy/dt = t * y + t^3 в C#, и я не могу получить результат, который я ожидаю. Я разделил свою программу на несколько классов,

-1зной

2ответ

python числовая ошибка: большое количество добавить небольшое число

Это упражнение из углубленного курса Udacity. Может ли кто-нибудь объяснить, почему окончательный ответ не является 1.0? v1 = 1e9 v2 = 1e-6 for i in range(int(1e6)): v1 = v1 + v2 print 'answer

0зной

1ответ

NumPy конечных разностей для Лапласа PDE с полигональными границами

NumPy сценария на http://www.physics.drexel.edu/courses/Comp_Phys/Physics-306/demo6.1.py демонстрируют, как решить Лаплас PDE для ПРЯМОУГОЛЬНОГО домена . Вопрос: есть ли аналогичный URL-адрес для POLY

0зной

1ответ

Алгоритм численной экономической вычислимости

Я знаю, что эти типы вопросов обескуражены как субъективные, но после получения нескольких «менее полезных» ответов на мою более конкретную версию этого вопроса я хотел бы попробовать еще раз. Моя про

0зной

1ответ

Использовать функциональный интеграл (f), который вычисляет интеграл от f по [0,1] для вычисления интеграла от f по любому интервалу

Как использовать (гипотетический) интеграл (f), который вычисляет интеграл от функции f над [0,1], чтобы вычислить интеграл от f по любому отрезку [a, b]? Если бы я хотел изменить интервал интегрирова

3зной

1ответ

Избегайте опустошения, используя ехр и минимальное положительное float128 в Numpy

Я пытаюсь вычислить следующее соотношение: w(i)/(sum(w(j)), где w обновляются с использованием экспоненциальной функции убывающей, т.е. w(i) = w(i) * exp(-k), k является положительным параметром. Все

-1зной

1ответ

Условия первого порядка с экспонентами в Mathematica

Я попытался использовать функции «Решить» и «NSolve» в Mathematica, чтобы найти значение p, которое максимизирует (p-1) (1 + exp (up) + exp (vq)) - exp (uq) = 0, где u, v, p, q больше 0, безрезультатн