function [v] = gauss2(f,a,b)
% Guassian integration using two nodes

a0 = 1.0;
a1 = 1.0;
x0 = -sqrt(3.0)/3.0;
x1 = -x0;
alpha = (a+b)/2;
beta = (b-a)/2;
f0 = f(alpha + beta*x0);
f1 = f(alpha + beta*x1);

v = beta*(a0*f0 + a1*f1);
end