added function that computes score equations using the constraint problem with the concentration matrix K