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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
, T- R" I' X/ R6 {" G    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
) p% l! y- w- w, m" i3 C  L    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
6 v0 j3 I, Z+ [/ _8 H7 x9 r
    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
! x  b! i  y) E
    [例子] 设一阶微分方程组及初值为：
0 l' S$ V' Y+ v& l& s; S        r'=2r-2rf, r(0)=1
        f'=-f+rf,  f(0)=3
, y, w/ T  f% s& c5 x0 E: }    计算t=1,2,...,10时的r、f的值。
! [/ d/ v1 X0 N: e
0 ]4 L5 E9 \9 N9 P% A+ w    程序如下：

2. !!!using["IMSL","math"];/ h1 n* n+ N& N) `1 q
   
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义2 Y# D* s( _6 M+ v+ W# f$ d/ M
   
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\\">/ h1 n* n+ N& N) `1 q</font><br /> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<font class=\\"jammer\\">2 Y# D* s( _6 M+ v+ W# f$ d/ M</font><br /> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   / f  d0 P3 i* u8 V/ T) \+ y
2. ! [7 B1 d5 n  E\" A4 ?8 G\" i
   
3. 1.  7.73453e-002  1.46445
   \" K0 a2 n; s( a+ g
4. 
   2 \$ @\" {3 M% w- `( d
5. 2.  8.49774e-002  0.577954* F\" R* H! I\" x' _
   
6. 0 I# ]: D* \, w$ E3 z
   
7. 3.  0.290891      0.249253\" u) t0 u\" Y$ m' R& K
   
8. $ m* ~9 @& l; p. p9 i. h2 g\" {7 _
   
9. 4.  1.4466        0.187219
   ' P& }# f5 I, W) `; a8 E
10. 
    3 R/ l- `# f5 e$ N& s\" {2 b! |
11. 5.  4.05146       1.43948
    ' a4 x+ Y: z9 c9 ]- k5 \8 m4 s$ d
12. 1 p; j1 Y; f) k0 _! X\" T: r
    
13. 6.  0.175618      2.25869 ]: {! q+ a9 C: s; n
    
14. 
    0 c6 I1 F\" G9 c3 M& c2 r% l/ q+ E
15. 7.  6.53112e-002  0.9088- _/ I/ W0 @+ l7 B  }2 `% |
    
16. 
    6 W5 D4 j7 W5 i; T, h0 d. n
17. 8.  0.147227      0.366718) D& F- N\" E. I2 A9 N3 a
    
18. 
    * t+ J3 t0 H' y' w5 K6 i
19. 9.  0.650596      0.187575
    5 y# Q1 E6 s7 t4 C8 u- B\" O
20. 
    9 P, Q0 I: q/ g, F/ K
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！