Compute Γ across strikes K∈{80,90,100,110,120} for T=0.5. Plot (or tabulate) and identify where gamma is maximised.
For K=100, compute Γ across expiries T∈{1.0,0.5,0.1,0.01}. What trend do you see?
Verify by finite difference. For S=100,K=100,T=0.5, compute Γ numerically via the central difference (Δ(S+h)−Δ(S−h))/(2h) with h=0.01. Compare to closed form.
Call–put gamma equality. Compute put gamma (by direct differentiation of put formula) at S=K=100,T=0.5. Verify it equals call gamma, as required by put-call parity.
Hint
Put formula: P(S,t)=Ke−r(T−t)Φ(−d2)−Se−q(T−t)Φ(−d1).