RuDas의 무능월드

function inter=newton_diff(x, y, new_x)% x,y: 입력데이터% new_x: 보간이 적용되는 x값 n=length(x); b=zeros(n,n); b(:,1)=y(:); %분할차분표 작성for i=2:n for j=1:n-i+1 b(j,i)=(b(j+1,i-1)-b(j,i-1))/(x(i+j-1)-x(j)); endend b%뉴턴보간다항식 계산product=1;inter=b(1,1);for j=1:n-1 product=product*(new_x-x(j)); inter=inter+product*b(1,j+1);end

function b=newtondd(x, y)% x,y: 입력데이터 n=length(x); b=zeros(n,n); b(:,1)=y(:); %분할차분표 작성for i=2:n for j=1:n-i+1 b(j,i)=(b(j+1,i-1)-b(j,i-1))/(x(i+j-1)-x(j)); endend

function [inter,error]=lagrange_interpol(x, y, new_x)% x,y: 주어진 좌표점% new x_ : 보간이 필요한 x값 n=length(x);sum=0;for i=1:n temp=y(i); %y(i)값 temp로 저장for j=1:n %y(i)*Li(x)계산 if(i~=j) temp=temp*(new_x-x(j))/(x(i)-x(j)); endendsum=sum+temp; %Pn(x)계산end inter=sum; %Pn(x)값 출력error=abs(tan(1.15)-inter)

function root=falseposition(func, x1, x2, threshold)% x1, x2 : 초기경계값% threshold: 반복여부를 결정하는 임계값% root: 가위치법에 의해 최종적으로 구해진 근의 근사값 if feval(func, x1)*feval(func, x2) > 0 % feval(func,x)는 함수값을 계산하는 함수 disp('해는 이 구간에 없습니다.') returnend %rel_error=100.0; % 상대오차 초기값%x3=x1;j=0;fprintf('| iteration | x | error |\n')while(1) j=j+1;% xold=x3; x3=x2-feval(func, x2)*(x1-x2)/(feval(func, x1) - feval(func, x2)..

function root=secent(func, x1, x2, threshold)fprintf('| iteration | x | error |\n')j=0;while(1) j=j+1; xn=x2-feval(func, x2)*(x2-x1)/(feval(func, x2)-feval(func,x1)); if xn ~=0 rel_error=abs((xn-x1)/xn)*100; end x1=x2; x2=xn; fprintf('| %d |% 10.5f | %9.5f |\n',j,x2,rel_error) if(rel_error < threshold) break endendroot=x2;

function root=newton2(func1,func2,x0, threshold)%func1: input function%func2: derivative of input%func2=diff('func1',1);fprintf('| iteration | x | error |\n')j=0;while(1) j=j+1; x=x0 - feval(func1, x0)/feval(func2, x0); rel_error=abs((x-x0)/x); x0=x; fprintf('| %d |% 10.5f | %9.5f |\n',j,x0,rel_error) if(rel_error < threshold) break endendroot=x0;

function root = bisection(func, x1, x2, threshold)%func = 함수%x1, x2 = 범위%threshold = 오차율 if feval(func, x1)*feval(func, x2) > 0 % feval(func,x) disp('해는 이 구간에 없습니다..') returnend j = 0;fprintf('| iteration | x | error |\n') while(1) % begin of while j=j+1; %iteration번호 x3=(x1+x2)/2; % new x3 계산 if feval(func, x2)*feval(func, x3) < 0 xold=x1; x1=x3; else xold=x2; x2=x3; end if x3 ~= 0 rel_error=ab..

n = input('n=');x = input('x='); sum = 1;disp('iteration | sum | Error'); for i = 1 : n sum = sum+(x^(i))/factorial(i); error = exp(x) - sum; fprintf('%d |%f |%f\n', i, sum, error);end

n = input('n=');x = input('x='); sum = 1;disp('iteration | F(x) 근사값 | Error');for i = 1 : n-1 sum = sum + ((-1)^(i) * x^(2*i))/factorial(2*i); error = abs(cos(x) - sum); fprintf('%d | %f |%f\n', i, sum, error);end

n = input('n=');x = input('x='); sum = 0;disp('iteration | F(x) 근사값 | Error');for i = 1 : n sum = sum + ((-1)^(i+1) * x^(2*i-1))/factorial(2*i-1); error = abs(sin(x) - sum); fprintf('%d | %f |%f\n', i, sum, error);end