ну вот если взять это(я выкинул 0 и 1 константы для простоты) In [42]: def ro0(x,v): ...: return math.exp(-x**2 - v**2)/(math.pi**2) то можно это переписать в таком виде In [48]: def ro0_vec(x,v): ...: return np.exp(-x**2 - v**2)/(math.pi**2) a = np.arange(0, 5) b = np.arange(0, 5).reshape((5, 1)) In [58]: ro0_vec(a, b) Out[58]: array([[ 1.01321184e-01, 3.72739804e-02, 1.85576221e-03, 1.25040274e-05, 1.14021971e-08], [ 3.72739804e-02, 1.37123311e-02, 6.82696765e-04, 4.59997462e-06, 4.19463390e-09], [ 1.85576221e-03, 6.82696765e-04, 3.39894705e-05, 2.29019251e-07, 2.08838525e-10], [ 1.25040274e-05, 4.59997462e-06, 2.29019251e-07, 1.54311957e-09, 1.40714291e-12], [ 1.14021971e-08, 4.19463390e-09, 2.08838525e-10, 1.40714291e-12, 1.28314824e-15]]) это то что нужно было?
Обсуждают сегодня