One approach is deriving the variance of Beta = (Beta_0, Beta_1, ..., Beta_p) in the matrix form, as

sigma^2 * (X' * X)^{-1}, **(1)**

where X = (1, x) is an N-by-(p+1) matrix. Then you would rewrite formula (1) for the case of simple linear regression (p = 1).