function sirkulaerpolarisering

clear all;
% Gir verdier til pi, lambda, k, T og omega:
pi = 3.141592654;
lambda = 1.0;
k = 2*pi/lambda;
T = 1.0;
omega = 2*pi/T;
% amplituder x0 og y0
y0 = 1.0;
z0 = -2.0;
% N = antall punkter, med N=36 tegnes en pil og
% et punkt for hver 10. grad
N = 36;
% plotter i yz-planet, dvs ved x = 0
x = 0.0;

% Animering av hvordan spissen av D-vektor roterer
% i yz-planet som funksjon av tiden
figure;
axis_range = 1.1*max(abs(y0),abs(z0));
axis([-axis_range axis_range -axis_range axis_range]);
hold on;
for i=1:N
    t = i*T/N;
    y = y0*cos(k*x-omega*t);
    z = z0*sin(k*x-omega*t);
    plot(y,z,'.b');  % b=blå farge 
    compass(y,z);
    pause(0.15);    % Forsinker fremvisningen av D-vektor (juster denne for din pc)
end