수치해석,테일러급수,매트랩코드

소개글

2011년 과제로 직접 짠 매트랩 코드입니다. 이상없이 잘 돌아갑니다.

목차

없음

본문내용

clear
clc
%이분법 풀이한 결과와 비교하기 위해 만들어보았습니다.
f=`x^6+exp(3*x)+3*x^3`
n=0;a=-1;b=1;
err=0.00001;
c=(a+b)/2;
while b-c>=err %지정한 오차범위 안에 들 경우 STOP
n=n+1;
x=b; fb=eval(f);
x=c; fc=eval(f);
if fb*fc<=0
a=c; c=(a+b)/2;
else
b=c; c=(a+b)/2;
end
fprintf(`%3.0f번 a:%5.6f b:%5.6f c:%5.6f 오차:%5.6e \n`,n,a,b,c,b-c)
end

<중 략>

clear
clc

f=`x^6+exp(3*x)+3*x^3`

n=0;a=-1;b=1;
err=0.00001;
x=a; fa=eval(f);
x=b; fb=eval(f);
c=b-fb*((a-b)/(fa-fb));
x=c; fc=eval(f);
while abs(fc)>=err
n=n+1;
x=a; fa=eval(f);
x=b; fb=eval(f);
c=b-fb*((a-b)/(fa-fb));
x=c; fc=eval(f);
if fa*fc<=0
b=c; c=b-fb*((a-b)/(fa-fb));
else
a=c; c=b-fb*((a-b)/(fa-fb));
end
fprintf(`%3.0f번 a:%5.6f b:%5.6f c:%5.6f 오차:%5.6e \n`,n,a,b,c,b-c)
end

참고 자료

없음