# -*- coding: utf-8 -*-
"""
Created on Thu Sep  8 18:37:08 2022

@author: Chris
Lab7 Q3
"""

import numpy as np
import matplotlib.pyplot as plt
import math as maths

def f(x):
    return np.sin(x)**2

pi = maths.pi

a,b = 0,pi

N = 10000

p = (b-a)*np.random.uniform(size=N)

I = (b-a)/N*sum(f(p))

print(I)