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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
) V. z2 H: l. t% V0 P' q, C    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。

1 E' y+ W. V' [$ T    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
  J  N% x" G! R6 I: [0 m" X' ?/ d
    [例子] 设一阶微分方程组及初值为：
        r'=2r-2rf, r(0)=1
# R) Y3 Y) o2 D& e' R5 y( P& ~9 h3 h        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。

    程序如下：

2. !!!using["IMSL","math"];
   
3. ; o+ V8 N/ q) i6 U. p) s! d\\" K
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义9 k* ~6 q' w. a
   
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\\">; o+ V8 N/ q) i6 U. p) s! d\\" K</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<font class=\\"jammer\\">9 k* ~6 q' w. a</font><br /> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.$ y/ N, _& m% d% [# Q
   
2. 6 ~) H0 S0 L( u8 @
   
3. 1.  7.73453e-002  1.46445: l5 g8 k+ n9 C+ |9 x1 [8 y
   
4. 
   6 K) z2 f: T7 L% O0 z* S- D
5. 2.  8.49774e-002  0.577954
   * Q! ~2 |( u0 x5 P
6. 3 d! F) x4 K. {8 d1 i7 g
   
7. 3.  0.290891      0.249253/ U\" P9 X) @# z) a1 I/ K8 v& o$ T
   
8. 
   7 v1 m0 j# k3 u: p
9. 4.  1.4466        0.187219
   , z# Q7 Z& _9 [9 X
10. : x1 a6 F2 D% W8 O- Z# j& W
    
11. 5.  4.05146       1.43948  _0 ~$ T- L\" E7 e6 Q
    
12. * O5 v6 Z( B; Q9 E\" O$ W
    
13. 6.  0.175618      2.2586
    , m1 N7 x7 {\" O/ O. u; T# e8 Y7 k
14. ; ?  z# o/ s( |
    
15. 7.  6.53112e-002  0.9088
    ' \9 S% z9 @2 q* t$ b
16. 
    - s) U* H  @$ N
17. 8.  0.147227      0.366718
    % n/ j0 k3 C; P0 ?/ F! l; d) U
18. 8 b7 C. N8 E+ M. Z$ r
    
19. 9.  0.650596      0.187575
    # G( D! ?' ^9 A/ g7 n) `( w+ v
20. 
    0 X. H8 V' A; _$ C% @
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！