function x=pMapInv(x,gamma)

delta = 0.25;


f = @(t,x)( -[x(2);x(1)-x(1).^3-delta*x(2)+ gamma*cos(t)] );

[t,y]=ode45(f,[0:(pi/10):2*pi],x);

x = y(end,:)';