def HaarTransform(N):
    wavelet = np.zeros((N,N))
    for i in range(N//2):
        wavelet[i,i*2] = sqrt(2)/2
        wavelet[i,i*2 + 1] = sqrt(2)/2
    for j in range(N//2):
        wavelet[N//2 + j,j*2] = -sqrt(2)/2
        wavelet[N//2 + j,j*2 + 1] = sqrt(2)/2
    return wavelet

wave600 = HaarTransform(600)
wave478 = HaarTransform(478)