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

wsc810 · 发表于 2012-11-7 08:58:18

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？欢迎注册

×

d = 11 - 4 I; pell = -1; P[0] = 0; Q[0] = 1;
x[0] = (P[0] + Sqrt[d])/Q[0];
a[0] = Round[x[0]];
i = 0;
While[(x[i] != 1/(x[0] - a[0]) && P[i] != pell) || i == 1,
P[i + 1] = Q[i] a[i] - P[i];
Q[i + 1] = (d - P[i + 1]^2)/Q[i];
x[i + 1] = (P[i + 1] + Sqrt[d])/Q[i + 1];
a[i + 1] = Round[x[i + 1]];
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___ · 发表于 2012-11-7 10:52:40

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

郭先抢 · 发表于 2012-11-7 15:16:27

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

wsc810 · 发表于 2012-11-7 15:38:50

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

wsc810 · 发表于 2012-11-20 13:18:22

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

账号		自动登录	找回密码
密码			欢迎注册

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。