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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
& H: h+ }! y. s# q0 `) t- |
    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：

9 q9 s- }1 `( \& ~% A" h8 z+ Z( k6 M    [例子] 设一阶微分方程组及初值为：
        r'=2r-2rf, r(0)=1
        f'=-f+rf,  f(0)=3
% m3 g7 ~/ X* \/ e1 |* i- [$ }    计算t=1,2,...,10时的r、f的值。
7 K' D) ]9 Q  ~
    程序如下：

2. !!!using["IMSL","math"];
   
3. 6 z# P3 A+ Q: p% {4 M
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
   
5. 0 i\\" _% w1 |, @0 s% Q0 {3 M
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\\">6 z# P3 A+ Q: p% {4 M</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<br /> <span style=\\"display:none\\">0 i\\" _% w1 |, @0 s% Q0 {3 M</span> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   * B+ R: P) T' p: x
2. % t; j' t' _3 _  l: S\" J: `
   
3. 1.  7.73453e-002  1.46445
   : r) Q* ~. b) u4 ]- q6 T
4. 
   4 I& P' V* j9 {+ t4 o
5. 2.  8.49774e-002  0.5779548 N1 q5 \! x$ i3 c4 W7 \1 z2 }
   
6. 
   + p( O2 R( k5 Z  O
7. 3.  0.290891      0.2492530 m3 {% M. R  l
   
8. 6 i9 f\" E\" n4 f3 R8 ]$ t
   
9. 4.  1.4466        0.187219
   , ?  `  t( l. ]; v
10. 
    : ?& k3 q1 k\" W4 k! M
11. 5.  4.05146       1.43948
    $ x! m% |0 B- L' [' Y6 w
12. \" ?+ G- C4 o( P. F
    
13. 6.  0.175618      2.2586
    8 D* N8 H1 x' P3 @
14. . ]( b\" c0 |6 o9 K5 ]
    
15. 7.  6.53112e-002  0.90888 f# A( v& n3 M
    
16. - K3 v8 X: Y' a6 H1 c5 r; {8 ]( G/ X
    
17. 8.  0.147227      0.366718
    / |8 \1 L& `: W' h! p6 j
18. 7 q4 _+ W) v7 t% p9 E! E  d; @/ q; e
    
19. 9.  0.650596      0.187575
    + J3 R% [4 |) a
20. & Z# _0 m9 ]; E' w+ D9 S
    
21. 10. 3.14433       0.348821

复制代码

详细说明：http://www.forcal.net/sysm/lu1/openlu_readme.htm
1 U' g( s7 O9 L! [5 Q    下载试用：http://www.forcal.net/xiazai/lu1/openlu32.rar

孤寂冷逍遥

liupeng723911

说的不错！