编写程序，把n×n矩阵a加上矩阵a的转置，存放在矩阵b中。

可以使用以下代码实现：

n = 3  # 矩阵的维度
a = [[1, 2, 3], [4, 5, 6], [7, 8, 9]]  # 原始矩阵

# 创建一个n×n的全0矩阵b
b = [[0 for j in range(n)] for i in range(n)]

# 计算a和a的转置的和，并存放在矩阵b中
for i in range(n):
    for j in range(n):
        b[i][j] = a[i][j] + a[j][i]

# 输出结果
print("矩阵a：")
for row in a:
    print(row)

print("矩阵b：")
for row in b:
    print(row)

输出结果为：

矩阵a：
[1, 2, 3]
[4, 5, 6]
[7, 8, 9]
矩阵b：
[2, 6, 10]
[6, 10, 14]
[10, 14, 18]

其中，矩阵b中的每个元素都是矩阵a中对应位置的元素和其转置位置的元素之和。