program main integer:: N,k integer, allocatable:: F(:,:),A(:,:) real:: eps,epsp1 real, allocatable:: x(:),y(:),b(:,:),r(:,:) real, parameter:: rad=3.1415926/180 print*, "please entry the N." read*, N allocate(F(1,1:N**2),A(1:N**4,1:N**4),x(1:N),y(1:N),b(1:N,1:N),r(1,N**2)) A=0 print*, " " do i=1,N**2 A(i,i)=4 A(i+3,i)=-1 A(i,i+3)=-1 end do do i=2,N**2,3 A(i-1,i)=-1 A(i,i-1)=-1 A(i+1,i)=-1 A(i,i+1)=-1 end do do i=1,N**2 print "(10I5)",(A(i,j),j=1,N**2) end do
