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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
2 d2 Q& P( b, h5 \  _! V    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
3 Y% Q" d3 M2 z0 ]0 \; G2 L  s6 c6 x
    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：

    [例子] 设一阶微分方程组及初值为：
7 m  d; {/ k* k/ P        r'=2r-2rf, r(0)=1
* j7 L$ O8 ?& ?+ e( e, w        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。
4 V3 i) ~) y4 ]7 ^6 D  w+ \
. A+ A) o, V' o8 h: _8 a1 X    程序如下：

2. !!!using["IMSL","math"];4 ~+ p7 B6 N8 I; o
   
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
   
5. : L0 V# I! _; r) b7 [2 [$ z
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\\">4 ~+ p7 B6 N8 I; o</font><br /> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<br /> <span style=\\"display:none\\">: L0 V# I! _; r) b7 [2 [$ z</span> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   7 z) A  L, V  q0 c
2. 0 e8 x& {1 k. Q) l5 M6 m1 {
   
3. 1.  7.73453e-002  1.464451 v- n! R) @' i. h6 K2 \
   
4. 
   / [+ ^! e2 O. |7 h/ R5 o
5. 2.  8.49774e-002  0.5779549 ?\" j0 G. H( F+ O; O$ l, Z
   
6. / f; c% u( n/ g; V
   
7. 3.  0.290891      0.249253* n# n7 X1 ~1 w, v
   
8. 
   % N: ~9 ^. }; {- @* s; L+ u6 [& |
9. 4.  1.4466        0.187219$ a& n: U5 `4 o6 D3 b: u
   
10. 
    0 F$ J; H* I: B1 R: f% j
11. 5.  4.05146       1.43948$ |& R6 V9 H& w) ]  ~# s, z
    
12. \" P* q+ Z7 M  g6 T* ]' |
    
13. 6.  0.175618      2.2586  b' ]% b! v) l) C3 s! {7 b
    
14. 4 y  [) d\" m\" A% k\" t
    
15. 7.  6.53112e-002  0.90888 Y+ `! Y( w! s7 m: i
    
16. 1 |+ c$ d4 d* K- _+ c9 O# D
    
17. 8.  0.147227      0.366718
    % z; v' p$ G5 a0 K
18. 
    0 u4 H1 V7 t1 K/ `5 B: H2 D  C8 ^
19. 9.  0.650596      0.187575
    $ C+ N/ I6 s- G( b4 {
20. 
    + d) `- `9 g/ Q* [7 M- l; S
21. 10. 3.14433       0.348821

复制代码

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

liupeng723911

说的不错！

孤寂冷逍遥