符号运算

mathe

不知道，我现在大部分功能还不会用

无心人

哦

现在要评估作为数值计算的功能是否强大

呵呵

mathe

你可以让它计算一下1000000!看看，好像挺慢的
－－－－－－－－－－－－
好像100000!已经不行了

gxqcn

原帖由 mathe 于 2008-6-19 12:00 发表
最后可以得到公式：
$a(n)=round((root{4}{54}+sqrt(6)+root{4}{24}+3)/{8 * (root{4}{3/2}-1)^n}) $
round表示4舍5入

看这个通项公式，其特征方程应该是两次左右，
而主帖中的特征方程次数很高，不知讨论的可还是同一问题？

mathe

原帖由 gxqcn 于 2008-6-19 15:24 发表


看这个通项公式，其特征方程应该是两次左右，
而主帖中的特征方程次数很高，不知讨论的可还是同一问题？

特征方程是4次，也就是$x^4-8x^3-12x^2-8x-2$,经过化简以后发现前面$(x+1)^8$是没有必要的。
而这个特征方程的4个根都可以解出来，也就是#34中的u,v,s,t,其中只有u的绝对值大于1（而且远远大于）。

无心人

倒，肯定不行的
毕竟是lisp

我想让你测数值计算阿
小规模的

gxqcn

谢谢 44#，我明白了。

这么说，其递推公式可为：$a_n = 8a_{n-1} + 12a_{n-2} + 8a_{n-3} +2a_{n-4} \quad (n>4)$

无心人

就测试下列问题吧

   求$x^{x + sqrt(x) - e^x} + sin(x)$在$[1/2 .. 1]$间的定积分

mathe

原帖由 无心人 于 2008-6-19 15:44 发表


倒，肯定不行的
毕竟是lisp

我想让你测数值计算阿
小规模的

呵呵，数值计算我还不会，现在只会用符号运算。可能使用浮点数就会自动使用数值计算吧

mathe

原帖由 gxqcn 于 2008-6-19 15:45 发表
谢谢 44#，我明白了。

这么说，其递推公式可为：$a_n = 8a_{n-1} + 12a_{n-2} + 8a_{n-3} +2a_{n-4} \quad (n>4)$

是的