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

forcal

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

    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
, L' |/ t( Z$ Z0 Z, U7 K
    [例子] 设一阶微分方程组及初值为：
        r'=2r-2rf, r(0)=1
$ u# @" J( r, B2 m/ N: s/ r) A        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。
" r% [! \" }& b# l! x6 U- a. [
2 b; a7 _% J& s" |0 Z& L+ p9 v    程序如下：

2. !!!using["IMSL","math"];
   
3. 2 x1 P# r* S8 ]
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义! E8 Z* m$ Y6 ^0 u- h1 S. V* i- R
   
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\\">2 x1 P# r* S8 ]</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<font class=\\"jammer\\">! E8 Z* m$ Y6 ^0 u- h1 S. V* i- R</font><br /> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3./ V0 P# R* P/ m\" H5 f$ l
   
2. - T1 i) C+ z) W7 o$ `7 N0 @% r
   
3. 1.  7.73453e-002  1.46445
   - s1 r# |2 x8 A, I: H
4. 6 n2 d4 E/ z$ X; b: {
   
5. 2.  8.49774e-002  0.5779546 b# T8 m3 a/ F$ _# ~# l4 K  c
   
6. 
   \" x+ s$ }# G/ Q+ `\" `8 `
7. 3.  0.290891      0.249253\" u1 }3 y% @0 D& v; J
   
8. 2 ^2 e* N5 v, P& m
   
9. 4.  1.4466        0.187219, b6 k3 A' ?\" {0 P3 s: E1 G
   
10. 6 E\" V+ U\" f- a/ v: c4 r
    
11. 5.  4.05146       1.439489 d) w1 r/ j; J9 L/ H: [5 Y
    
12. ; \' i1 @- Z( l% \( `
    
13. 6.  0.175618      2.2586
    5 i: h0 c( b1 K  l- V8 X
14. ) O* j0 }. f/ Q# m! T* C
    
15. 7.  6.53112e-002  0.9088/ J1 U& H% @. ~$ v1 }  O7 Y6 m9 E
    
16. # i0 W8 ~1 T4 u) O1 y% k0 W- E9 P- ^3 {
    
17. 8.  0.147227      0.3667184 ?) y& t, \7 `\" `) E
    
18. 0 w\" w9 s) }  |
    
19. 9.  0.650596      0.187575
    - ^. o# u4 K( M
20. 
    4 Z9 `/ n* ^# {7 C' g& g
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！