方向相的同步图和程序---程兴超.doc

（一）前三幅图的程序program coupleimplicit nonereal,parameter:e1=0.6,u=3.7integer,parameter:con=1real,dimension(100,100):x2,x3integer,dimension(100,100):x1integer i,j,t,time,k k=time() open(10,file=output10.txt)do j=1,100 do i=1,100 x2(i,j)=ran(k) end do end do do t=1,20000 do j=1,100 do i=1,100 1 if(i=1.and.j=1)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(i+1,j)*(1-x2(i+1,j)+x2(100,j)*(1-x2(100,j)+x2(i,j+1)*(1-x2(i,j+1)+x2(i,100)*(1-x2(i,100)-4*x2(i,j)*(1-x2(i,j)/4 2 elseif(i=1.and.j=100)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(i+1,j)*(1-x2(i+1,j)+x2(100,j)*(1-x2(100,j)+x2(i,1)*(1-x2(i,1)+x2(i,j-1)*(1-x2(i,j-1)-4*x2(i,j)*(1-x2(i,j)/4 3 elseif(i=100.and.j=1)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(1,j)*(1-x2(1,j)+x2(i-1,j)*(1-x2(i-1,j)+x2(i,j+1)*(1-x2(i,j+1)+x2(i,100)*(1-x2(i,100)-4*x2(i,j)*(1-x2(i,j)/4 4 elseif(i=100.and.j=100)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(1,j)*(1-x2(1,j)+x2(i-1,j)*(1-x2(i-1,j)+x2(i,1)*(1-x2(i,1)+x2(i,j-1)*(1-x2(i,j-1)-4*x2(i,j)*(1-x2(i,j)/4 5 elseif(i=1.and.j/=1.and.j/=100)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(i+1,j)*(1-x2(i+1,j)+x2(100,j)*(1-x2(100,j)+x2(i,j+1)*(1-x2(i,j+1)+x2(i,j-1)*(1-x2(i,j-1)-4*x2(i,j)*(1-x2(i,j)/4 6 elseif(i=100.and.j/=1.and.j/=100)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(1,j)*(1-x2(1,j)+x2(i-1,j)*(1-x2(i-1,j)+x2(i,j+1)*(1-x2(i,j+1)+x2(i,j-1)*(1-x2(i,j-1)-4*x2(i,j)*(1-x2(i,j)/4 7 elseif(i/=1.and.i/=100.and.j=1)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(i+1,j)*(1-x2(i+1,j)+x2(i-1,j)*(1-x2(i-1,j)+x2(i,j+1)*(1-x2(i,j+1)+x2(i,100)*(1-x2(i,100)-4*x2(i,j)*(1-x2(i,j)/4 8 elseif(i/=1.and.i/=100.and.j=100)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(i+1,j)*(1-x2(i+1,j)+x2(i-1,j)*(1-x2(i-1,j)+x2(i,1)*(1-x2(i,1)+x2(i,j-1)*(1-x2(i,j-1)-4*x2(i,j)*(1-x2(i,j)/4 9 elseif(i/=1.and.i/=100.and.j/=1.and.j/=100)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(i+1,j)*(1-x2(i+1,j)+x2(i-1,j)*(1-x2(i-1,j)+x2(i,j+1)*(1-x2(i,j+1)+x2(i,j-1)*(1-x2(i,j-1)-4*x2(i,j)*(1-x2(i,j)/4 end if if(t=20000)then if(x3(i,j)-x2(i,j)0)then x1(i,j)=con elseif(x3(i,j)-x2(i,j)0)then ! x1(i,j)=con ! elseif(x3(i,j)-x2(i,j)0)then ! x1(i,j)=