# -*- coding: utf-8 -*-
"""
Created on Sun Sep  4 14:18:14 2022

@author: Chris

Lab 6, Q5
"""

import numpy as np
import matplotlib.pyplot as plt
from scipy.special import p_roots

def f(x):
    f = np.exp(x**2) #x**5+4*x**3 + x**4+ 2*x**2
    return f
    


for n in range(1,12):
    # set up zeros and weights
    [xj,wj]=p_roots(n)
    
    I = sum(wj*f(xj))
    print(n,I)