To reproduce the graphs in example #5, you can use this code.

You will also need general model solving code as well as this software.