n1=1 ;
    n2=1.5;

    a = pi/6;
    a1 = linspace(-pi/2,pi/2);
    a2 = asin(((n1.*sin(a1))./n2));
    a3 = asin((n2.*sin(a-a2))/n1);
    a4 = (a1+a3-a);
    a5 = (a4(imag(a4)==0));
    x = linspace(-pi/2, pi/2, 60);



    plot(x, a5);
    grid on


    %% zinkdoppat kronglas


    z = @(x, a1) asin(((1.517+10.72./x).*sin(pi./6-asin((sin(a1)./(1.517+10.72./x))))))

    a  = pi./6;
    a1 = linspace(-pi/2,pi/2);

    a2 = @(x, a1) (a1+z(x, a1)-pi./6)

    z1 = a2(480, a1);
    z11 =(z1(imag(z1)==0));
    z2 = a2(530, a1);
    z22 =(z2(imag(z1)==0));
    z3 = a2(580, a1);
    z33 = (z3(imag(z1)==0));
    z4 = a2(630, a1);
    z44 = (z4(imag(z1)==0));

    x = linspace(-pi/2, pi/2, 59);


    grid on

    hold on
    plot(x, z11, 'blue');
    plot(x, z22, 'green');
    plot(x, z33, 'yellow');
    plot(x, z44, 'red');



    %% tungt flintglas


    z = @(x, a1) asin(((1.653+10.27./x).*sin(pi./6-asin((sin(a1)./(1.653+10.27./x))))))

    a  = pi./6;
    a1 = linspace(-pi/2,pi/2);

    a2 = @(x, a1) (a1+z(x, a1)-pi./6)

    z1 = a2(480, a1);
    z11 =(z1(imag(z1)==0));
    z2 = a2(530, a1);
    z22 =(z2(imag(z1)==0));
    z3 = a2(580, a1);
    z33 = (z3(imag(z1)==0));
    z4 = a2(630, a1);
    z44 = (z4(imag(z1)==0));

    x = linspace(-pi/2, pi/2, 56);


    grid on

    hold on
    plot(x, z11, 'blue');
    plot(x, z22, 'green');
    plot(x, z33, 'yellow');
    plot(x, z44, 'red');