Octave / task_10

x0 = 0;
x1 = 10;
f1 = @(x,y)(x^2 * y^3 * sin((x+y)^3));

f2_2 = @(x,u,v)(-2 * u * v + sqrt(u^3 + x^2));
f2_1 = @(x,u,v) (v);

figure;
hold on
[xrk3, yrk3] = task_10_r3(f1, [x0, x1], 1);
[xrk4, yrk4] = task_10_r4(f1, [x0, x1], 1);
[xod45, yod45] = ode45(f1, [x0, x1], 1);
plot(xrk3, yrk3, 'r');
plot(xrk4, yrk4, 'g');
plot(xod45, yod45, 'b');
xlabel('X')
ylabel('Y')
grid on
legend('Метод Рунге-Кутта 3 порядка',
'Метод Рунге-Кутта 4 порядка',
'ode45');

F = @(x, y)[y(2); -2 * y(1) * y(2) + sqrt(y(1)^3 + x^2)];
figure;
hold on
[xrks3, yrks3] = task_10_r3_system(f2_1, f2_2, [x0, x1], [0.3,1]);
[xrks4, yrks4] = task_10_r4_system(f2_1, f2_2, [x0, x1], [0.3,1]);
[T, Y] = ode45(F, [x0, x1], [0.3,1]);
plot(xrks3, yrks3(:, 1), '-r');
plot(xrks3, yrks3(:, 2), '-g');
plot(xrks4, yrks4(:, 1), '-b');
plot(xrks4, yrks4(:, 2), '-y');
plot(T, Y(:, 1), '-m');
plot(T, Y(:, 2), '-c');
xlabel('X')
ylabel('Y')
grid on
legend('Метод Рунге-Кутта 3 порядка (1)',
'Метод Рунге-Кутта 3 порядка (2)',
'Метод Рунге-Кутта 4 порядка (1)',
'Метод Рунге-Кутта 4 порядка (2)',
'ode45 (1)',
'ode45 (2)');