用Python代码写一个程序判断4阶矩阵是否对称，并统计矩阵中素数的个数。

可以使用以下代码实现：

import math

def is_prime(num):
    if num <= 1:
        return False
    for i in range(2, int(math.sqrt(num))+1):
        if num % i == 0:
            return False
    return True

def count_primes(matrix):
    count = 0
    for row in matrix:
        for num in row:
            if is_prime(num):
                count += 1
    return count

def is_symmetric(matrix):
    for i in range(len(matrix)):
        for j in range(i+1, len(matrix)):
            if matrix[i][j] != matrix[j][i]:
                return False
    return True

# example matrix
matrix = [
    [1, 2, 3, 4],
    [2, 5, 6, 7],
    [3, 6, 8, 9],
    [4, 7, 9, 10]
]

if is_symmetric(matrix):
    print("Matrix is symmetric")
else:
    print("Matrix is not symmetric")

print("Number of primes in matrix: ", count_primes(matrix))

这个程序使用了两个函数，一个用于计算一个数是否为素数，另一个用于计算矩阵中素数的个数。判断矩阵是否对称的函数使用了两个嵌套的for循环来遍历矩阵中的元素，并检查对应的元素是否相等。