Let’s start with a simpler problem. Select 26 cards at random. If the number of HCPs in the 26 cards is exactly 7, what is the chance they include Ace and King of the same suit?
Start by enumerating all the ways to get seven HCP. There aren’t that many. Let’s break them down according to whether there are 2, 3, 4 or 5 honor cards. (The indented lines are the lines I input to the bc calculater.)
A,K
p = c(36,24) * 16
A,Q,J or K,K,J or K,Q,Q
q = c(36,23) * (64+24+24)
A,J,J,J or K,Q,J,J or Q,Q,Q,J
r = c(36,22) * (16 + 96 + 16)
K,J,J,J,J or Q,Q,J,J,J
s = c(36,21) * (4 + 24)
p counts the number of ways to have 24 non-honors (there are 36 non-honors altogether) multiplied by the number of ways to have an Ace and a King (4 Aces, 4 Kings yields 4x4). Similarly, q is the number of ways to have 23 non-honors and either A,Q,J (444 ways) or K,K,J (6 ways to select a pair of Kings times 4 Jacks) or K,Q,Q (also 6*4). After computing these numbers, calculate
p / (p + q + r + s) / 4
(We divide by 4 to change “probability of A and K” to “probability of A and K in same suit.”)
This gives .00543816 if I’ve made no mistake, or about 0.5%.
This assumes no other information. (Just the fact that your hands are balanced changes the availability of non-honors.)