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

%main.m

radian = newton("func", "diff_func", pi/4, 0.01);
result = radian * 180 / pi;
disp(result);

 

%func.m

function f = func(x)
f  = (0.45*20*9.81)/(cos(x)+0.45*sin(x)) - 92;
end

 

%diff_func.m

function f = diff_func(x)
f  = -(0.45*20*9.81)*(-sin(x)+0.45*cos(x))/(cos(x)+0.45*sin(x))^2;
end