Люди з обчислювальної науки:
Я спочатку опублікував це питання на біржі Math Stack, і хтось прокоментував, що я можу отримати "набагато кращі" відповіді тут:
Я новачок у числових методах та Matlab. Я намагаюся оцінити наступну суму двох потрійних інтегралів (це, очевидно, можна записати простіше, але ви все ще не можете його оцінити символічно (?)). У мене виникають проблеми з отриманням працюю тут, тому я неохоче розбив його на шматки: я хочу знайти суму
і
де
і
EDIT (2 березня 2013 р.): Хтось відповів, що вони мають Mathematica робити інтеграли символічно. Я просто спробував це зробити (за спрощеними версіями інтегралів), і Mathematica міг зробити лише зовнішні два з перших, і зупинився на другому. Буду вдячний за допомогу. Ось що я зробив .:
Я намагався оцінити
via
Інтегруйте [r1 ^ 3 / r2 ^ 3 * t ^ 2 * Exp (-t ^ 2), {t, 0, r2 - r1}, {r1, 1, r2}, {r2, 1, 2}]
і Mathematica повертається (у мене були проблеми з тут, оскільки результат довгий. Я розбив його на два рівняння. якщо хтось знає хороший спосіб відобразити це, будь ласка, скажіть мені:
Потім я спробував оцінити
використовуючи
Інтегруйте [(r1 + r2 - t) ^ 4 * (t ^ 2 + 2 * t * (r1 + r2) - 3 * (r2 - r1) ^ 2) ^ 2 * Exp [-t ^ 2] / r1 ^ 3 / r2 ^ 3, {r2, 1, 2}, {r1, 1, r2}, {t, r2-r1, r2 + r1}]
саме зараз, і Mathematica не повернув відповіді приблизно через півгодини (але у мене зараз проблеми з комп'ютерною мережею, і вони можуть бути винні).
[Кінець 2 березня редагувати]
Я знаю, що Matlab - це гарне програмне забезпечення, але я чув, що чисельні потрійні інтеграли важко точно зробити, і математики повинні бути скептичними, тому я хочу певним чином перевірити точність цієї відповіді. Інтеграли дають очікувану цінність певного експерименту (якщо хто хоче, я можу відредагувати це запитання, щоб описати експеримент): я реалізував експеримент у Matlab, використовуючи відповідні випадково згенеровані числа, мільйон разів, і усереднював результати. Я повторив цей процес чотири рази. Ось результати (прошу вибачення, якщо я неправильно використав слово "випробування"):
Хоча для кожного випробування було використано мільйон вибірок, значення моделювання узгоджуються лише у першій значній цифрі. Вони для мене недостатньо близькі один до одного, щоб визначити, чи чисельний потрійний інтеграл точний.
Тож чи може хто-небудь сказати мені, чи можу я довіряти результату "трійки" тут і за яких обставин можна довіряти йому взагалі?
Одним із пропозицій, який я отримав на Math Stack Exchange, було спробувати інше програмне забезпечення, наприклад, Mathematica, Octave, Maple та SciPy. Це хороша порада? Чи справді люди виконують чисельну роботу в Mathematica та Maple? Octave - це клон Матлаба, тож чи можу я припустити, що він використовує ті ж алгоритми інтеграції? Я навіть не чув про SciPy раніше і буду вдячний за будь-яку думку з цього приводу.