python:画由抛物线: y^2=2x 与直线 y=x-4 所围成的图形

《高等数学》同济大学版 P339

编写 test_diff_3_area.py 如下

python 复制代码
# -*- coding: utf-8 -*-
""" 画由抛物线: y^2=2x 与直线 y=x-4 所围成的图形 """
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.patches import Polygon

def fun(x):
    return x-4

def fun1(x):
    return np.sqrt(2*x)

def fun2(x):
    return -1*np.sqrt(2*x)


x = np.linspace(0, 10, num=100)
y = fun(x)
y1 = fun1(x)
y2 = fun2(x)

fig, ax = plt.subplots()
plt.plot(x, y,  'r', linewidth=2)
plt.plot(x, y1, 'r', linewidth=2)
plt.plot(x, y2, 'r', linewidth=2)

a = 0.0
b = 8.0
# 坐标轴设置
ax.set_xticks([0, 2, 4, 6, 8])
ax.set_yticks([])
#ax.set_xticklabels(['$a$', '$b$']) # 换成公式字体
plt.figtext(0.98, 0.05, '$x$')
plt.figtext(0.01, 0.98, '$y$') #0~1代表在图的比例处

# 绘制灰色多边形 y轴>0
ix = np.linspace(4, 8.0)
iy = fun(ix)
ixy = zip(ix, iy)
verts  = [(0, 0)] + list(ixy) + [(0, 4)]
# 多边形 Polygon
poly = Polygon(verts, facecolor='0.9', edgecolor='0.3')
ax.add_patch(poly)

ix = np.linspace(0, 8.0)
iy1 = fun1(ix)
ixy1 = zip(ix, iy1)
verts1 = [(0, 0)] + list(ixy1) + [(0, 4)]
# 多边形 Polygon
poly1 = Polygon(verts1, facecolor='1.0', edgecolor='0.3')
ax.add_patch(poly1)

# 绘制灰色多边形 y轴<0
ix = np.linspace(0.0, 2.0)
iy2 = fun2(ix)
ixy2 = zip(ix, iy2)
verts2 = [(0, 0)] + list(ixy2) + [(2, 0)]
# 多边形 Polygon
poly2 = Polygon(verts2, facecolor='0.9', edgecolor='0.3')
ax.add_patch(poly2)

ix = np.linspace(2, 4.0)
iy = fun(ix)
ixy = zip(ix, iy)
verts  = [(2, 0)] + list(ixy) + [(4, 0)]
# 多边形 Polygon
poly = Polygon(verts, facecolor='0.9', edgecolor='0.3')
ax.add_patch(poly)

# 添加 LaTex数学公式
x_math = 3
y_math = 3
latex = r'$\int_{-2}^4 (y+4 - \frac{1}{2} y^2)dy $'
plt.text(x_math, y_math, latex, fontsize=14, horizontalalignment='center')
plt.grid()
plt.show()

运行 python test_diff_3_area.py

相关推荐
郭庆汝4 小时前
pytorch、torchvision与python版本对应关系
人工智能·pytorch·python
思则变7 小时前
[Pytest] [Part 2]增加 log功能
开发语言·python·pytest
漫谈网络7 小时前
WebSocket 在前后端的完整使用流程
javascript·python·websocket
try2find8 小时前
安装llama-cpp-python踩坑记
开发语言·python·llama
博观而约取9 小时前
Django ORM 1. 创建模型(Model)
数据库·python·django
精灵vector11 小时前
构建专家级SQL Agent交互
python·aigc·ai编程
Zonda要好好学习11 小时前
Python入门Day2
开发语言·python
Vertira11 小时前
pdf 合并 python实现(已解决)
前端·python·pdf
太凉11 小时前
Python之 sorted() 函数的基本语法
python
项目題供诗12 小时前
黑马python(二十四)
开发语言·python