| 查看: 516 | 回复: 1 | ||
[求助]
请教个mathematica的问题
|
|
mathematica输出表达式后，能否定义规则让它化简？ 如以下代码： sol = Solve[{x^2 + y^2 + z^2 == 1, x + y + z == a , x*y*z == b}, {x, y, z}]; xx = x /. sol yy = y /. sol zz = z /. sol Print["Writing Fortran Code . . . : / "]; SetDirectory["F:\\tang\\mathtest"]; strm = OpenWrite["test.f90", FormatType -> FotranForm, PageWidth -> 70]; (* write subroutine of invisopar*) WriteString[strm, "subroutine test(x,y,z,a,b)\n"]; WriteString[strm, "implicit none\n"]; WriteString[strm, "real*8::x,y,z,a,b\n"]; nroot = Length[xx]; For[ii = 1, ii <= nroot, WriteString[strm, "x = " <> ToString[FortranForm[xx[[ii]]]] <> "\n"]; ii++]; For[ii = 1, ii <= nroot, WriteString[strm, "y = " <> ToString[FortranForm[yy[[ii]]]] <> "\n"]; ii++]; For[ii = 1, ii <= nroot, WriteString[strm, "z = " <> ToString[FortranForm[zz[[ii]]]] <> "\n"]; ii++]; WriteString[strm, "end subroutine\n"]; Close[strm]; Print["Finished Writing Fortran Code . . . : / "]; 输出到test.f90后，第一个解为： x = a/3. + (-6 + 2*a**2)/(3.*2**0.6666666666666666*(-36*a + 20*a**3 + Sqrt(4*(-6 + 2*a**2)**3 + (-36*a + 20*a**3 - 216*b)**2) - 216*b)**0.3333333333333333) - (-36*a + 20*a**3 + Sqrt(4*(-6 + 2*a**2)**3 + (-36*a + 20*a**3 - 216*b)**2) - 216*b)**0.3333333333333333/(6.*2**0.3333333333333333) 希望有以下输出： tmp0 = 216*b tmp1 = -36*a + 20*a**3 tmp2 = (tmp1 - tmp0 )**2 tmp3 = -6 + 2*a**2 tmp4 = 4*tmp3**3 tmp5 = Sqrt(tmp4 + tmp2 ) tmp6 = tmp1 + tmp5 - tmp0 tmp7 = tmp6**0.3333333333333333 tmp8 = 2**0.6666666666666666 x = a/3. + tmp3 /(3.*tmp8 *tmp7 ) - tmp7/(6.*tmp8) 请问如何实现？多谢 P.S: maple有关键字可以实现上面类似的功能，如以下代码： with(codegen, fortran) A := array(1 .. 2, 1 .. 2, symmetric); A[1, 1] := log(x); A[1, 2] := 1-log(x); A[2, 2] := 2-log(x); print(A); fortran(A, optimized, mode = double) |
回复此楼» 猜你喜欢
江苏省自然基金 什么时候出结果
已经有6人回复
-
面上有专家说实验设备不是我们单位的
已经有3人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有11人回复
-
闲聊
已经有3人回复
-
如何查询2026年获批项目名称？
已经有3人回复
-
两块石头
已经有10人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有19人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有15人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
初秋的晨风
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用Wofram Mathematica解方程的问题
已经有9人回复
- Mathematica编程高级指引---通读之后可以解决虫友80%的关于Mathematica的问题
已经有158人回复
- 急切求助！一个mathematica求解一个简单一元非线性方程的问题
已经有4人回复
- Mathematica 简单 四元一次方程组 ，求解问题~
已经有9人回复
- 请教个mathematica画图的问题
已经有5人回复
- 请教一个Mathematica积分的问题（重金求解）
已经有6人回复
- mathematica请教计算结果问题
已经有4人回复
- 【求助】请教一个关于Mathematica中计算小量的问题
已经有4人回复
- 【求助】请问MATHEMATICA输出数据格式问题
已经有5人回复
- 【求助】求教mathematica7的积分问题和根号下代数运算问题【已完结】
已经有3人回复
- 【求助】mathematica化简问题【已完结】
已经有15人回复
- 【求助】【求助】mathematica中的一种数值积分问题
已经有5人回复
- 【讨论】mathematica绘图问题
已经有7人回复
- 【求助】用mathematica解决问题【已完结】
已经有7人回复
- 【求助】关于Mathematica绘图问题【已完结】
已经有3人回复
- 【求助】Mathematica画图问题
已经有7人回复
- 【讨论】Mathematica分段函数绘图问题
已经有8人回复
- 【求助】请教mathematica中关于Plot函数的问题【已完成】
已经有5人回复
- 【求助】请教Mathematica程序问题
已经有9人回复
- 【求助】再次关于Mathematica画图颜色问题（急用）
已经有3人回复
- 【求助】关于Mathematica画图颜色问题（急用）
已经有7人回复
- 【求助】mathematica5.0运行问题，各位大师兄，帮忙啊
已经有3人回复
- 【求助】origin和mathematica的问题！
已经有6人回复
- 【求助】紧急求助有关MATHEMATICA的问题
已经有9人回复
2楼2013-04-17 10:26:52