找回密码
 快速注册
搜索
查看: 120|回复: 4

这个程序为什么只能正确运行一次?

[复制链接]

119

主题

446

回帖

3179

积分

积分
3179

显示全部楼层

TSC999 发表于 2022-5-8 20:08 |阅读模式
本帖最后由 TSC999 于 2022-5-17 20:58 编辑 下面是一个证明两线段平行的简单程序,程序是用 mathematica 写的。

程序中 g 和 h 是两个自由变量,其余各点的复数坐标都是 g 、h 的函数。此程序中需要求出 g 的共轭复数表达式,

为此由已知条件 BE // CF 和 EF 平行于实轴列出一个方程,解这个方程可以得到 的表达式,并且首次运行程序可以

得到正确的结果。但是不知什么原因,再次运行程序时就出错了。要想避免这个毛病,只能把首次运行结果补充到

程序中去,同时把解方程语句注销掉。但是这办法显得笨笨的,也让读者不解。

有什么办法解决?

今天能上传图片了。

证明.png

解方程这一段程序有问题,但是哪里出了毛病?找不出来:

  1. Clear["Global`*"];
  2. \!\(\*OverscriptBox[\(b\), \(_\)]\) = b = 0; \!\(\*OverscriptBox[\(c\), \(_\)]\) = c = 1; e = (g (1/2 - I h))/\!\(\*OverscriptBox[\(g\), \(_\)]\);
  3. \!\(\*OverscriptBox[\(e\), \(_\)]\) = ((I h + 1/2) \!\(\*OverscriptBox[\(g\), \(_\)]\))/g;
  4. f = (2 I h g + g - 2 I h - 2 \!\(\*OverscriptBox[\(g\), \(_\)]\) + 1)/(2 - 2 \!\(\*OverscriptBox[\(g\), \(_\)]\));
  5. \!\(\*OverscriptBox[\(f\), \(_\)]\) = (-2 g + 2 I h - 2 I h \!\(\*OverscriptBox[\(g\), \(_\)]\) + \!\(\*OverscriptBox[\(g\), \(_\)]\) + 1)/(2 - 2 g);
  6. Simplify@Solve[{(e - b)/(\!\(\*OverscriptBox[\(e\), \(_\)]\) - \!\(\*OverscriptBox[\(b\), \(_\)]\)) == (f - c)/(\!\(\*OverscriptBox[\(f\), \(_\)]\) - \!\(\*OverscriptBox[\(c\), \(_\)]\)), (e - f)/(\!\(\*OverscriptBox[\(e\), \(_\)]\) - \!\(\*OverscriptBox[\(f\), \(_\)]\)) == 1}, {\!\(\*OverscriptBox[\(g\), \(_\)]\)}]
复制代码

730

主题

1万

回帖

9万

积分

积分
93593
QQ

显示全部楼层

kuing 发表于 2022-5-17 21:51
复制你最后给出的这段程序,运行没问题呀,输出

  1. \!\(\*OverscriptBox["g", "_"]\) -> (g - 2 I g h)/(-1 + 2 g - 2 I h)
复制代码


运行多次都一样。

点评

噢,那有可能是我的 mathematica 版本低了? 我是 9.0 版的。  发表于 2022-5-18 10:25

730

主题

1万

回帖

9万

积分

积分
93593
QQ

显示全部楼层

kuing 发表于 2022-5-18 14:09
我的 mathematica 版本更低,我还在用 7.0

119

主题

446

回帖

3179

积分

积分
3179

显示全部楼层

 楼主| TSC999 发表于 2022-5-25 09:06
kuing 发表于 2022-5-18 14:09
我的 mathematica 版本更低,我还在用 7.0

问题找出来了。这一段小程序如果单独运行,多次运行都正常。但是如果之前运行了一个有关 g 的有毛病的程序,这段小程序就 “中毒”  了,只能运行一次。按说程序前面都加了一个删除全部变量的指令 Clear["Global`*"];,相当于打了新冠病毒防疫针,但是它还是感染了 “奥密克戎”,不知道为啥这条 Clear["Global`*"]; 指令没有起作用?

手机版|悠闲数学娱乐论坛(第3版)

GMT+8, 2025-3-4 11:55

Powered by Discuz!

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