因此,我想绘制简单的伽玛函数,但是有一些问题。我的代码是:

#!/usr/bin/env python
# -*- coding: cp1250 -*-
#import math
from scipy.special import *
#from scitools.std import *
from pylab import *

def f1(x):
    return gamma(x)


x = linspace(-6, 6, 512)
y1 = f1(x)

# Matlab-style syntax:
plot(x, y1)

xlabel('x')
ylabel('y')
legend(r'$\Gamma(x)$')
grid(True)

show()


我尝试从数学和scipy.special导入gamma函数,但出现以下错误:


  追溯(最近一次通话最近):文件“ D:/faxstuff/3.godina/kvantna/plotgamma.py”,第13行,y1 = f1(x)文件“ D:/faxstuff/3.godina/kvantna/plotgamma .py“,第1行,在f1中返回gamma(x)文件“ mtrand.pyx”,行1599,在mtrand.RandomState.gamma中(numpy \ random \ mtrand \ mtrand.c:8389)ValueError:shape

怎么做?这应该很容易,但是我似乎失败了:(

最佳答案

其中一个模块(我认为是pylab)是由gamma随机变量函数遮盖了gamma函数。这行得通,但是我不得不关闭对图例的调用(我不确定为什么)。

from scipy.special import gamma as Gamma
#from scitools.std import *
from pylab import *

def f1(x):
    return Gamma(x)


x = linspace(-6, 6, 512)
y1 = f1(x)
gca().set_autoscale_on(False)

# Matlab-style syntax:
plot(x, y1)

xlabel('x')
ylabel('y')
# legend(r'$\Gamma(x)$')
axis([-6, 6, -100, 100])
grid(True)

show()

关于python - Python中的Gamma函数图,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/5685576/

10-12 23:35