numerical-integration

2зной

1ответ

Я решаю систему ODE с использованием RK4. Я создаю прямой график, который, по-видимому, связан с тем, что k3_1 ограничено -3.1445e + 24. Я не понимаю, почему это ограничено. function RK4system_MNModel

2зной

1ответ

Как использовать Rcpp для численного интегрирования в C++ внутри R

Мне было интересно, как Rcpp может использоваться для выполнения численного интегрирования, вызывая C++ в R. Моя текущая настройка занимает очень много времени и очень подвержена ошибкам. Я думаю, что

-2зной

1ответ

Интеграция окружности круга с помощью Matlab

Это может быть простой вопрос. Я пытаюсь понять численное интегрирование в Matlab. Для единичной окружности, fun = @(x) sqrt(cosd(x).^2+sind(x).^2); q = integral(fun,0,360); Ответ 360.0000, в то вр

0зной

1ответ

Численное интегрирование двойного интеграла с использованием квадратуры Гаусса-Лежандра

У меня уже есть код для генерации весов и абсцисс для квадратурного правила Гаусса-Лежандра для интеграции, но мне трудно использовать его для двойных интегралов. Для создания весов и абсцисс я исполь

0зной

1ответ

Двойная интеграция с использованием квадратур Гаусса-Лежандра

Я хочу использовать квадратуру Гаусса-Лежандра для вычисления двойного интеграла. Я придумал следующий код m=5000; %generate weights and abscissas [wx,xx]=leg(-1,1,m); [wy,xy]=leg(-1,1,m); %defi

0зной

1ответ

Оптимизация python quad с использованием ctypes

Мне нужно использовать функции ctypes, чтобы сократить время работы quad в python. Вот мой оригинальный вопрос original question, но теперь я знаю, какой путь мне нужно следовать. Мне нужно выполнить

1зной

2ответ

Расчета трапециевидных интеграций от происхождения каждый момент времени в R

У меня есть кадр данных x: head(x) # time Qfr #1 1 0.004751271 #2 2 0.005405618 #3 3 0.005785781 #4 4 0.006028213 #5 5 0.006179973 #6 6 0.006263814 Я пытаюсь вычислить численное интегрирова

-3зной

1ответ

Можно ли написать программу, которая может найти нули определенной функции?

В интеграции Gauss-legendre нам нужно найти нули функции legendre, но я не могу найти способ написать код, который позволит мне это сделать? Насколько я понимаю, есть список «xi», в котором эта функци