degrees_arr = [0.087266463 0.174532925 0.261799388 0.34906585 0.436332313 0.523598776 0.610865238 0.698131701 0.785398163 0.872664626 0.959931089];
visibility_arr = [0.936170213 0.914893617 0.870967742 0.784946237 0.677419355 0.548387097 0.397849462 0.268817204 0.086956522 0.078651685 0.288888889];
colors_arr = [[1,0,0] [0,1,0] [0 ,0, 1] [0, 1, 1] [1, 0 ,1] [1, 1, 0] [0, 0 ,0] [0.8500 ,0.3250 ,0.0980] [0.4940 ,0.1840, 0.5560] [0.4660, 0.6740 ,0.1880] [0.6350, 0.0780, 0.1840]];
length_arr = length(visibility_arr);

theta = linspace(0, 2*pi);

for i = 1:length_arr
    x = visibility_arr(i) * cos(degrees_arr(i)) * cos(theta);
    y = visibility_arr(i) * sin(degrees_arr(i)) * sin(theta);

    plot(x,y, 'Color' ,colors_arr(3 * i - 2 : 3 * i), 'LineWidth',3);
    title('Ellipses of the measured angles of  Quarter-wave plate and the Visibility calculated');
    xlabel('X = calculated visibility \cdot cos(angle) \cdot cos(t)')
    ylabel('Y= calculated visibility \cdot sin(angle) \cdot sin(t)')

    axis equal;
    hold on;
end