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

forcal

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

    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
$ t7 t. Y. s- n
! |( D8 g. M4 Z/ M- M    [例子] 设一阶微分方程组及初值为：
" u% R3 a7 T' }4 g6 e        r'=2r-2rf, r(0)=1
        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。

5 s/ r6 _; Z- z8 u0 E* _, v    程序如下：

2. !!!using["IMSL","math"];2 L' H0 P4 e& D3 G3 X8 F
   
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义6 b/ j% h5 Z) A* {( |  @3 `7 E- @( n8 y
   
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\\">2 L' H0 P4 e& D3 G3 X8 F</font><br /> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<font class=\\"jammer\\">6 b/ j% h5 Z) A* {( |  @3 `7 E- @( n8 y</font><br /> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.. E& m' X+ E7 X' K
   
2.   A* s- I% A, \+ k3 k5 }' B
   
3. 1.  7.73453e-002  1.46445# Y+ Y# y1 A  v4 V
   
4. 0 U- o# D( t$ a# j
   
5. 2.  8.49774e-002  0.577954
   + f( F' I: T7 z
6. 
   \" T* q0 ]# E( l8 l  Y
7. 3.  0.290891      0.249253
   . T* G8 y- N\" H! o9 k1 L
8. $ w0 w3 S' T) t  I( T
   
9. 4.  1.4466        0.187219+ [, P4 e/ F5 u- [5 S
   
10. 7 s+ X! ~2 k& p- h
    
11. 5.  4.05146       1.43948
    ; f8 |\" H9 r- F) C8 R
12. 
    ' A& C/ W0 x+ `' h% L\" J4 d
13. 6.  0.175618      2.2586
    1 f/ ~0 O7 ?+ R4 ~1 M
14. 
    , g3 J+ d( Y; @; o+ a8 J2 K* `
15. 7.  6.53112e-002  0.9088; Q/ l2 g/ L% r3 X5 G$ E
    
16. 
    $ G5 J2 J/ K* f0 q7 C1 P, {2 N
17. 8.  0.147227      0.366718
    7 c( G8 [9 r/ O0 {7 ?4 C' h
18. 
    \" O, I- M$ z. R) ], X
19. 9.  0.650596      0.187575
    ! U7 j# z5 S7 R9 c4 P# ^8 r
20. 4 U5 `: `) @) h/ X: z
    
21. 10. 3.14433       0.348821

复制代码

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