Stabilization matrix method

The answer to the SMM exercise can be found at

smm_ans.c
smm_mc_ans.c