1、20406080102040608010 -1.0-0.750-0.50-0.2500.00.2500.500.7501.020406080102040608010 -1.0-0.750-0.50-0.2500.00.2500.500.7501.03.20406080102040608010 0.900.92500.9500.97501.01.0251.051.0751.0.20406080102040608010 0.7600.7340.76080.7810.8150.84290.87030.89760.9250列列列.（一）前三幅图的程序program coupleimplicit non
2、ereal,parameter:e1=0.6,u=3.7integer,parameter:con=1real,dimension(100,100):x2,x3integer,dimension(100,100):x1integer i,j,t,time,kk=time()open(10,file=“output10.txt“)do j=1,100do i=1,100x2(i,j)=ran(k)end doend dodo t=1,20000 do j=1,100do i=1,1001 if(i=1.and.j=1)then x3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*
3、(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)/42 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)/43 elseif(i=1
4、00.and.j=1)thenx3(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)/44 elseif(i=100.and.j=100)thenx3(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
5、(i,j-1)-4*x2(i,j)*(1-x2(i,j)/45 elseif(i=1.and.j/=1.and.j/=100)thenx3(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)/46 elseif(i=100.and.j/=1.and.j/=100)thenx3(i,j)=u*x2(i,j)*(1-x2(i,j)+e1*u*(x2(1,j)*(1
6、-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)/47 elseif(i/=1.and.i/=100.and.j=1)thenx3(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)/48 elseif(i/=1.
7、and.i/=100.and.j=100)thenx3(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)/49 elseif(i/=1.and.i/=100.and.j/=1.and.j/=100)thenx3(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)+
8、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)/4end ifif(t=20000)thenif(x3(i,j)-x2(i,j)0)thenx1(i,j)=conelseif(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! end if! end ifend doend do x2=x3end do write(10,“(100f12.9)“)x3write(*,“(100f12.9)“)x3close(10)pauseend
