Correlation functions by cluster variation method for Ising model with NN, NNN, and plaquette interactions