%make a tabulated angle_harmonic potential

howmany=1000
k=10
phi0=3*pi/4
phi=0:pi/howmany:pi';
V=0.5*k*(phi-phi0).^2;
F=k*(phi0-phi);
plot(phi,V,phi,F)
data(:,1)=phi;
data(:,2)=F;
data(:,3)=V;