Pari/GP 编程例子及注释

elim

打开 Pari/GP 照上面截图逐行键入指令，或者拷贝一下代码一次性粘贴到 Pari/GP:

1. \p 100
   
2. s(x)={ my(s=0.); my(t=x); my(k=1);
   
3.                   while((abs(t) > 10^(-100))&&(s=s+t)&&(k=k+2), t= -t*x^2/(k*(k-1)));
   
4.                   return(s);}
   
5. c(x)=s(Pi/2-x);
   
6. ac(x)={ my(y=sqrt(2*(1-x)));
   
7.                    while((y1=y+(c(y)-x)/s(y))&&(abs(y1-y) > 10^(-100)), y=y1);
   
8.                    return(y);}
   
9. as(x)=Pi/2-ac(x);

复制代码

命令提示符 <09:28> gp > 表示待命时刻
\\ 是注释行前缀
\p 100 告诉系统计算结果要有 100 位有效数字. 100 可换成其他正整数，上限大得吓人.
Pi 是Pari/GB 的\(\pi\) 值. 其精度由上述 \p 决定. Pari/GP 的缺省精度有计算机 bit 决定
参见手册.
函数定义在分隔符 {} 之间可占多行，以便阅读。
用 %k 表示第 k 个用户命令的计算结果. 用 % 表示上一个用户命令的计算结果.