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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
7 B2 `1 n: @6 c/ h/ i+ {    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
4 F& z% O4 f  I5 u; [+ P
" ]: O5 g9 t# X- {, S8 q    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
: J6 l1 {+ Y2 v. X5 w
    [例子] 设一阶微分方程组及初值为：
% j% J# q0 r$ A3 }$ n        r'=2r-2rf, r(0)=1
( M$ z: Q# r) U8 a        f'=-f+rf,  f(0)=3
1 z* C6 D- O% d# k. k9 `% Y4 R/ [6 m    计算t=1,2,...,10时的r、f的值。

    程序如下：

2. !!!using["IMSL","math"];% j- b! R6 E- y4 d
   
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义* Q# E6 h5 t! O4 U+ Q6 h$ D9 x
   
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\\">% j- b! R6 E- y4 d</font><br /> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<font class=\\"jammer\\">* Q# E6 h5 t! O4 U+ Q6 h$ D9 x</font><br /> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   ( j* l; g& T6 b, u/ l0 V
2. 
   + J& E7 c3 X( N2 e: q# H
3. 1.  7.73453e-002  1.46445
   3 i! R0 ~1 ~/ w0 L
4. 4 {4 f# r+ @5 J. R# C+ O
   
5. 2.  8.49774e-002  0.577954- K\" l/ V/ B% @/ W! Y3 ]
   
6. % I/ x, u2 _/ _6 A$ p
   
7. 3.  0.290891      0.249253
     B5 q4 z1 i6 f* L! Q6 }
8. . \2 e5 `* ^1 _! Y4 `
   
9. 4.  1.4466        0.187219% w3 B# t6 r! i; \  V' V
   
10. 
    ! L% c$ T) V# S7 o2 O
11. 5.  4.05146       1.439483 P9 y+ e, }/ n9 p
    
12. 
    \" l( p3 i- Y% S) _9 u
13. 6.  0.175618      2.2586, a) C. x5 g5 Z! x
    
14. 3 v# T+ p8 e# a' b
    
15. 7.  6.53112e-002  0.9088
    9 l, l7 R2 O* I; T& B
16. 3 c8 J7 I\" k  s0 w- a+ T, U$ s
    
17. 8.  0.147227      0.366718
      \* X2 ]/ m3 R
18. 9 b; G. i- a/ g0 L2 Y, J( T0 C3 g
    
19. 9.  0.650596      0.187575
    - Y! G5 i* l, Y: @\" J; q+ h
20. 
    . d\" S2 c& b+ }, T, E
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！