Newton method（正確版）

[sitting (臭小開吃屎)]

x=input('請輸入一整數值：');    % 輸入一整數
if x>0                           % x>0才能繼續往下做
  delta = 0.001;                    % set convergence tolerance=0.001
  r=1;rold=2;                      % 設r和rold的初始值
  it=0;                            %設定計數
  maxit=30;                        %設定最大計數次數
  while  abs(rold-r)> delta & it<maxit

      rold = r;                    %save old value for convergence test
      r =0.5*(rold + x/rold);      %牛頓近示式
      it = it + 1;
        if abs(rold-r)< delta      %當abs(rold-r)< delta 則跳出程式,不再計算
              break;
            end
    end
else
   fprintf('你輸入的數非正整數！');
  return; %當輸入負數,則return to calling function
end
fprintf('\n %d^(1/2)=%g,共計數=%d次',x,r,it);


--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.123.224.93