开放式计算软件OpenLu－工程计算助手

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
- [. d7 I" a4 ~0 T    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
$ y6 S! e/ g2 I% Y. {4 e
    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
. C0 P6 O# q5 c% H( L
    [例子] 设一阶微分方程组及初值为：
! F# f1 `4 F, ^: A& Z- H& P& f        r'=2r-2rf, r(0)=1
! b0 v: t$ F8 K6 D  h        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。
7 i' A/ L2 U- R% Y
- `" ?% m7 ~. n6 R$ D& d    程序如下：

2. !!!using["IMSL","math"];
   
3. , u# [! ^- i8 W4 |. |' F+ H# X
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
   
5. * E/ f& [: x, b/ b/ D\\" M
6. ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

!!!using["IMSL","math"];<br /> <span style=\\"display:none\\">, u# [! ^- i8 W4 |. |' F+ H# X</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<br /> <span style=\\"display:none\\">* E/ f& [: x, b/ b/ D\\" M</span> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   * s' U7 @( A9 g$ q. L
2. # p9 O8 t  E6 M) B5 @
   
3. 1.  7.73453e-002  1.46445/ s! ~! V* o6 E; o
   
4. 
   # @6 ?! g* T+ K4 a, }6 c
5. 2.  8.49774e-002  0.577954- G) v1 D( _: V3 K
   
6.   \\" }. W. T5 p8 B
   
7. 3.  0.290891      0.249253  s. t+ P8 x* r- d/ A
   
8. * N' n% a8 l6 W: o9 J; D
   
9. 4.  1.4466        0.1872199 k9 @0 [& K2 c  `& q
   
10. 
    3 p  f/ {0 b( C
11. 5.  4.05146       1.43948
      L4 L& n* K2 _1 H
12. 5 o5 }1 M/ t9 a, i
    
13. 6.  0.175618      2.2586! o& E' Z: Y9 q! _2 y3 z0 W7 }
    
14. . ~3 I/ a- S& J2 C/ ]
    
15. 7.  6.53112e-002  0.9088
    8 d1 u\" w  X8 A2 F5 n9 `
16. 8 A# \: M* n7 r7 W3 e
    
17. 8.  0.147227      0.366718
    - W& M9 |# B  `* {
18. 
    ( R2 D! E: a: O2 R/ p
19. 9.  0.650596      0.187575
    * r% z, d3 p; _/ A3 G
20. 
    ( V* @8 A7 H1 S+ E% P) D2 \$ N
21. 10. 3.14433       0.348821

复制代码

详细说明：http://www.forcal.net/sysm/lu1/openlu_readme.htm
. p* {" L7 X! `) R, z2 B    下载试用：http://www.forcal.net/xiazai/lu1/openlu32.rar

孤寂冷逍遥

liupeng723911

说的不错！