找回密码
 注册
Simdroid-非首页
查看: 126|回复: 6

就一个简单的代入而已,maple的结果不对?【已经解决,是楼主自己看错了】

[复制链接]
发表于 2009-9-14 23:10:14 | 显示全部楼层 |阅读模式 来自 黑龙江哈尔滨
本帖最后由 TBE_Legend 于 2009-9-15 12:49 编辑

maple程序:

  1. restart;
  2. k1:=t0*g0^(-1/n1);  
  3. gama1:=(tao/k1)^n1;
  4. gama2:=((tao-t0)/k2)^n2+g0;
  5. GAMA1:=-subs({tao=-tao,t0=-t0,g0=-g0},gama1);
  6. GAMA2:=-subs({tao=-tao,t0=-t0,g0=-g0},gama2);
复制代码
mmtc程序:

  1. k1=t0*g0^(-1/n1)
  2. gama1=(tao/k1)^n1
  3. gama2=((tao-t0)/k2)^n2+g0
  4. GAMA1=gama1/.{tao->-tao,t0->-t0,g0->-g0}
  5. GAMA2=gama2/.{tao->-tao,t0->-t0,g0->-g0}
复制代码
各自的结果见图片,mmtc的结果应该是正确的。但为什么maple的结果是错误的呢?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2009-9-15 00:11:58 | 显示全部楼层 来自 湖北武汉
Simdroid开发平台
本帖最后由 maplelab 于 2009-9-15 00:13 编辑

1# TBE_Legend

应该用eval做数学代换,而不是用subs做语法代换。
这并不是什么bug。
回复 不支持

使用道具 举报

 楼主| 发表于 2009-9-15 01:33:07 | 显示全部楼层 来自 黑龙江哈尔滨
1# TBE_Legend

应该用eval做数学代换,而不是用subs做语法代换。
这并不是什么bug。
maplelab 发表于 2009-9-15 00:11


发问时,我还是比较负责的,我看过help中的subs何eval的区别,我也试过本题用eval,但结果仍然有问题。

手头机器没maple,明早到办公室时,再试试。

PS:如果替换中只用一个替换 “tao=-tao” 结果就对了。
回复 不支持

使用道具 举报

发表于 2009-9-15 09:41:12 | 显示全部楼层 来自 浙江杭州
我看了几遍,也没看出来哪错了 :L
回复 不支持

使用道具 举报

发表于 2009-9-15 12:29:39 | 显示全部楼层 来自 湖北武汉
我用eval昨晚就试过了,没有问题.
回复 不支持

使用道具 举报

 楼主| 发表于 2009-9-15 12:37:47 | 显示全部楼层 来自 黑龙江哈尔滨
我用eval昨晚就试过了,没有问题.
maplelab 发表于 2009-9-15 12:29


哦? 如果愿意的话贴下你的代码吧。

下面是我的:

eval的代码:

  1. restart; k1 := t0*g0^(-1/n1); gama1 := (tao/k1)^n1; gama2 := ((tao-t0)/k2)^n2+g0; GAMA1 := -(eval(gama1, {g0 = -g0, t0 = -t0, tao = -tao})); GAMA2 := -(eval(gama2, {g0 = -g0, t0 = -t0, tao = -tao}))
复制代码


能给出正确结果的代码:

  1. restart; k1 := t0*g0^(-1/n1); gama1 := (tao/k1)^n1; gama2 := ((tao-t0)/k2)^n2+g0; GAMA1 := -subs({tao = -tao}, gama1); GAMA2 := -subs({tao = -tao}, gama2)
复制代码

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 不支持

使用道具 举报

 楼主| 发表于 2009-9-15 12:46:45 | 显示全部楼层 来自 黑龙江哈尔滨
好的,搞定了,结果都是对的,我看错误了。 不好意思,各位!:L
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-20 00:05 , Processed in 0.051831 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表