f(z):=z^3-z+5;              /* 関数の定義 */
df(z):=''(diff(f(z),z));    /* 導関数 */        
pi:3.14;                    /* 定数の定義 */
N:20;
s1:0;                       /* 和の初期値 */
t:-pi/2;
dt:pi/N;                 /* 積分変数の増分 */
z1:float(2*cos(t)+2*%i*sin(t));
for j: 1 thru N
  do (w:rectform(df(z1)/f(z1)),t:t+dt,z2:float(2*cos(t)+2*%i*sin(t)),
      s1:rectform(s1+w*(z2-z1)), z1:z2);
s2:0;
z1:2*%i;
dz:4.0*%i/N;
for j: 1 thru N
do (w:rectform(df(z1)/f(z1)),z2:z1-dz,
    s2:rectform(s2+w*(z2-z1)), z1:z2); 
s:rectform((s1+s2)/2/pi/%i);  /* 実部が根の個数の近似値となる */
