자기 거울 내부의 입자의 운동을 설명하고자 합니다. 자기 거울 내부에 multi particle이 존재하는 경우에는 multi particle을 plasma fluid로 처리한다고 합니다. steady state인 background plasma가 깔려있는 상태에서 test particle만 충돌에 의해 궤적이 바뀌는 방식의 접근을 취하고자 하는데, 유체 방정식이 너무 어렵네요. 저는 test particle이 구속되는 시간을 구하고자 하는데, 혹시 이와 관련된 matlab 프로그램이나 포트란 코드를 가지고 계시면 공유 부탁드립니다.