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

forcal

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

0 G* ~+ ]! b1 P$ o3 T    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：

" q" u3 Q) d6 M' s( B: e$ R    [例子] 设一阶微分方程组及初值为：
        r'=2r-2rf, r(0)=1
        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。

' c  ?/ Y! |9 m- t% o    程序如下：

2. !!!using["IMSL","math"];
   
3. 0 d1 @# R, o) ~% ~- I
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义! e) P; w, [8 o% o0 A* `3 U, Y
   
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\\">0 d1 @# R, o) ~% ~- I</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<font class=\\"jammer\\">! e) P; w, [8 o% o0 A* `3 U, Y</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 _+ ?7 e3 I8 `
   
2. 2 u+ D/ ^! e\" x4 X, H# H
   
3. 1.  7.73453e-002  1.46445# [* F4 d' {! Q, f9 P, x
   
4. 
   \" G6 o7 w& l. G2 N' |4 }4 `4 U
5. 2.  8.49774e-002  0.5779547 R# c' r2 i5 U3 E/ i
   
6. 2 T) L  q- ^6 |0 h1 V5 \) S
   
7. 3.  0.290891      0.249253. p3 N! z/ i8 D2 h
   
8. 
   ( W$ k% T$ {/ t7 p' a2 N& ^8 F
9. 4.  1.4466        0.1872194 `; a. m: U, J% l
   
10. 6 }4 P9 M1 d5 {* m6 D! p4 @9 o
    
11. 5.  4.05146       1.439482 h3 U; U! r, |4 M( d\" G6 G& C. ?
    
12. 
    2 ^2 ?7 b) `3 G* {' x2 t) y
13. 6.  0.175618      2.25867 L) R; L+ J, g\" ~* H9 b* M
    
14. 7 p  a! ?$ S5 d) `, d: K: V# A
    
15. 7.  6.53112e-002  0.90881 }; @) l6 d: |( y# r1 s0 @
    
16. , X1 J$ K/ W\" D' X3 V4 P1 ^+ U
    
17. 8.  0.147227      0.366718
      \3 `8 S3 u/ m\" j) \& c
18. 
    ) M3 V4 T: Z0 {6 |2 ~2 J
19. 9.  0.650596      0.187575
    4 R* |0 E; {. X0 h( x
20. 
    4 }; v8 F% O! e+ v  Z7 L3 Y
21. 10. 3.14433       0.348821

复制代码

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

孤寂冷逍遥

liupeng723911

说的不错！