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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
$ m1 Z- j" y' G5 ?$ H1 S, O7 M. X    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
/ ]1 Z2 q8 }. @- K3 H8 G
' v# H7 P4 ^# @( {" ?9 w: \8 A    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：

    [例子] 设一阶微分方程组及初值为：
* o: c7 J9 K; w( E; `        r'=2r-2rf, r(0)=1
        f'=-f+rf,  f(0)=3
; w, K+ _7 @# F* d# [    计算t=1,2,...,10时的r、f的值。

    程序如下：

2. !!!using["IMSL","math"];
   
3. . v( X# F  j# A9 Z
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
   
5. ! `! F4 }, R% N- i9 _/ `
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\\">. v( X# F  j# A9 Z</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<br /> <span style=\\"display:none\\">! `! F4 }, R% N- i9 _/ `</span> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   : |+ _: d1 }2 n* Q, B: F7 u9 ?\" P$ L: n, E
2. : s9 w. f0 P3 l* c* j7 x\" k
   
3. 1.  7.73453e-002  1.464452 \/ j! E) o# J. q, t/ a
   
4. : T4 U- g; [9 o$ ~6 T- l0 P
   
5. 2.  8.49774e-002  0.5779542 _+ h7 g7 I( \5 M
   
6. 
   % C# \+ r0 d* l
7. 3.  0.290891      0.249253
   3 y* w5 k! B* m8 ]
8. 5 m, }' ]; b. J6 x6 i\" s$ M
   
9. 4.  1.4466        0.187219
   $ Q6 ~# w9 y! l) W3 D
10. 
    ' r' H! }# L8 H' ]3 q! Y  |
11. 5.  4.05146       1.43948
    + E% C7 z/ ^; q- s
12. * y4 ~# h/ f* r3 K0 R3 X( {8 N- i6 Z
    
13. 6.  0.175618      2.25867 \' e- i4 l\" P3 {- z: |
    
14. 
    - }, G. [5 J6 d
15. 7.  6.53112e-002  0.9088% M& u5 K7 E/ J5 e; ~
    
16. 1 |! S' F: P9 [
    
17. 8.  0.147227      0.366718
    5 {2 d! A9 @* k1 ?/ O' e0 R7 C/ x4 c
18. 
      k0 P( Y+ h. n\" ^9 H0 C
19. 9.  0.650596      0.187575
    \" C: Y9 H1 C5 ~) ~
20. 
    0 n. x+ O! X1 {3 O% p; u
21. 10. 3.14433       0.348821

复制代码

详细说明：http://www.forcal.net/sysm/lu1/openlu_readme.htm
    下载试用：http://www.forcal.net/xiazai/lu1/openlu32.rar

liupeng723911

说的不错！

孤寂冷逍遥