I came across a MATLAB app called "Aircraft Intuitive Design" by Zac Lietzau:
https://www.mathworks.com/matlabcentral/fileexchange/66770-aircraft-intuitive-design-aid
It allows the user to see the changes made to the geometry reflected in the stability and aerodynamic characteristics in real time.
It also computes the profile drag coefficient, as seen in the image below (Drag: CD = 0.047 (profile) + 0.057CL^2 (induced)):
I know induced drag can be calculated using vortex lattice methods, but I want to understand how the profile drag is being calculated. (I am unable to access the code)