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

forcal

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

# V+ b% S$ Y3 V0 n; s) c5 `  ~    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
3 }* u0 E0 }* d, G  X
    [例子] 设一阶微分方程组及初值为：
0 M( u* f9 c, M) V* m6 `9 M. ?& C        r'=2r-2rf, r(0)=1
. z0 C; k- N' K" o; k/ P        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。
# C7 Z) Q$ q0 v0 j( r4 c6 k0 Y
6 h2 D- {4 |' N    程序如下：

2. !!!using["IMSL","math"];, Y7 g! {. T  ?- a
   
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义, t) Z0 x& C! a7 r0 G) t
   
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\\">, Y7 g! {. T  ?- a</font><br /> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<font class=\\"jammer\\">, t) Z0 x& C! a7 r0 G) t</font><br /> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   8 S; T( X7 r8 \3 W- C0 i
2. 
   ) M2 R' u2 W- \% n! @
3. 1.  7.73453e-002  1.464457 {  A4 x4 h8 w$ E' b+ I
   
4. 
   / |1 A3 h6 [  c# }- a' s
5. 2.  8.49774e-002  0.5779540 g% b+ j0 m6 P
   
6. ' p$ ^' R9 g( ?, U& ^
   
7. 3.  0.290891      0.249253
   & R% r2 D5 D  `. L/ F  u4 k( J
8. 
   9 V, w. Z) G* b
9. 4.  1.4466        0.187219
   ) N: q! t- u4 L' e- R2 E. A/ q
10. 
      d. ?, [1 l7 ]0 o- A4 x- d$ l
11. 5.  4.05146       1.43948
    # m+ o7 e! C\" n
12. + d6 R2 ~8 H5 G4 ]
    
13. 6.  0.175618      2.2586; l% \* o) |. w* X. T
    
14. 
    0 c5 V/ e/ s7 o; F9 a6 x8 v
15. 7.  6.53112e-002  0.9088
    5 L) n* I/ n: S6 u' E2 {0 i1 o, r
16. , C4 y& G) W) c+ X8 s* l! |4 T
    
17. 8.  0.147227      0.366718* H. Y! y. l- \\" L+ r  _. j
    
18. 
    . U& A2 ]1 C% W
19. 9.  0.650596      0.187575. v% `7 O. G, u4 W! c7 g
    
20. 1 ?9 q- m  L2 a+ T
    
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！