function levelCurvesTrajectory(pThreshold,pThreshold2,xSteps,t,y0)
%Plota curvas de nível de acordo com a trajetória
%O Valor de threshold vai de ditar se aparecerão mais curvas na parte
%interna ou externa ao cilco limite. Quanto maior o threshold, mais curvas
%externas

%Trajetórias passando pelas curvas de nível
x2_i = sqrt(5); 
x1_i = 10^0.25;
v=@(x1,x2) (x1.^4 + 2*x2.^2 - 10).^2;
[~,Y]=ode23(@rigid,t,y0); 
x=Y(:,1);
y=Y(:,2);
[x1,x2] = meshgrid(-(pThreshold+x1_i):xSteps:(pThreshold+x1_i),-(pThreshold2+x2_i):xSteps:(pThreshold2+x2_i));
z = v(x1,x2);

figure1 = figure;
axes1 = axes('Parent',figure1);
set(axes1,'FontName','Arial','FontSize',12,'FontWeight','bold');
contour(x1,x2,z,'b','ShowText','on');
xlabel({'x_1'},'FontWeight','bold','FontSize',12,'FontName','Arial');
zlabel({'V(x)'},'FontWeight','bold','FontSize',14,'FontName','Arial');
ylabel({'x_2'},'FontWeight','bold','FontSize',12,'FontName','Arial');
hold on;
phasePlane(t,y0);
plot(x,y,'color','green','linewidth',2.5);
hold off;