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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
$ K: P" m2 C% U5 @, d/ u" y) S    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
: D+ _* Z8 q4 X; b5 W9 o- ~  X    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。

# ]6 M" D, t0 p1 h, T. @: j    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：

    [例子] 设一阶微分方程组及初值为：
* }& @7 K5 p- |. n8 J        r'=2r-2rf, r(0)=1
7 n! W+ T) y, I        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。
$ K9 B* _% V$ q" g
    程序如下：

2. !!!using["IMSL","math"];
   
3. : D0 O  F: p/ e3 L! [# w
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
   
5. # |. @' d& F5 w$ L$ @
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\\">: D0 O  F: p/ e3 L! [# w</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<br /> <span style=\\"display:none\\"># |. @' d& F5 w$ L$ @</span> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3., G7 Y1 S3 v' w3 ^
   
2. + U3 ^% ~5 K( q9 Y* s3 I3 T6 R
   
3. 1.  7.73453e-002  1.46445
     d6 _+ B, c+ t0 i' z& Z
4. 
   9 x$ w% D: E6 D0 |: L6 N' R
5. 2.  8.49774e-002  0.5779549 Z' |- n$ t6 o+ Z& ?/ ^
   
6. 
   : Z& Z. R\" X: p9 ~; G& r6 O( h2 v
7. 3.  0.290891      0.2492536 y% c/ n' q  g) E  o8 J
   
8. 7 }7 P4 v% A# C2 u
   
9. 4.  1.4466        0.1872196 B3 D0 p- \% ^4 y# l. O# K\" U; U\" I
   
10. 7 ?1 M8 v$ h: J% }, k$ S
    
11. 5.  4.05146       1.439486 k& h- u# p  |9 V
    
12. 8 L1 X; @' V! o- Y$ N. Z! d) H$ e6 I+ K, D
    
13. 6.  0.175618      2.2586: ]\" _3 y$ Q9 X  H  f
    
14. 
    - e' Z2 Q0 C# u0 x8 l$ v
15. 7.  6.53112e-002  0.9088, k9 k5 S/ C: ]3 U  E
    
16. 1 m8 B* J* l: v+ Y! K; a8 b3 u4 ?
    
17. 8.  0.147227      0.366718
    : M5 I! V1 v3 K0 \' ~
18. 
    $ t8 C% m; s6 z& P5 r\" j
19. 9.  0.650596      0.187575
    0 J\" k  n3 w7 ?# l2 H
20. 7 i1 R. I$ H. g) s. i0 D
    
21. 10. 3.14433       0.348821

复制代码

详细说明：http://www.forcal.net/sysm/lu1/openlu_readme.htm
% z7 \4 X: p6 z5 ?  U    下载试用：http://www.forcal.net/xiazai/lu1/openlu32.rar

孤寂冷逍遥

liupeng723911

说的不错！