我有一个函数,有两个变量f(x,y)=x+2y,我想把它与x积分,然后对y求导,然后用y绘制结果,这是我的代码:
import numpy as np
import matplotlib.pyplot as plt
import math
import scipy.integrate
from scipy.misc import derivative
x = raw_input("enter x: "); x=float(x)
y = raw_input("enter y: "); y=float(y)
def f1(x,y):
r=x+2*y
return r
r=f1(x,y)
print "r=",r
def int(y):
return scipy.integrate.quad(f1, 0,x,args=(y))
vint=np.vectorize(int)
ylist=np.linspace(0,y,10)
results,errs=vint(ylist)
plt.plot(ylist,results)
plt.show()
def dI(y):
return derivative(int,y)
vdI=np.vectorize(dI)
sol,err=vdI(ylist)
plt.plot(ylist,sol)
plt.show()
但我收到了一条错误信息:
^{pr2}$我做错什么了?提前谢谢。在
目前没有回答
相关问题 更多 >
编程相关推荐