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
