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

forcal

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

3 F- L# T, t6 V1 n* S" T; M) }    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：

( l6 i) Y  \- _& M    [例子] 设一阶微分方程组及初值为：
: g. g- h, K, c1 }- h; F9 m        r'=2r-2rf, r(0)=1
) r+ _3 u" Z1 Y* F- j# \        f'=-f+rf,  f(0)=3
8 N" T+ t" c- D8 m: Q  `    计算t=1,2,...,10时的r、f的值。
3 J7 h9 X/ [4 c/ i0 R
    程序如下：

2. !!!using["IMSL","math"];
   
3. 3 @\\" u4 |9 J- a5 ]9 J  b& Z1 P
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
   
5. 7 Z; i8 |+ X5 a: K
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\\">3 @\\" u4 |9 J- a5 ]9 J  b& Z1 P</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<br /> <span style=\\"display:none\\">7 Z; i8 |+ X5 a: K</span> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.& R) v: ^$ o; B  \# K
   
2. 0 S4 a6 c8 ?7 [+ b
   
3. 1.  7.73453e-002  1.46445
   9 l\" ~\" x5 y% D: {- ~
4. 
   # S- A; u7 q) i% f  m! s
5. 2.  8.49774e-002  0.5779543 O3 `+ q& G, B, g2 t% R' V
   
6. 
   8 [3 J' _' j  {/ E, f
7. 3.  0.290891      0.249253
   # \- z. T* I! v2 ^8 I
8. 
   2 I3 h& j( O. [$ p, D: P. i  M
9. 4.  1.4466        0.187219
   6 ^/ S' t' Q; i0 S$ B- L4 Y! L6 F
10. 
    . J5 f9 y6 C- E2 V; i% q7 u2 R: k& q
11. 5.  4.05146       1.43948
    : \) a5 p& k1 R0 P% i7 q5 f7 `0 v
12. 
    8 v% ?; S; \. ^0 p' Q' {
13. 6.  0.175618      2.2586
    2 I) g( ~- J/ [* e4 a% [% `
14. 
    + g7 A( q( K+ W+ l$ S
15. 7.  6.53112e-002  0.9088
    6 k. ^4 O2 P% F5 t* c1 F# w$ x; i
16. 
    1 i& D; A1 c0 M. y
17. 8.  0.147227      0.366718
    4 R4 t! |- P, i9 K* R# \
18. 
    + C* T6 O0 ^2 X5 B3 Z) w. k$ j
19. 9.  0.650596      0.187575
    8 b0 \( \2 {; H\" f6 T/ P\" ?& }
20. / F. k- K# |+ W6 V0 n8 V( \3 `
    
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！