双线性变换部分代码

import numpy as np
import math
def Cacul_power(pnumA,pnumB):
    a=np.zeros((pnumA+1,1))
    for i in range(pnumA+1):
        a[i]=math.factorial(pnumA)*np.power(-1,i)/(math.factorial(i)*math.factorial(pnumA-i))
        print a[i]
    b=np.zeros((pnumB+1,1))
    for j in range(pnumB+1):
        b[j]=math.factorial(pnumB)/(math.factorial(j)*math.factorial(pnumB-j))
        print b[j]
    c=np.zeros((pnumA+1,pnumB+1))
    print c
    for i in range(pnumA+1):
       for j in range(pnumB+1):
        c[i][j]=a[i]*b[j]
    d=np.zeros((pnumA+1,pnumA+pnumB+1))
    for i in range(pnumA+1):
        for j in range(pnumB+1):
            d[i][j+i]=c[i][j]
    print d
    e=np.zeros((pnumA+pnumB+1,1))
    for i in range(pnumA+1):
        for j in range(pnumA+pnumB+1):
            e[j]+=d[i][j]
    for i in range(pnumA+pnumB+1):
        e[i]=np.power(3.14/2,i)*e[i]
    print e
Cacul_power(2,6)