RK45andRK54

gmres 方法。视图求解线性方程组A*x=b的解x。nXn的稀疏矩阵A必须是方程且应是大型稀疏矩阵。列向量b的长度必须为n。(gmres method. View solving linear equations A* x = b the solution x. nXn sparse matrix A must be equation and should be large sparse matrix. The length of the column vector b must be to n.)

用刷选法求出100之内的素数。并输出数有多少个。(100 prime number is calculated by the method of selection. And the output count how many.)

实现矿井采空区的孔隙率和粘性阻力系数，主要用于采空区瓦斯抽采(The porosity and viscosity resistance coefficient of mine goaf can be achieved.)

一个简单的计算块的运动矢量源码，使用四步搜索算法(A simple calculation of the block motion vector source, use the four-step search algorithm)

FFT with fix point 2*N

关于焊接残余应力的Abaqus子程序-DFLUX(The Abaqus subprogram-DFLUX on welding residual stress)

一维等离子体的plrc—fdtd算法程序(one-dimensional plasma plrc- fdtd Algorithm)

2维TE波的FDTD的源程序。周期边界条件，有限宽度的吸收边界条件(Two-dimensional TE wave FDTD source. Periodic boundary conditions, the limited width of the absorbing boundary conditions)

使用龙格-库塔间断Galerkin方法（RKDG）求解无粘流动Euler方程的程序，适用于计算流体力学RKDG方法的学习和研究。(Using the discontinuous Galerkin method (DG) for solving inviscid flow Euler equations procedures applicable to DG methods of study and research.)

幻方 (magic cube )