from math import *
import matplotlib.pyplot as plt
import numpy as np
from matplotlib.patches import Polygon

def f(x):
    return 2*x + 3

def mTrapez(a,b,n):

    wysokosc = (b-a)/n
    suma = 0
    podstawa_a = f(a)


    for i in range(1, n):
        podstawa_b = f(a + i*wysokosc)
        p = 0.5 * (podstawa_a + podstawa_b) * wysokosc
        suma = suma + p
        podstawa_a = podstawa_b
    return round(suma,2)
a = int(input("Podaj dolna granice przedzialow: "))
b = int(input("Podaj gorna granice przedzialow: "))
n = int(input("Podaj ilosc trapezow: "))

#print('Pole = ', mTrapez(a,b,n))

x1 = np.arange(a, b, 0.1)
x2 = np.arange(a, b, 0.1)

plt.figure('Pole wynosi:   ' + str(mTrapez(a, b , n)))
plt.xlabel('X')
plt.ylabel('Y')
plt.plot(x1, f(x1), 'r')
ix = np.linspace(a,b)
plt.show()