RuDas의 무능월드

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

n=input('n=');x=input('x=');a=input('a='); sum=sqrt(a)+(x-a)*1/(2*sqrt(a)); mul=1;number=1;for i=2:nmul=mul*number;sum=sum+(x-a)^i/factorial(i)*(-1)^(i+1)*(1/2)^i*a^(1/2-i)*mul;number=number+2;endsumerror=abs(sum-sqrt(x))

function x=gauss_jordan2(A)% A: 계수행렬% I: 단위행렬% x: inv(A)의 계산값 [m1,k1]=size(A);% [m2,k2]=size(b); x=zeros(m1,k2);I=eye(m1); % 상삼각 행력 형성for i=1:m1-1 m=-A(i+1:m1,i)/A(i,i); A(i+1:m1,:)=A(i+1:m1,:)+m*A(i,:); I(i+1:m1,:)=I(i+1:m1,:)+m*I(i,:);end%A%I % 하삼각 행렬 형성for i=1:m1-1 for j=i+1:k1 m=-A(i,j)/A(j,j); A(i,:)=A(i,:)+m*A(j,:); I(i,:)=I(i,:)+m*I(j,:); end %A %Iend%A%Ifor i=1:m1 x(i,:)=I(i,:)/A(i,i); A..

Fr_bin2decBinary => Decimal Fr_dec2binDecima => Binary