Решение задачи Коши методом Рунге-Кутта 4 порядка в MatLab

@Maikl2020 · Регистрация: 13.04.2014

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

Здравствуйте подскажите, пожалуйста как можно вывести численные расчета на экран

Дано задача Коши $https://www.cyberforum.ru/cgi-bin/latex.cgi?y'=y(1-x);y(0)=1$ в интервале [0; 1], h=0.2

Код М-файл

Matlab M

function f=pr(x,y)
f=zeros(1,1);
f=y*(1-x);
end

Основной код

Matlab M

1	[x,y]=ode45(@pr,[0,1],[0.2]);

И еще подскажите пожалуйста что обозначает функция zeros и символ @

Зосима · 21.08.2014, 11:04

Чтобы вывести результаты расчетов на экран нужно использовать команду, например, disp:

Matlab M

[x,y] = ode45...
disp([x,y])

или просто не ставить в конце точкозапятую

zeros создает массив нулей нужного размера, в твоем случае это не обязательно, потому как ДУ первого порядка.
@ указывает, что идет обращение к функции pr
однако как я понимаю h - это шаг времени, а не начальное условие, поэтому нужно писать так:
[x,y]=ode45(@pr, [0:0.2:1], 0);

@Maikl2020 · 21.08.2014, 11:21 **[ТС]**

Зосима, Спасибо за ответ. Еще вопрос почему значения y всегда ноль показывает

Matlab M

>> [x,y]=ode45(@pr,[0:0.2:1],0);
disp([x,y])
         0         0
    0.2000         0
    0.4000         0
    0.6000         0
    0.8000         0
    1.0000         0

Зосима · 21.08.2014, 11:25

Maikl2020, попробуй взять другое начальное условие

[x,y]=ode45(@pr, [0:0.2:1], 3);

@Maikl2020 · 21.08.2014, 12:04 **[ТС]**

Зосима, спасибо!

Зосима 5236 / 3564 / 378 Регистрация: 02.04.2012 Сообщений: 6,471 Записей в блоге: 17
	21.08.2014, 11:04	2
	Чтобы вывести результаты расчетов на экран нужно использовать команду, например, disp: Matlab M [x,y] = ode45... disp([x,y]) или просто не ставить в конце точкозапятую zeros создает массив нулей нужного размера, в твоем случае это не обязательно, потому как ДУ первого порядка. @ указывает, что идет обращение к функции pr однако как я понимаю h - это шаг времени, а не начальное условие, поэтому нужно писать так: [x,y]=ode45(@pr, [0:0.2:1], 0); 0

Зосима 5236 / 3564 / 378 Регистрация: 02.04.2012 Сообщений: 6,471 Записей в блоге: 17
	21.08.2014, 11:25	4
	Maikl2020, попробуй взять другое начальное условие [x,y]=ode45(@pr, [0:0.2:1], 3); 1

@Maikl2020 Всегда онлайн 49 / 49 / 10 Регистрация: 13.04.2014 Сообщений: 1,409
	21.08.2014, 12:04 [ТС]	5
	Зосима, спасибо! 0

Опции темы