function mylu(A)
	m,n=size(A)
	if m!=n
		println("Need a square matrix!")
		throw(exit())
	end
	U=zeros(size(A))
	L=zeros(size(A))
	Ak=copy(A)
	for k=1:n
		U[k,:]=Ak[k,:]
		L[:,k]=Ak[:,k]/U[k,k]
		Ak=Ak-L[:,k]*U[k,:]'
	end
	return L,U
end