function [x,w] = gausscos();
%GAUSSCOS [X,W] = GAUSSCOS  returns the nodes  X  and weights  W
%	for 16-point Gauss quadrature on [-PI/2, PI/2] 	with weight
%	function  COS(X).

%	Reference:  Evans, G. (1993). Practical numerical integration.
%	John Wiley, Chichester.
%
%	GKS  29 May 98

x=[
-1.5327507132362304779
-1.4446014873666514608
-1.3097818904452936698
-1.1330068786005003695
-0.92027786206637096497
-0.67861108097560545347
-0.41577197673418943962
-0.14003444424696773778
 0.14003444424696773778
 0.41577197673418943962
 0.67861108097560545347
 0.92027786206637096497
 1.1330068786005003695
 1.3097818904452936698
 1.4446014873666514608
 1.5327507132362304779
];

w=[
0.0024194677567615628193
0.014115268156854008264
0.040437893946503669410
0.083026647573217742131
0.13834195526951273359
0.19741148870253455567
0.24763632094635522403
0.27661095764826050408
0.27661095764826050408
0.24763632094635522403
0.19741148870253455567
0.13834195526951273359
0.083026647573217742131
0.040437893946503669410
0.014115268156854008264
0.0024194677567615628193
];
