function I = Gauss(a,b,n,p,w)

r=zeros(5,5);
c=zeros(5,5);

r(2,1:2)=[-.5773502692,.5773502692];
r(3,1:3)=[-.7745966692,0,.7745966692];
r(4,1:4)=[-.8611363116,-.3399810436, .3399810436, .8611363116];
r(5,1:5)=[-.9061798459,-.5384693101,0,.5384693101,.9061798459];

c(2,1:2)=[1,1];
c(3,1:3)=[.5555555556,.8888888889,.5555555556];
c(4,1:4)=[.3478548451,.6521451549,.6521451549,.3478548451];
c(5,1:5)=[.2369268850,.4786286705,.5688888889,.4786286705,.2369268850];

I=0;
for i=1:n
    x=((b-a)*r(n,i)+(b+a))/2;
    I = I + c(n,i)*func(x,p,w)*(b-a)/2;
end