#!/usr/bin/octave

# Rács definiálása, x koord -2*pi-től 2*pi-ig, 0.2-es lépésekben,
# xp (xpont) koord -3-tól 3-ig, szintén 0.2-es lépésekben
[x,xp]=meshgrid(-2*pi:.2:2*pi,-3:.2:3);
# Felcímkézzük az ábra tengelyeit
xlabel("{/Symbol f}");
ylabel("{/Symbol f}^.");
hold on;
# A fázisdiagrammon a sebesség megadása: 
v=xp;
# gyorsulás megadása: az ingánál x a szög, F=-m*g*l*sin(x)
# alkalmas mértékegységválasztással az állandók 1-nek vehetők
a=-sin(x);
# kis vektorocskákkal ezt kirajzoljuk
quiver(x,xp,v,a);
#a következő rajz előtt nem törlünk:
hold on;
# A harmonikus oszcillátor energiája:
E = 1/2*xp.^2 + (1-cos(x));
contour(x,xp,E);
# Elmentjük az ábrákat:
print("ingafazis.ps","-dps","-color","-landscape");
