탄젠트 pi/4가 1이니까  아크탄젠트함수 이용해서 원주율 구하려고하는데


function y=prob1_2_ex(x,n)

y=0;

for i=0:n


    if mod(i,2)==1

         y=y+(-1)^((i+1)/2+1)*x^i/i;

  

        

    end


end

%y=abs((y-sin(x))/sin(x))*100;

%y=abs((y-single(sin(x)))/single(sin(x)))*100;

end

clear;
for j=0:100
    rl(1,j+1)=prob1_2_ex(1,j);
    if(rl(1,j+1)==0)
        disp('n= ');
        disp(j);
        %fprintf('n=%d',j);
        break
    end
end
n=0:1:j;
semilogy(n,rl,'o','markerfacecolor','b');
set(gca,'fontsize',15);
xlabel('n[terms]');
ylabel('relative error[%]');
hold on

이렇게 돌리면


 


왜이러냐... 나는 원주율 구하고 상대오차 구하려고한건데 상대오차가 10에 수렴해버리노