yp) integer,parameter :: L=10 integer,parameter :: N=20 real,parameter :: e=0.1 integer a(L,L) real x(N),y(N),xp(N),yp(N) real :: dr=0.1 real p,q,rnd1,rnd2 integer i,j integer :: k=1 do i=1,L do j=1,L a(i,j)=0 end do end do do while(k10) then y(k)=y(k)-10.0 end if a(ceiling(p*L),ceiling(q*L))=1 xp(k)=x(k) yp(k)=y(k) k=k+1 end if end do return end !在y方向有小的扰动 !标记已用过的格点 计算力的子程序 subroutine Force(x,y,dx,dy,Fx,Fy) implicit none integer,parameter :: N=20 integer,parameter…
Words 602 - Pages 3