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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
4 q2 |* }8 N5 `, {" L; @" Y+ x    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。

7 H2 G7 X# [  n5 s) J2 R# r2 M3 Z    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：

. Y7 ?+ x8 `/ D6 Z/ B, F    [例子] 设一阶微分方程组及初值为：
        r'=2r-2rf, r(0)=1
        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。

" b9 l% k8 u% T0 Z    程序如下：

2. !!!using["IMSL","math"];$ i; v- [5 w: w* u+ Y- J* P! c
   
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
   
5. 3 E5 A( l2 e6 r3 b0 g\\" M* D
6. ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

!!!using["IMSL","math"];<font class=\\"jammer\\">$ i; v- [5 w: w* u+ Y- J* P! c</font><br /> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<br /> <span style=\\"display:none\\">3 E5 A( l2 e6 r3 b0 g\\" M* D</span> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.! g  r  _8 @- Y  [- d- C
   
2. 
   . T& u4 N) |. p  t; X* O
3. 1.  7.73453e-002  1.46445
   ; Y: [8 y: w! }* W
4. 1 n/ l+ {# t! M* s; |/ o! d0 Z
   
5. 2.  8.49774e-002  0.5779549 g3 u5 Z9 a- _5 p. ^
   
6. 
   8 F3 B\" z7 C7 o- e6 ]! R
7. 3.  0.290891      0.249253) y5 _3 M, A) h5 H+ @0 n+ U) W\" h7 K
   
8. & z5 F3 J9 M* R0 N
   
9. 4.  1.4466        0.187219$ h$ l\" K1 S; h3 i' K. Z
   
10. 
    3 F9 W) w- e9 B: @! q
11. 5.  4.05146       1.43948
    9 \/ X  j4 J5 ]' g4 d, l
12. , @, z& q! u5 h. @
    
13. 6.  0.175618      2.2586
    $ h+ f7 A1 S+ f6 S\" s
14. 5 V+ w/ I9 k% F
    
15. 7.  6.53112e-002  0.9088
    * J$ R( s3 S# \8 c
16. ! b  v0 i. X8 j: ^) w9 r- U
    
17. 8.  0.147227      0.3667188 Q$ j( B: v* }( N- ?( n
    
18. : Y# r- e7 F& R9 q. @% L
    
19. 9.  0.650596      0.1875753 r' ^8 e3 ]2 I' I
    
20. , j; ?) _  _4 h\" m& P\" [
    
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！