% trajectory in momentum space.

global X Y
%X=[]; Y=[]; n=0;

cP2plotpack(n);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];
hold on
cP2plotpack(n+100);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];
cP2plotpack(n+200);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];
cP2plotpack(n+300);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];
cP2plotpack(n+400);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];
cP2plotpack(n+500);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];
cP2plotpack(n+600);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];
cP2plotpack(n+700);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];
cP2plotpack(n+800);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];
cP2plotpack(n+900);
X = [X expect(xop(pdata),pdata)];
Y = [Y expect(yop(pdata),pdata)];

plot(X,Y);
hold off
title(sprintf('particle propagation (momentum space)'));

n=n+1000;
