Subroutine solver(ismax,jsmax,ksmax)
		Use global
	implicit none
	integer :: i,j,k
	integer,intent(in) :: ismax,jsmax,ksmax
	real(kind=8) :: dx,dy,dz
	real(kind=8) :: xx,yy,zz

	dx = (xu-xl)/(ismax-1.0)
	dy = (yu-yl)/(jsmax-1.0)
	dz = (zu-zl)/(ksmax-1.0)

		do k = 1,ksmax
			zz = zl + (k-1.0)*dz
			do j = 1,jsmax
				yy = yl+(j-1.0)*dy
				do i = 1,ismax
				xx = xl+(i-1.0)*dx
					w(i,j,k) = zz*zz + yy*yy + xx*xx
					x(i,j,k) = xx
					y(i,j,k) = yy
					z(i,j,k) = zz

				enddo
			enddo
		enddo
	
	end subroutine solver