本文目录一览：
如何在lammps中用蒙特卡洛方法进行分子模拟？
LAMMPS是专门做分子动力学（MD）模拟的程序。它跟蒙卡完全是两种截然不同的方法。简单的蒙卡code规模很小，你可以自己写 —— 需要一个随机数发生器，一个抽样算法（比如metropolis，接近临界时可用swendson-wang），一个构型空间（如离散化的模型）。针对专门的问题，一般会有专门优化过的程序甚至做好的构型供抽样。
“ LAMMPS is a classical molecular dynamics code that models an ensemble of particles in a liquid, solid, or gaseous state. It can model atomic, polymeric, biological, metallic, granular, and coarse-grained systems using a variety of force fields and boundary conditions. ”
当然LAMMPS也提供了蒙卡功能：
巨正则系综蒙卡：
direct simulation Monte Carlo：
lammps该怎么入门
这个...感觉Lammps入门还是要有段时间的，不过我个人不是很擅长学新的东西。
命令的话不要去背，太多了。
首先看你题目做那个方向，做传热相关还是力学相关? 那种材料？材料的势能场参数很关键，要确保能有好的势能参数，输到Lammps里头才能跑出好结果。
然后要先从例子开始。如果你有师兄师姐以前做过，那就好办了，先拿人家的input file（也就是给Lammps输入的命令序列）过来看，里面会已经有个能跑得框架就好，之后一些常用命令，比如怎么设置仿真盒子大小，读入原子坐标，设置压强温度控制，跑多少步什么的，直接抄过来（最好到网站上搜这些命令的用法，然后自己对照着自己想法改，或者有书也不错，查查用法。）
至于建模比较麻烦，因为这款软件不像有限元集成了建模功能，你得自己编程，matlab，fortran，c语言都可以，想办法把你要仿真的一堆原子坐标都算出来，然后打出到纯文本文件（里面记录每个原子的序号，类别，坐标和速度），再和你的input file放到一起给Lammps算。之后跑完Lammps给的结果还需要些后处理，有可能还要编点小程序来处理。
另外你是要用自己电脑跑吗？还是有学校超级计算机让你用。反正最好是能先跑起来，哪怕把人家例子（或者网站上的例子）跑会了就成功一大半了。光看书资料自己没感觉的。
我看没人回答就随便写上几句，希望你能成功。
c语言高手急救：计算以下公式的值：s=1\1^2+1\2^2+1\3^2+……+1\n^2
这个显然不能用整数除法来做，如果用整数除法的话，大概结果总是1.
我也没有什么好的思路，想到一个笨办法，用两个数模拟分数加法，比如 1/2 + 1/3,你就模拟通分运算，== 3/6 + 2/6 = 5/6，分别记录当前分子和分母，然后一直这样计算下去，因为求分母的最小公倍数有现成算法，而且从一到n可以用循环来实现，所以实际编程的难度并不像表面上看起来这么繁琐复杂。最后求得分数后可以再用double型出发来求小数解，这样在除法以前是没有计算误差的。
另外你的式子给的是 \ 而实际除法应该为 /, C语言我没有试过两者差别，但是matlab里两者是不一样的，这里请多加注意。
我的方法比较一般，希望有好的方法的话楼主贴出来，大家共同学习。
