parameter N=100
                  real BDA(N+1),BDB(N+1),BDC(N+1),BDD(N+1)
                  real F(N+1,N+1)
                  real,allocatable::W(:)
                  character*20 frmt 
                  open (10, file='laplace.dat')
		  do i=1,n
		  do j=1,n
                    f(i,j)=0
		  enddo
		  enddo
                  id=4*(N+1)+(13+int(log(float(N+1))/log(2.)))*(N+1)
		    
		  do i=1,n+1
		  bda(i)=0.
		  bdb(i)=0.
		  bdc(i)=0.
		  bdd(i)=0.
		  enddo

                  do i=n/2-10,n/2+10
		  do j=n/2-10,n/2+10
                    f(i,j)=-5
                  enddo
		  enddo
		  
                  allocate(W(id))

                  CALL HWSCRT(0.,1.,N,1,BDA,BDB,0.,1.,N,1,BDC,BDD,0.,F,
     *            n+1,PERTURB,IERROR,W)
    		  write(frmt,*) "(",n+1,"f15.5)"
!    		  print *, frmt
                  write(10,frmt) F
				
		  stop
		  end