1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
% Copyright ©2016 The Gonum Authors. All rights reserved.
% Use of this source code is governed by a BSD-style
% license that can be found in the LICENSE file.
clc
clear all
close all
% Generate herm points
min = 216;
max = 216;
fprintf('xCache = [][]float64{\n')
for i = min:max
[x,w] = hermpts(i);
fprintf('{')
for j = 1:i-1
fprintf('%1.16e, ',x(j))
end
fprintf('%1.16e},\n',x(i))
end
fprintf('}\n')
fprintf('wCache = [][]float64{\n')
for i = min:max
[x,w] = hermpts(i);
fprintf('{')
for j = 1:i-1
fprintf('%1.16e, ',w(j))
end
fprintf('%1.16e},\n',w(i))
end
fprintf('}\n')
|
