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).