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).