Batch Normalization Formula

Spot the missing square root and epsilon in this batchnorm.

import numpy as np

def batchnorm(x):
    mean = x.mean(axis=0)
    var = x.var(axis=0)
    return (x - mean) / var

What is the bug in this batch normalization?