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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
) R7 w# m0 X: Q! n: C7 n    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
; R  d6 d2 M  I0 F2 q- p; ]7 w) n
! Z- [9 c- p3 X+ w$ J; x! O7 M+ o, y/ s    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
) [4 G0 v+ z& L0 X
5 d' h3 p; Y) M    [例子] 设一阶微分方程组及初值为：
        r'=2r-2rf, r(0)=1
0 P1 ?5 ], W0 I- d8 t: ?. q! J        f'=-f+rf,  f(0)=3
7 H+ @3 Y# |3 U0 g6 J( C- [    计算t=1,2,...,10时的r、f的值。

3 b1 Y5 S, @. I5 k4 @    程序如下：

2. !!!using["IMSL","math"];
   
3. $ }0 g  U2 C+ v7 a7 @
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义2 ~: h+ G5 N/ j- `
   
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\\">$ }0 g  U2 C+ v7 a7 @</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<font class=\\"jammer\\">2 ~: h+ G5 N/ j- `</font><br /> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   ) L' Y7 N- S: K* e  O
2. ( Y\" X\" m6 o: J- ^0 C
   
3. 1.  7.73453e-002  1.464456 f% x6 m\" f5 P/ f, H8 N, \
   
4. 
   ) D8 E, n0 ~2 W/ B# U6 [+ _- o
5. 2.  8.49774e-002  0.5779549 A3 z# a5 e\" }; h, W: u. `6 g( p
   
6. 
   . p) |+ d& L$ w+ S( {9 A\" U
7. 3.  0.290891      0.249253& N/ C$ ]  t: `6 p4 Q+ _! j0 L3 Z
   
8. 
   % v2 ^- O2 ~# E5 `& Z
9. 4.  1.4466        0.187219
   / m0 P$ b6 V. c0 r( R
10. - o! \9 k8 C& P
    
11. 5.  4.05146       1.43948
    6 \  N4 z6 k8 J) j) P* P
12. 
    $ [# k; Q/ d, j0 }/ k; {
13. 6.  0.175618      2.2586* _& |9 D* y0 C* Q
    
14. , a* ^& p9 }# O& C9 |
    
15. 7.  6.53112e-002  0.9088$ W0 D1 H/ ?  k+ q- g3 P1 F2 P' @
    
16. 1 h  m. C0 [$ l/ S, _+ o
    
17. 8.  0.147227      0.3667182 f# [3 h  O- o( q
    
18. 
    0 J. V8 ?8 ^- Z) o  S! U
19. 9.  0.650596      0.187575% F  H- {9 e0 [4 m0 l  c
    
20. 
    # U2 H' O% d\" P3 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

说的不错！