python中积分的导数

2024-06-01 12:07:41 发布

您现在位置:Python中文网/ 问答频道 /正文

我有一个函数,有两个变量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}$

我做错什么了?提前谢谢。在


Tags: importinputrawreturndefasnpplt