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

forcal

    OpenLu是OpenFC的升级版本，OpenLu用Lu作为支持脚本。OpenLu力求消除冗余，简化计算。基本计算内容包括表达式计算、解非线性方程（组）、多元积分、微分方程求解、参数优化拟合、矩阵运算等等。
    OpenLu工作模式有三种，其即时编译计算模式特别适合于常用计算：在代码窗口写代码时，即时给出代码中的错误，若没有错误，将进行计算给出结果。
% G4 k- G, k% g    OpenLu不仅支持整数、实数、复数、三维向量等基本运算，还可以自定义常量、函数、命令菜单等。按需定制，简化计算是OpenLu的基本特点。
9 y. ^- ?9 N7 w! v  N4 R  T% J
    本次升级，OpenLu简化了LuIMSL、LuOpt等扩展库，保留求方程（组）全部解、参数优化拟合等基本功能，函数用法也有诸多简化，以微分方程求解为例：
" e$ U  ?, A  U' e% v
. o: e& h* e. t: o4 `    [例子] 设一阶微分方程组及初值为：
8 _; x+ S( B; J) O5 g/ w8 r        r'=2r-2rf, r(0)=1
' i( b7 w/ W+ Q        f'=-f+rf,  f(0)=3
    计算t=1,2,...,10时的r、f的值。
. E- L9 i" D, T6 d1 T4 t
    程序如下：

2. !!!using["IMSL","math"];
   
3. 9 E' ]- s6 G& f! `
4. f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义
   
5. ) L' e7 |$ X3 z5 F
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\\">9 E' ]- s6 G& f! `</span> f(t,r,f,dr,df)={dr=2*r-2*r*f, df=-f+r*f}; //函数定义<br /> <span style=\\"display:none\\">) L' e7 |$ X3 z5 F</span> ode[@f,ra1[0,1,2,3,4,5,6,7,8,9,10],ra1[1,3]].outa[];

结果：

1. 0.  1.            3.
   - ~* i  a* W; s: ^\" L
2. 
   $ u9 t2 Q& J% I8 r: W9 J, o
3. 1.  7.73453e-002  1.464450 G+ v3 ~5 _0 C' I7 y
   
4. 
   9 B8 J\" L4 X. @& {  v8 W, d
5. 2.  8.49774e-002  0.577954\" E' p$ R3 ]) v# z
   
6. 
   / k) a* R1 M( D! y$ \- o, G
7. 3.  0.290891      0.249253
   ! Z) o6 m6 r3 `* |' E! h
8. 5 Q' B* r: Z' W/ ~7 w
   
9. 4.  1.4466        0.187219
   1 l. x1 h' k, k- `\" {2 z
10. 1 Y) U/ @# e& p8 _4 b& o
    
11. 5.  4.05146       1.43948* o6 ~: @# Z3 w' G0 l, `* E
    
12. 
    ) \, R1 N) N5 j( B& c6 B
13. 6.  0.175618      2.2586
    ; ^6 f- |8 y3 Z$ s
14. 
    % i$ n% n; X\" k( r+ a: T# K1 u8 F
15. 7.  6.53112e-002  0.9088) q2 P4 O/ v5 F
    
16. / U! [5 z, e9 [
    
17. 8.  0.147227      0.3667183 S/ m; t2 u4 B) D- J
    
18. 
    3 t. x# i/ |# Q( C! ~1 T& f
19. 9.  0.650596      0.187575
    : K: b# T6 e3 U\" X, a  u9 I
20. 
    ) ^# k& y: \' @; s7 d: d7 w* m
21. 10. 3.14433       0.348821

复制代码

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

liupeng723911

说的不错！

孤寂冷逍遥