f(z):=z^3-z+5;           /* 関数の定義 */
g(z):=''(diff(f(z),z));  /* 導関数を新関数として定義 */        
eps:1.0*10^(-5);        /* 終了条件 */
z:1.0+1.0*%i;              /* 適当な初期値を設定 */
do (w:rectform(f(z)/g(z)),z:rectform(ev(z-w)),     /* 無限ループ */
    if float(cabs(w))<eps then return(z));  /* 抜け出す */
print(z);
