from math import cos
while True:
a = float(input('a: '

)
Z1 = cos(a)+cos(2*a)+cos(6*a)+cos(7*a)
Z2 = 4*cos(a/2)*cos(2.5*a)*cos(4*a)
print(Z1, Z2, Z1-Z2, sep = '\n'

Вот что получается при целых а от 1 до 10:
Разница в значениях Z1 и Z2 получается или нулевая, или порядка машинной точности для восьмибайтных вещественных чисел типа float, так что формула с точки зрения практических компьютерных вычислений при данной выборке значений переменной а оказывается верной. Кстати, результаты не должны совпадать абсолютно точно!
А строго данную формулу можно доказать лишь чисто математически и компьютерное программирование здесь вообще ни при чём!