MATLAB을 이용한 알기 쉬운 수치해석) 2장 7번

% x - 1 = 1 이므로,
zero = func(1);
first = func(1) + first_diff_func(1);
second = func(1) +  first_diff_func(1) + second_diff_func(1)/2;
third = func(1) +  first_diff_func(1) + second_diff_func(1)/2 + 25;
rel = func(2);

first_error = abs((rel-first)/rel*100);
second_error = abs((rel-second)/rel*100);
third_error = abs((rel-third)/rel*100);

fprintf("%f %f %f %f %f \n", zero, first, second, third, rel);
fprintf("0차 근사값 : %f %% \n", abs((rel-zero)/rel*100));
fprintf("1차 근사값 : %f %% \n", abs((rel-first)/rel*100));
fprintf("2차 근사값 : %f %% \n", abs((rel-second)/rel*100));
fprintf("3차 근사값 : %f %% \n", abs((rel-third)/rel*100));

function f = func(x)
f = 25 * x^3 - 6*x^2 + 7*x -88;
end

function f = first_diff_func(x)
f = 25 * 3 * x^2 - 6 * 2 * x + 7;
end

function f = second_diff_func(x)
f = 25 * 3 * 2 * x - 6 * 2;
end