例子：

sumsqint（1）％返回[0,1]
sumsqint（65）％返回[1,8; 4,7]

如果要考虑的话，sumsqint可以与大量合作，例如解决x^2+y^2 =（25e6）^2：
sumsqint（repmat（因子（25E6），1,2））％返回9x2阵列

如果您拥有John D'Errico的可变精度整数工具箱（2009年7月20日更新），则SUMSQINT可与可变精度整数一起使用 - 请参见下面的链接。

sumsqint（vpi（5）^100）％返回51x2 VPI数组
sumsqint（vpi（65）^100）％返回5101x2 VPI数组

该算法使用高斯整数分解。请参阅帮助以获取更多详细信息。