数学研发网设为首页收藏本站
切换到窄版

 找回密码
 欢迎注册
数学研发论坛»论坛 › 【算法研发】 › 编程擂台 › Mathematica7.0问题在哪里
返回列表 发新帖
查看: 6186|回复: 8

[讨论] Mathematica7.0问题在哪里

[复制链接]
dlsh
发表于 2021-11-22 19:53:45 | 显示全部楼层 |阅读模式

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

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

×

第二次运行出错

第二次运行出错
为什么第二次出现提示？很长时间运行不出结果？第一次很快出来？

第一次正确

第一次正确

  2. \!\(\*OverscriptBox["i", "_"]\) = i = 0;
  3. \!\(\*OverscriptBox["m", "_"]\) = 1/m;
  4. a = (2 p m)/(m + p);
  5. \!\(\*OverscriptBox["a", "_"]\) = 2/(m + p); b = (2 p n)/(n + p);
  6. \!\(\*OverscriptBox["b", "_"]\) = 2/(n + p); c = (2 n m)/(n + m);
  7. \!\(\*OverscriptBox["c", "_"]\) = 2/(m + n);
  8. d = (a + c)/2;
  9. \!\(\*OverscriptBox["d", "_"]\) = (
  10. \!\(\*OverscriptBox["a", "_"]\) +
  11. \!\(\*OverscriptBox["c", "_"]\))/2; f = (a - \[Lambda] d)/(
  12. 1 - \[Lambda]);
  13. \!\(\*OverscriptBox["f", "_"]\) = (
  14. \!\(\*OverscriptBox["a", "_"]\) - \[Lambda]
  15. \!\(\*OverscriptBox["d", "_"]\))/(1 - \[Lambda]);(*假设
  16. \!\(\*OverscriptBox["FA", "\[RightVector]"]\)：F
  17. \!\(\*OverscriptBox["D", "\[RightVector]"]\)=\[Lambda]*)

  19. k[a_, b_] := (a - b)/(
  20. \!\(\*OverscriptBox["a", "_"]\) -
  21. \!\(\*OverscriptBox["b", "_"]\));
  22. \!\(\*OverscriptBox["k", "_"]\)[a_, b_] := 1/k[a, b];(*复斜率定义*)

  24. \!\(\*OverscriptBox["Jd", "_"]\)[k1_, a1_, k2_, a2_] := -((a1 - k1
  25. \!\(\*OverscriptBox["a1", "_"]\) - (a2 - k2
  26. \!\(\*OverscriptBox["a2", "_"]\)))/(k1 - k2));
  27. (*复斜率等于k1,过点A1与复斜率等于k2,过点A2的直线交点*)
  28. Jd[k1_, a1_, k2_, a2_] := -((k2 (a1 - k1
  29. \!\(\*OverscriptBox["a1", "_"]\)) - k1 (a2 - k2
  30. \!\(\*OverscriptBox["a2", "_"]\)))/(k1 - k2));
  31. FourPoint[a_, b_, c_, d_] := ((
  32. \!\(\*OverscriptBox["c", "_"]\) d - c
  33. \!\(\*OverscriptBox["d", "_"]\)) (a - b) - (
  34. \!\(\*OverscriptBox["a", "_"]\) b - a
  35. \!\(\*OverscriptBox["b", "_"]\)) (c - d))/((a - b) (
  36. \!\(\*OverscriptBox["c", "_"]\) -
  37. \!\(\*OverscriptBox["d", "_"]\)) - (
  38. \!\(\*OverscriptBox["a", "_"]\) -
  39. \!\(\*OverscriptBox["b", "_"]\)) (c - d));
  40. (*过两点A和B、C和D的交点*)

  42. \!\(\*OverscriptBox["FourPoint", "_"]\)[a_, b_, c_, d_] := -(((c
  43. \!\(\*OverscriptBox["d", "_"]\) -
  44. \!\(\*OverscriptBox["c", "_"]\) d) (
  45. \!\(\*OverscriptBox["a", "_"]\) -
  46. \!\(\*OverscriptBox["b", "_"]\)) - ( a
  47. \!\(\*OverscriptBox["b", "_"]\) -
  48. \!\(\*OverscriptBox["a", "_"]\) b) (
  49. \!\(\*OverscriptBox["c", "_"]\) -
  50. \!\(\*OverscriptBox["d", "_"]\)))/((a - b) (
  51. \!\(\*OverscriptBox["c", "_"]\) -
  52. \!\(\*OverscriptBox["d", "_"]\)) - (
  53. \!\(\*OverscriptBox["a", "_"]\) -
  54. \!\(\*OverscriptBox["b", "_"]\)) (c - d)));
  55. XiangjiaoxuanLianxin[o1_, a_, o2_, b_] := 1/(2 (
  56. \!\(\*OverscriptBox["o2", "_"]\) -
  57. \!\(\*OverscriptBox["o1", "_"]\))) (a
  58. \!\(\*OverscriptBox["a", "_"]\) - b
  59. \!\(\*OverscriptBox["b", "_"]\) +
  60. \!\(\*OverscriptBox["b", "_"]\) o2 + b
  61. \!\(\*OverscriptBox["o2", "_"]\) -
  62. \!\(\*OverscriptBox["a", "_"]\) o1 - a
  63. \!\(\*OverscriptBox["o1", "_"]\) +
  64. \!\(\*OverscriptBox["o2", "_"]\) o1 - o2
  65. \!\(\*OverscriptBox["o1", "_"]\));(*圆 （O1，A）与圆 （O2，B)连心线与公共弦的交点*)

  67. \!\(\*OverscriptBox["XiangjiaoxuanLianxin", "_"]\)[o1_, a_, o2_, b_] :=
  68.    1/(2 (o2 - o1)) (a
  69. \!\(\*OverscriptBox["a", "_"]\) - b
  70. \!\(\*OverscriptBox["b", "_"]\) +
  71. \!\(\*OverscriptBox["b", "_"]\) o2 + b
  72. \!\(\*OverscriptBox["o2", "_"]\) -
  73. \!\(\*OverscriptBox["a", "_"]\) o1 - a
  74. \!\(\*OverscriptBox["o1", "_"]\) + o2
  75. \!\(\*OverscriptBox["o1", "_"]\) -
  76. \!\(\*OverscriptBox["o2", "_"]\) o1);
  77. Chuizu[a_, b_, p_] := (
  78. \!\(\*OverscriptBox["a", "_"]\) b - a
  79. \!\(\*OverscriptBox["b", "_"]\) + p (
  80. \!\(\*OverscriptBox["a", "_"]\) -
  81. \!\(\*OverscriptBox["b", "_"]\)) +
  82. \!\(\*OverscriptBox["p", "_"]\) (a - b))/(2 (
  83. \!\(\*OverscriptBox["a", "_"]\) -
  84. \!\(\*OverscriptBox["b", "_"]\)));(*=(1/2)[p+(
  85. \!\(\*OverscriptBox["a", "_"]\)b-a
  86. \!\(\*OverscriptBox["b", "_"]\)+
  87. \!\(\*OverscriptBox["p", "_"]\)(a-b))/(
  88. \!\(\*OverscriptBox["a", "_"]\)-
  89. \!\(\*OverscriptBox["b", "_"]\))]P到直线AB的垂足*)

  91. \!\(\*OverscriptBox["Chuizu", "_"]\)[a_, b_, p_] := (a
  92. \!\(\*OverscriptBox["b", "_"]\) -
  93. \!\(\*OverscriptBox["a", "_"]\) b +
  94. \!\(\*OverscriptBox["p", "_"]\) (a - b) + p (
  95. \!\(\*OverscriptBox["a", "_"]\) -
  96. \!\(\*OverscriptBox["b", "_"]\)))/(2 (a - b));
  97. (*Duichendian[a_,b_,p_]:=(
  98. \!\(\*OverscriptBox["a", "_"]\)b-a
  99. \!\(\*OverscriptBox["b", "_"]\)+
  100. \!\(\*OverscriptBox["p", "_"]\)(a-b))/(
  101. \!\(\*OverscriptBox["a", "_"]\)-
  102. \!\(\*OverscriptBox["b", "_"]\));P关于直线AB的对称点

  104. \!\(\*OverscriptBox["Duichendian", "_"]\)[a_,b_,p_]:=(
  105. \!\(\*OverscriptBox["b", "_"]\)-
  106. \!\(\*OverscriptBox["a", "_"]\)b+
  107. \!\(\*OverscriptBox["p", "_"]\)(a-b))/(a-b);*)
  108. Wx[a_, b_, c_] := (a (-b + c)
  109. \!\(\*OverscriptBox["a", "_"]\) + b (a - c)
  110. \!\(\*OverscriptBox["b", "_"]\) + (-a + b) c
  111. \!\(\*OverscriptBox["c", "_"]\))/((-b + c)
  112. \!\(\*OverscriptBox["a", "_"]\) + (a - c)
  113. \!\(\*OverscriptBox["b", "_"]\) + (-a + b)
  114. \!\(\*OverscriptBox["c", "_"]\));
  115. \!\(\*OverscriptBox["Wx", "_"]\)[a_, b_, c_] := ((b - c)
  116. \!\(\*OverscriptBox["b", "_"]\)
  117. \!\(\*OverscriptBox["c", "_"]\) +
  118. \!\(\*OverscriptBox["a", "_"]\) ((a - b)
  119. \!\(\*OverscriptBox["b", "_"]\) + (-a + c)
  120. \!\(\*OverscriptBox["c", "_"]\)))/((-b + c)
  121. \!\(\*OverscriptBox["a", "_"]\) + (a - c)
  122. \!\(\*OverscriptBox["b", "_"]\) + (-a + b)
  123. \!\(\*OverscriptBox["c", "_"]\));

  125. e = Chuizu[i, b, a];
  126. \!\(\*OverscriptBox["e", "_"]\) =
  127. \!\(\*OverscriptBox["Chuizu", "_"]\)[i, b, a];
  128. j = Jd[k[e, f], e, -p n, d];
  129. \!\(\*OverscriptBox["j", "_"]\) =
  130. \!\(\*OverscriptBox["Jd", "_"]\)[k[e, f], e, -p n, d];
  131. o = Wx[d, e, f];
  132. \!\(\*OverscriptBox[
  133. StyleBox["o",
  134. FontSize->14], "_"]\) =
  135. \!\(\*OverscriptBox["Wx", "_"]\)[d, e, f]; k = p  n (
  136. \!\(\*OverscriptBox["d", "_"]\) -
  137. \!\(\*OverscriptBox["o", "_"]\)) + o;
  138. \!\(\*OverscriptBox["k", "_"]\) = (d - o)/(p n) +
  139. \!\(\*OverscriptBox["o", "_"]\); q = Wx[j, k, f];
  140. \!\(\*OverscriptBox[
  141. StyleBox["q",
  142. FontSize->24], "_"]\) =
  143. \!\(\*OverscriptBox["Wx", "_"]\)[j, k, f];(*(d-o)/p n*)
  144. t = XiangjiaoxuanLianxin[q, f, i, m];
  145. \!\(\*OverscriptBox["t", "_"]\) =
  146. \!\(\*OverscriptBox["XiangjiaoxuanLianxin", "_"]\)[q, f, i, m];
  147. Simplify[{d,
  148. \!\(\*OverscriptBox["d", "_"]\), e,
  149. \!\(\*OverscriptBox["e", "_"]\), f,
  150. \!\(\*OverscriptBox["f", "_"]\), 1, j,
  151. \!\(\*OverscriptBox["j", "_"]\), o,
  152. \!\(\*OverscriptBox["o", "_"]\), k,
  153. \!\(\*OverscriptBox["k", "_"]\), 2, q,
  154. \!\(\*OverscriptBox["q", "_"]\), 3, t,
  155. \!\(\*OverscriptBox["t", "_"]\), , t  
  156. \!\(\*OverscriptBox["t", "_"]\)}]
  157. Factor[{d,
  158. \!\(\*OverscriptBox["d", "_"]\), e,
  159. \!\(\*OverscriptBox["e", "_"]\), f,
  160. \!\(\*OverscriptBox["f", "_"]\), 1, j,
  161. \!\(\*OverscriptBox["j", "_"]\), o,
  162. \!\(\*OverscriptBox["o", "_"]\), k,
  163. \!\(\*OverscriptBox["k", "_"]\), 2, q,
  164. \!\(\*OverscriptBox["q", "_"]\), 3, t,
  165. \!\(\*OverscriptBox["t", "_"]\), , t  
  166. \!\(\*OverscriptBox["t", "_"]\)}]



复制代码

毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复

使用道具 举报

markfang2050
发表于 2021-11-24 12:35:14 | 显示全部楼层
有高版本了，还在搞这低版本有意义？

点评

markfang2050
注册机的问题那是你的水平问题，12.0我都破解了。  发表于 2021-11-24 23:06
dlsh
10.3注册机不行  发表于 2021-11-24 19:46
uk702
我试了高版本的，问题也一样。  发表于 2021-11-24 18:26
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 1 反对 0

使用道具 举报

TSC999
发表于 2021-11-24 19:03:27 | 显示全部楼层
问题出在程序最前面缺少了一条清除所有变量的指令。加上这条指令就没有这个毛病了。

点评

dlsh
谢谢老朋友  发表于 2021-12-2 19:58
TSC999
以前没有遇到过是因为，不是每个程序都必须在前面有这条指令。但是程序开头就写上这个指令是个好习惯。  发表于 2021-11-25 08:11
dlsh
谢谢，确实如此，但是以前没遇过  发表于 2021-11-24 20:07
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-23 16:23 , Processed in 0.026418 second(s), 21 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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