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

%main.m

i = newton("current", "diff_current", 1, 0.01);
disp(i);

 

%current.m

function f = current(t)
f = 9*exp(-t)*sin(2*pi*t) -3.5;
end

 

%diff_current.m

function f = diff_current(t)
f = -9*exp(-t)*sin(2*pi*t) + 9 *exp(-t) * cos(2*pi*t) * 2*pi;
end