Moved simulation and plotting into separate functions. Working on matlab and octave.
Parts 1&2 working on matlab and octave after lots of debugging.

Parts of main deterministic.m have been moved to seperate functions
and deterministic.m can be started using doMesh=0|1,
doPart1=0|1, doPart2=0|1 to toggle what gets run.

Still a problem plotting semilog in matlab, for part 2,
possibly because of zeros? good enough for now.