稍修改下我的Mathematica代码，关于高斯整数开平方

wsc810

1. d = 11 - 4 I; pell = -1; P[0] = 0; Q[0] = 1;
   
2. x[0] = (P[0] + Sqrt[d])/Q[0];
   
3. a[0] = Round[x[0]];
   
4. i = 0;
   
5. While[(x[i] != 1/(x[0] - a[0]) && P[i] != pell) || i == 1,
   
6.   P[i + 1] = Q[i] a[i] - P[i];
   
7.   Q[i + 1] = (d - P[i + 1]^2)/Q[i];
   
8.   x[i + 1] = (P[i + 1] + Sqrt[d])/Q[i + 1];
   
9.   a[i + 1] = Round[x[i + 1]];
   
10.   Print[{i, P[i], Q[i], a[i]}]; i++];

复制代码

其中pell为事先要迭代停止的P值，等于-1则不用管它，P[0]为P的初值，这里等于0也不用管它，现在的有两个问题，
一, 怎样修改代码，使得 P[ i ],Q[ i ],a [ i ]分别按照一维表的格式输出。
二, 怎样修改代码，在复平面上分别画出P[ i ],Q[ i ],a[ i ]点值分布图。
Mathematica我还不太熟悉，Any help will be appreciated !
程序运行的结果可以自己看

zgg___

呵呵，没有明白什么叫做高斯整数开平方呢。
所以弱弱的问一下，如果输入的是11-4I，那么要求输出的结果是什么呢？

郭先抢

我猜他是为了求解pell方程而写的这个程序，不过我觉得没什么意思！
况且与其修改他写的程序，还不如我自己写一个呢

wsc810

不错，是求解pell方程的函数，不过该程序是暂时先得到高斯整数的不尽根的连分式展开。

wsc810

我想把把它用于大数分解，想动态利用P，Q，a的值，即即每次迭代后，不保留原来的值。用什么语句,帮帮我啊