function y = DA_Funktion(x)

y = 12/(1+tan(x).^2) + sqrt((12/(1+tan(x).^2).^2) - (108/(1+tan(x).^2)))

for x = fliplr(1:60)
    
   y = fzero(12/(1+tan(x).^2)+sqrt((12/(1+tan(x).^2).^2) - (108/(1+tan(x).^2))),0);
   
   plot (y);
   
   hold on
end
