Scipy

1.二维曲面画图

*#二维曲线画图
#plot(表达式,变量取值范围，属性=属性值）
#plot(表达式1，表达式2，变量取值范围，属性=属性值）
#plot((表达式1，变量取值范围1），（表达式2，变量取值范围2))

#在同一个图形界面画出 y1=2sinx,x∈[-5，5],y2=cos(x+π/4),x∈[-7，7]

from sympy.plotting import plot
from sympy.abc import x,pi  #引进符号变量x,常量pi
from sympy.functions import sin,cos

plot((2*sin(x),(x,-5,5)),(cos(x+pi/4),(x,-7,7)))

运行结果：

2.画三维曲面

#三维曲面画图
#画出三维曲面z=sin((x^2+y^2)^(1/2))

from pylab import rc  #pylab为matplotlib的接口
from sympy.plotting import plot3d
from sympy.abc import x,y
from sympy.functions import sin,sqrt
rc('font',size=16)
plot3d(sin(sqrt(x**2+y**2)),(x,10,10),xlabel='$x$',ylabel='$y$')

运行结果

3.隐函数作图

#隐函数作图
#绘制(x-1)^2+(y-2)^3-4=0的图形

#方法1
from pylab import rc
from sympy import plot_implicit as pt,Eq
from sympy.abc import x,y  #引进符号变量
rc('font',size=16)
#rc('text',usetex=True)
pt(Eq((x-1)**2+(y-2)**3,4),(x,-6,6),(y,-4,4),xlabel='$x$',ylabel='$y$')


#Eq(x,y):x=y

运行结果：

隐函数作图，第二种方法

#匿名函数做隐函数
from sympy import plot_implicit as pt
from sympy.abc import x,y
ezplot=lambda expr:pt(expr)
ezplot((x-1)**2+(y-2)**3-4)