gaussjac.m

jacobi quadrature, to resolve linear system with jacobi.

%%%%%%%%%%%%%%
%Top level script for calling getting the weights and abscissae of
%n-point Gauss-Jacobi quadrature

%First get upper and lower integration limits
lower_integration_limit = 0;
upper_integration_limit = 1;
number_points = 3;

%Then call the function to determine weights and abscissae
[abs,wgts] = Jacobi(number_points,lower_integration_limit,upper_integration_limit);

%Output results to files
fid = fopen('output.txt','w');
fprintf(fid,'%% Abscissae   Weights\n');
for i=1:length(abs)
  fprintf(fid,'%10.7f',abs(i));
  fprintf(fid,'%10.7f\n',wgts(i));
end
fclose(fid);