Страница 1 из 1

интеграл от решения дифференциальных уравнений

Добавлено: Вт июл 17, 2018 7:13 pm
fenix85
Подскажите пожалуйста, как получить интеграл от численного решения диф. уравнений. У меня есть решение - функция Z6(x), нужно ее проинтегрировать на интервале и получить график этого интеграла. Я не знаю как это сделать...(((

Re: интеграл от решения дифференциальных уравнений

Добавлено: Чт июл 19, 2018 9:21 pm
Markiyan Hirnyk
Добавьте в конце вашего кода команду

Код: Выделить всё

plot(x->int(Z6(s), s = 0 .. x, numeric, epsilon = 0.1e-1) , 0 .. tau, numpoints = 5, colour = [black], thickness = [4, 4], labeldirections = [horizontal, vertical], labels = [s*M, `dθ`/ds]);

Re: интеграл от решения дифференциальных уравнений

Добавлено: Пт июл 20, 2018 3:03 pm
fenix85
Зделал и получил вот такое сообщение Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct
Смотрел на странице документации по этой ошибке но так ничего и не понял(
У меня Maple 13.

Re: интеграл от решения дифференциальных уравнений

Добавлено: Пт июл 20, 2018 3:21 pm
Markiyan Hirnyk
Пожалуйста, представьте через Dropbox или иной файлообменник mw file с выполненным кодом. Возможно также, что вы не выполнили все ваши предыдущие команды. Конечно же, программа должна начинаться с команды restart.

Re: интеграл от решения дифференциальных уравнений

Добавлено: Сб июл 21, 2018 7:43 am
Markiyan Hirnyk
fenix85, спасибо за файл, переданный через Гугл. Действительно, при выполнении команды в Мэйпле 16.02 ( версии 13 не имею) получается та же ошибка. Работает в Мэйпле 2018.1.Возможно, изменился синтаксис Мэйпла. Подумаю над этим.

Re: интеграл от решения дифференциальных уравнений

Добавлено: Сб июл 21, 2018 8:11 am
Markiyan Hirnyk
Команда

Код: Выделить всё

plot(evalf(Int('Z6(s)', s = 0 .. x)), x = 0 .. tau, numpoints = 5);
работает в Мэйпле 16.02, но медленно: выполняется примерно 6 минут на моем компе.

Re: интеграл от решения дифференциальных уравнений

Добавлено: Сб июл 21, 2018 7:36 pm
fenix85
Благодарю за помощь! У меня тоже заработало, на моем компьютере выполняется примерно 7 минут. Правда меня удивляет что интеграл получается прямая линия, тогда как сама функция осциллирует.

Re: интеграл от решения дифференциальных уравнений

Добавлено: Сб июл 21, 2018 8:25 pm
Markiyan Hirnyk
Амплитуда колебаний порядка 10^(-7) около среднего значения 0.14e-4, поэтому график первообразной и близок к прямой линии.

Re: интеграл от решения дифференциальных уравнений

Добавлено: Сб июл 21, 2018 10:28 pm
Markiyan Hirnyk
Еще один подход изложен здесь. Идея состоитв следующем: добавить

Код: Выделить всё

diff(z9(x), x) = z6(x)
в , в и в c понятными дальнейшими изменениями. Проверял, работает.

Re: интеграл от решения дифференциальных уравнений

Добавлено: Пн июл 23, 2018 8:23 pm
fenix85
спасибо большое!