k1= [pi/3];
k2=[0.7105, 1.2972];
k3=[];

kv=k2;

xx=[0:pi/100:pi];
plot(xx,zeros(length(xx),1));hold on;

plot(kv,zeros(length(kv),1),'*')
