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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
9 E+ P( L$ t. ^$ w% D: F    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
1 N( [# Z4 d/ s4 G/ d
/ Q+ D' f- F/ X9 X: I% c4 s! |+ f    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
8 a" Z) |& r% Z- G; \
    [例子] 设一阶微分方程组及初值为：
% z  N1 |# A/ O* l3 ^" `! e2 B        r'=2r-2rf, r(0)=1
0 e5 f4 M& h6 V0 f        f'=-f+rf,  f(0)=3
0 Q, w7 ]/ U' Q5 X, f# A3 r" M+ T    计算t=1,2,...,10时的r、f的值。

    程序如下：

2. !!!using["IMSL","math"];
   
3. : g# J, i\\" K8 ]( S0 ^5 d& _% ?
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义( C\\" d( g6 X6 E0 R
   
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\\">: g# J, i\\" K8 ]( S0 ^5 d& _% ?</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<font class=\\"jammer\\">( C\\" d( g6 X6 E0 R</font><br /> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   1 c4 ?1 k0 T. Y; c5 u! t
2. 
   . L, C5 v: _' h& ^# x0 [' o
3. 1.  7.73453e-002  1.46445
   / _/ g# x3 ~# i2 o7 W2 h
4.   d0 f1 S* L* p2 U2 v\" D+ I. \
   
5. 2.  8.49774e-002  0.577954
   # \* w0 l; g3 |8 J8 E- d
6. 
   2 J; \. V& a% ^( {7 J# \2 F+ E
7. 3.  0.290891      0.249253
   - t; O2 w( A\" G% }9 D7 y\" Q
8. , ^# _( P& m4 L5 D9 N' u
   
9. 4.  1.4466        0.187219
   0 ~$ i% c6 U4 _- Y
10. 
    9 m' a2 K* u2 }2 F5 X
11. 5.  4.05146       1.439486 \5 V8 ~8 X: G\" I6 y, F
    
12. : |- }! \4 }! Y, O- C
    
13. 6.  0.175618      2.2586  @' K$ X\" e6 r: J- ^! O
    
14. 
    2 Q\" i. r+ r& B% a7 ?\" X
15. 7.  6.53112e-002  0.9088) K' J: H3 z; U: F( H) [+ a% H
    
16. 
    1 H) l0 ~. B1 n7 i- m
17. 8.  0.147227      0.366718
    % x3 U9 r\" Y9 q; ?* N! I; |
18. # e& ^! B9 V' Q) G/ u9 [4 M
    
19. 9.  0.650596      0.187575
    \" B$ g: k6 d9 N# z; P0 B
20. ! r' t# V4 [% N, S: d
    
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！