Write a program that determines the Riemann integral of a function f on the interval [0,1] in different ways.

Write a program that determines the Riemann integral of a function f on the interval [0,1] in different ways.
Output the results of the calculations and the plot for the function f= sin(πx) as an example.  First divide the interval into n equal sub intervals 

 and then determine the values of the corresponding upper and lower sums 

 Now calculate the integral of f symbolically using SymPy and numerically using np.trapz() and again output the results using the example f= sin(πx).Plot the function f and the staircase functions belonging to the upper and lower sum using Matplotlib (all together in one plot).