%% EM Wave polarization
close all; clear all; clc
T=linspace(0,100,1000);
x=linspace(0,4*pi,100);

for j=1:length(T)
    t=T(j);

    for i = 1:9
        phi = pi*(i-1)/4;
        y=sin(x-t);
        z=sin(x-t+phi);
        subplot(3,3,i);
        plot3(x,y,z, 'linewidth', 2); axis equal tight
        hold on
        plot3(zeros(size(y)),y,z,'y','linewidth',4); hold off
        %title(sprintf('%d',phi))
        set(gcf,'color',[1,1,1]);
        view([-80,16])
        pause(.00000000001)
        
    end
end