因此,我想绘制简单的伽玛函数,但是有一些问题。我的代码是:
#!/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/