Add Marian's code from last week
Paste in code, plot local