Метод Рунге-Кутта

@parkito · Регистрация: 22.03.2010

Студворк — интернет-сервис помощи студентам

Здравствуйте, помогите пожалуйста решить проблему.
Решаю простой дифур

Matlab M
1
dy=-2*x*y

Создаю m-файл

Matlab M
1
2
3
function dy=fu(x,y)
dy=-2*y*x^2;
end

Запускаю

Matlab M
1
[T,Y]=ode45(@fu,[0,2],1

Получаю ошибку

Matlab M
1
2
Error in ode45 (line 222)
    rh = norm(f0 ./ max(abs(y),threshold),inf) / (0.8 * rtol^pow);

Что я делаю неправильно?

@S_el · 26.02.2015, 02:22

Сообщение от parkito

Что я делаю неправильно?

Ошибки я не вижу.попробуйте запустить мой код:

Matlab M
1
2
3
A=dsolve('Dy = -2*x^2*y','y(0)==1','x')
[x,Y]=ode45(@fu,[0,2],1);
plot(x,Y,x,subs(A))

как видите численное решение полностью совпадает с аналитическим.

@parkito · 26.02.2015, 03:09 **[ТС]**

S_el, Такая же ситуация

Зосима · 26.02.2015, 17:10

parkito, странно, у меня работает

Новые блоги и статьи Все статьи Все блоги /
EmBitz, создание проекта, отладка, прошивка locm 15.07.2025 Создание проекта для Blue Pill (STM32F103C8T6) в EmBitz 2. 30, написания кода blink, запуск отладки в ОЗУ, заливка релизной прошивки во flash используя ST-Link и др. . . .	Трассировка корутин Kotlin с OpenTelemetry mobDevWorks 14.07.2025 Асинхронное программирование меняет правила игры, особенно когда речь заходит о трассировке операций. В Kotlin с его корутинами эта проблема приобретает особый оттенок, который я хотел бы детально. . .	Облачные приложения на Rust: руководство по архитектуре микросервисов golander 13.07.2025 Когда я впервые взялся за проектирование облачной платформы для одного из наших клиентов, выбор стоял между привычными Go и Java. Но после нескольких месяцев разработки микросервисной системы,. . .	Как Node.js выполняет асинхронные операции Reangularity 13.07.2025 Каждый раз, когда я рассказываю про Node. js, возникает один и тот же вопрос: "Как эта штука может быть быстрой, если JavaScript — однопоточный язык?" И это действительно кажется парадоксом. Ведь в. . .	Как писать чистый, тестируемый и качественный код на Python py-thonny 12.07.2025 Помню свой первый проект на Python. Работал тогда быстро, грязно, лишь бы работало. Код был похож на запутанный клубок - переменные по одной букве, функции на 200 строк, комментарии отсутствовали как. . .
Blazor и контроллер сервопривода IoT Meadow Maple Wired 11.07.2025 Я решил разобраться, как можно соединить современные веб-технологии с миром "железа". Интересная комбинация получилась из Blazor в качестве веб-интерфейса и микроконтроллера Meadow с его веб-сервером. . .	Генерация OpenQASM из кода Q# EggHead 10.07.2025 Летом 2024-го я начал эксперименты с библиотекой Q# Bridge, и знаете что? Она оказалась просто находкой для тех, кто работает на стыке разных квантовых экосистем. Основная фишка этой библиотеки -. . .	Изучаем новый шаблон ИИ-чата .NET AI Chat Web App stackOverflow 10.07.2025 В . NET появилось интересное обновление - новый шаблон ИИ-чата под названием . NET AI Chat Web App. Когда я впервые наткнулся на анонс этого шаблона, то сразу понял, что Microsoft наконец-то. . .	Результаты исследования от команды ARP (июль 2025 г.) Programma_Boinc 10.07.2025 Результаты исследования от команды ARP (июль 2025 г. ) Африканский проект по дождям (ARP) World Community Grid снова запущен! Мы рады поделиться обновленной информацией о нашем прогрессе с осени. . .	Angular vs Svelte - что лучше? Reangularity 09.07.2025 Сегодня рынок разделился на несколько четких категорий: тяжеловесы корпоративного уровня (Angular), гибкие универсалы (React), прогрессивные решения (Vue) и новая волна компилируемых фреймворков. . .

@parkito 15 / 14 / 10 Регистрация: 22.03.2010 Сообщений: 695
	26.02.2015, 03:09 [ТС]
	S_el, Такая же ситуация 0

Зосима 5243 / 3571 / 379 Регистрация: 02.04.2012 Сообщений: 6,475 Записей в блоге: 17
	26.02.2015, 17:10
	parkito, странно, у меня работает 0

Опции темы

Метод Рунге-Кутта

Решение