function [cd]=Drag_coef(v,c)
    M=v/c;
        if M<=0.9
            cd=0.75/sqrt(1-(v/c)^2);
        elseif M>0.9 && M<=1.09
             cd=1.73;
        else
             cd=0.75/sqrt(((v/c)^2)-1);
        end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% KOD PRO GRAFY %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% close all; clear; clc
% [t_C,p,rho,c] = USSA(0);
% for v=1:0.1:1000
%     M=v/c;
%         if M<=0.866
%             cd=0.75/sqrt(1-(v/c)^2);
%         elseif M>0.866 && M<=1.12
%              cd=1.5;
%         else
%              cd=0.75/sqrt(((v/c)^2)-1);
%         end
%     plot(M,cd,'k',Marker='.')
%     hold on
% end
% title("Koeficient odporu")
% xlabel("M")
% ylabel("c_D")