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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
) D. J0 @) m6 U! _3 t# {    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
) a2 b. H& W; m. y% L3 c2 u    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
5 ?4 Y  V$ _9 M
) @* y0 o# l2 {/ B    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
, x  X5 Z" z9 G5 Q  ^, _& q
& m- n6 r# O  F3 S3 {; N    [例子] 设一阶微分方程组及初值为：
        r'=2r-2rf, r(0)=1
* g# _; t3 U) |, c* }+ ?        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。
, k5 W9 y! p  g0 V6 [+ h( J& i
, t+ R$ @  M. l8 R    程序如下：

2. !!!using["IMSL","math"];9 A( _/ G7 U% z' H: r# \1 D
   
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
   
5. 2 {; C$ B' U! g7 f. M
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\\">9 A( _/ G7 U% z' H: r# \1 D</font><br /> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<br /> <span style=\\"display:none\\">2 {; C$ B' U! g7 f. M</span> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   ! a- E) n4 Y: a  o
2. 
   9 X: ~0 S\" h9 }) e
3. 1.  7.73453e-002  1.46445
   - _% V; z) }& x3 Y, U6 a3 }/ u& f' y
4. ( o; k) D- F% C3 l\" Q
   
5. 2.  8.49774e-002  0.5779549 v  Y+ \\" N3 ~8 R& K$ H
   
6. 
   & A' S\" ~4 n% J
7. 3.  0.290891      0.249253: k$ m, B' O. w8 {
   
8. # X# D  a( v5 g  C7 p$ p
   
9. 4.  1.4466        0.1872190 f7 g! F! j7 f# a! n
   
10. / y9 n; @) ^5 Y! ^$ {* L! G  }' X
    
11. 5.  4.05146       1.43948
    ( K6 ?7 n$ G& \. b
12. , n% f* d0 M\" w6 D0 n
    
13. 6.  0.175618      2.2586
    6 Y, r. i; Q7 y  w( j/ {
14. 
    ; G! w. Z; A* j! T* V
15. 7.  6.53112e-002  0.90886 \! \+ d\" }. y# N6 f\" ^7 q
    
16. & s* @) A9 |5 b) I: j3 ~
    
17. 8.  0.147227      0.366718
    1 M, v7 K8 l/ M
18. 
    5 y7 \+ D  J( n* k4 O9 J
19. 9.  0.650596      0.187575/ C& e, c% ]* k/ z6 c+ J  ]3 A\" M
    
20. 
    \" l8 D& P. c7 N. s4 V. O' J/ K
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！