数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
楼主: 天山草

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

[复制链接]
发表于 2013-10-24 22:31 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

好,大家一起交流
发表于 2013-10-24 23:57 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

下面引用由天山草2013/10/24 05:28pm 发表的内容:
风花飘飘提出了一个【一元三次方程的求解公式】,根据此公式编程如下。
注:方程的一般形式是 x^3 + ax^2 + bx + c = 0。

A1BED2BBD4AAC8FD_1382630250.gif

134.89 KB, 下载次数: 1

 楼主| 发表于 2013-10-25 09:45 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

下面引用由denglongshan2013/10/24 10:31pm 发表的内容:
好,大家一起交流
denglongshan 学习 mathematica 比本人早,期待您提供更多更好的学习经验。
发表于 2013-10-27 22:31 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

Mathematica我主要是在证明几何结论时用于符号计算。有时仍然难以计算出结果,还期望楼主帮助。Mathematica经常用就熟悉了。压缩文件是一段证明概率考的一道很难的几何题的源程序,没有算出结论。

AbonisB7D6B5E3D5_1382884292.rar

63.74 KB, 下载次数: 1

发表于 2013-10-27 22:57 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

楼主老师:可以把你在 http://www.mathchina.com/cgi-bin/topic.cgi?forum=5&topic=18665 的源程序发过来吗?站内短信也可以。
 楼主| 发表于 2013-10-28 17:08 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

[这个贴子最后由天山草在 2013/10/28 05:11pm 第 1 次编辑]

直接发 mathematica 文件不知是否能行,将其转换成了 txt 文件:

D3C320mathematic_1382951297.txt

1.14 KB, 下载次数: 1

 楼主| 发表于 2013-10-28 17:21 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

[这个贴子最后由天山草在 2013/10/28 05:36pm 第 1 次编辑]

奇了怪了,今天运行这个文件,为何不是先前那个结果了?
显示了下述信息:
A very large output was generated. Here is a sample of it:《1》
现在明白了,原先是给 x 和 y 赋过具体的值,之后虽然消去了这些值,但是其作用还存在。如何解决这个问题?
-------------------------------------------------
BC = 1; a = 1; AB = Sqrt[
x^2 + y^2]; c = AB; AC = Sqrt[(1 - x)^2 + y^2]; b = AC; p = (
a + b + c)/2;
r = Sqrt[((p - a) (p - b) (p - c))/p]; JB =
ArcCos[(a^2 + c^2 - b^2)/(2 a*c)]; JC =
ArcCos[(a^2 + b^2 - c^2)/(2 a*b)];
XD = r/Tan[JB/2]; YD = 0;
XE = 1 - (1 - XD)*Cos[JC]; YE = (1 - XD)*Sin[JC];
XF = XD*Cos[JB]; YF = XD*Sin[JB];
k = YE/XE; k1 = YF/(1 - XF);
YQ = (k^2*r + k*XD)/(1 + k^2) - (
  k Sqrt[k^2*r^2 + 2 k*r*XD - k^2*XD^2])/(1 + k^2);
XQ = (k*r + XD - Sqrt[k^2*r^2 + 2 k*r*XD - k^2*XD^2])/(1 + k^2);
YR = k1 - k1^3/(1 + k1^2) + (r*k1^2)/(1 + k1^2) - (k1*XD)/(
  1 + k1^2) - (
  k1 Sqrt[-k1^2 + 2 k1*r + k1^2 r^2 + 2 k1^2*XD - 2 k1*r*XD -
    k1^2*XD^2])/(1 + k1^2);
XR = (k1^2 - k1*r + XD +
  Sqrt[-k1^2 + 2 k1*r + k1^2 r^2 + 2 k1^2 XD - 2 k1*r*XD -
   k1^2 XD^2])/(1 + k1^2);
YM = 0; XM = (XR*YE - XE*YR)/(YE - YR);
k2 = YE/(XE - XD); k3 = YF/(XF - XM); k4 = (YR - YQ)/(XR - XQ);
XG = (k2*XD - k3*XM)/(k2 - k3); YG = (k2*k3*XD - k2*k3*XM)/(k2 - k3);
XG1 = (k2*XD - k4*XQ + YQ)/(k2 - k4); YG1 = (k2*k4*XD - k2*k4*XQ +
   k2*YQ)/(k2 - k4);
If[XG == XG1 && YG == YG1, Print["G点与G1点重合"], Print["G点与G1点不重合"]]
 楼主| 发表于 2013-10-28 21:05 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

[这个贴子最后由天山草在 2013/11/08 09:22pm 第 5 次编辑]

1111_1383916958.gif

17.21 KB, 下载次数: 1

CFE0B5C8CFE0B5C8_1383105602.gif

9.85 KB, 下载次数: 1

0000_1383916642.gif

13.67 KB, 下载次数: 1

 楼主| 发表于 2013-10-28 21:10 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

[这个贴子最后由天山草在 2013/11/08 09:24pm 第 1 次编辑]


“mathematica 好像只能对于具体数字进行判断。” 这种说法是不正确的!

1111111111111111_1382965849.gif

25.12 KB, 下载次数: 1

发表于 2013-10-28 22:51 | 显示全部楼层

【分享】交流一下 Mathematica 编程方法,期望学习该软件的网友积极参与

[这个贴子最后由denglongshan在 2013/10/28 11:04pm 第 1 次编辑]
下面引用由天山草2013/10/28 05:21pm 发表的内容:
奇了怪了,今天运行这个文件,为何不是先前那个结果了?
显示了下述信息:
A very large output was generated. Here is a sample of it:《1》
现在明白了,原先是给 x 和 y 赋过具体的值,之后虽然消去了这些值 ...
只好关闭程序后又重新,你的源程序证明概率考的结论时,程序提示“产生了一个很大的输出”。
我的程序有问题吗?为什么计算不出结果?
    另外,张景中院士和合作者在研究几何自动证明器最近的一篇论文中,已经提到共轭比概念,这项研究是国家自然基金项目,参考:http://crad.ict.ac.cn/zhuantipaper/130918.pdf
,是在 http://bbs.cnool.net/cthread-104808052.html中找到的,但我看不懂。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2019-11-12 11:51 , Processed in 0.193187 second(s), 17 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

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