[Octave]250318_RYU_poly,dff

clc; clear;

r = [3, 2, -5];

p = poly(r);
x = linspace(-10, 10, 200);
y = polyval(p, x);


%다항식 그래프
figure(1)
plot(x, y, 'b', 'linewidth', 2);
##xlim([-10, 10]);
##ylim([-1000, 1000]);
xlabel('x');
ylabel('y');
legend('다항식의 곡선');
title('다항식');
grid on;

%미분 그래프
dp = polyder(p);
dy = polyval(dp, x);

figure(2)
plot(x, dy, 'r', 'linewidth', 2);
##xlim([-10, 10]);
##ylim([-1000, 1000]);
xlabel('x');
ylabel('dy/dx');
legend('미분의 곡선');
title('미분');
grid on;

%차분의 정의에 의한 곡선

h = x(2) - x(1);
y_delta = (y(2:end)-y(1:end-1));
y_diff = y_delta/h;
x_diff = x(1:end-1);

figure(3)
plot(x_diff, y_diff, 'g', 'linewidth', 2);
xlabel('x');
ylabel('{\Delta}y/{\Delta}x');
legend('차분의 곡선');
title('차분');
grid on;

% 차분 계산함수 diff

figure(4)
y_diff2 = diff(y)/h;
plot(x_diff, y_diff2, 'm', 'linewidth', 2);
xlabel('x');
ylabel('{\Delta}y/{\Delta}x');
legend('차분의 곡선');
title('차분');
grid on;
% 모든 곡선을 하나의 축에 표현

figure(5)
plot(x, y, 'b', 'linewidth', 2); hold on
plot(x, dy, 'r', 'linewidth', 4); hold on
plot(x_diff, y_diff, 'g', 'linewidth', 3); hold on
plot(x_diff, y_diff2, 'm', 'linewidth', 2); hold off
xlim([-10, 10]);
ylim([-1000, 1000]);
grid on;
legend('다항식', '미분', '정의에 의한 차분', 'diff함수의 차분');
title('다항식, 미분, 차분 비교');