load resu.mat
figure(1)
plot(dx, dFi_MP)
xlabel('\delta \tau, m');
ylabel('\phi_{MP}');