this is a simple program for initial magnetization simulation 