수치해석 hw1 - 매트랩이용

소개글

수치해석 hw1 - 매트랩이용

목차

1) %HW1 Bisection method code
2) %HW2 Secant method code

본문내용

1)
%HW1 Bisection method code

%1-(a)
clear all
f=@(x) -1.5*x^2+3.5*x+4.5;
xl=0; xu=10; imax=100; tol=1e-8;
if f(xl)*f(xu)>0
disp('Error : xl점과 xu점에서 같은 부호를 가집니다.')
else
disp(' i xl xu (xr)Solution f(xr) tolerance')
for i=1:imax
xr=(xl+xu)/2;
toli=(xu-xl)/xu;
fprintf('%2i %10.6f %10.6f %12.6f %12.6f %12.9f\n', i,xl,xu,xr,f(xr),toli)
if f(xr) == 0
fprintf('해석해 x=%11.6f was found',xr)
break
end
if toli<tol
break
end
if f(xl)*f(xr) < 0
xu=xr;
else
xl=xr;
end
xrr(i)=xr;
y(i)=i;
tolerance(i)=toli;
end
end

참고 자료

없음