I am trying to code a Gaussian Quadrature funtion. My function as of now is...
function I = GQ5_Riley(fun,a,b)
w = [0.5688888888888889 0.4786286704993665 0.4786286704993665 0.2369268850561891 0.2369268850561891];
x = [0.0000000000000000 -0.5384693101056831 0.5384693101056831 -0.9061798459386640 0.9061798459386640];
I = ((b-a)/2)*sum(w.(1:1:N))*fun(((b-a).*x+b+a)/2);
I need MatLab to sum the w values one at a time until it reaches the last and 5th w value. as well as multiply the b-a by the 5 x values.
Any help would be appreciated.