Решение дифференциального уравнения второго порядка методом Рунге-Кутта

30.05.2013, 00:36. **Показов** 13289. **Ответов** 1

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

Доброй ночи, форумчани! Помогите наисать программу в матлабе для решения следующего уравнения:

x²y''-xy'=0

y(1)=0; y'(1)=1; X принадлежит [1,2]

Зосима · 30.05.2013, 10:08

Это не сложно

Matlab M

function diffsolve
clc
Tx = [1 2];
y0 = [0 1];
[x y] = ode45(@diffs,Tx,y0); % услиенно решаем методом Рунге-Кутта 4-го порядка
plot(x, y) 
grid on,
legend('y(x)','dy/dx')
xlabel('x')
end
 
% система 
function dydx = diffs(x,y)
% y(1) -> y(x)
% y(2) -> y'(x)
% dydx(1) -> y'(x)
% dydx(2) -> y''(x)
dydx = zeros(2,1); % создаем вектор-столбец
dydx(1) = y(2);
dydx(2) = ( x*y(2) )/x^2;
end

ferado
		1
	Решение дифференциального уравнения второго порядка методом Рунге-Кутта 30.05.2013, 00:36. Показов 13289. Ответов 1 Метки нет (Все метки) Доброй ночи, форумчани! Помогите наисать программу в матлабе для решения следующего уравнения: x²y''-xy'=0 y(1)=0; y'(1)=1; X принадлежит [1,2]

Опции темы