assert_almost_equal(
            np.var(s1),
            np.dot(weights, vars) + np.dot(weights, means**2) - np.dot(weights, means)**2,
            1