Playing around with Matlab plot for DeltaP-vs-ln(r/DX) calculation