compute_contrast_bcf_model function

Compute Contrast for BCF Model