Diagnostics/DiagUtility/ini_poly3.m

function [poly3] = ini_poly3()
% P=INI_POLY3
%
% Reads the file 'POLY3.COEFFS' and returns coefficients in P

fid=fopen('POLY3.COEFFS','r');

n=fscanf(fid,'%i',1);

for k=1:n,
 a=fscanf(fid,'%g',3);
 T(k)=a(1);
 S(k)=a(2);
 D(k)=a(3);
 poly3(k).t=a(1);
 poly3(k).s=a(2);
 poly3(k).dens=a(3);
end
for k=1:n,
 a=fscanf(fid,'%g',9);
 P(k,:)=a';
 poly3(k).coeffs=a';
end
P(:,10)=T';
P(:,11)=S';
P(:,12)=D';
1	function [poly3] = ini_poly3()
2	% P=INI_POLY3
3	%
4	% Reads the file 'POLY3.COEFFS' and returns coefficients in P
5
6	fid=fopen('POLY3.COEFFS','r');
7
8	n=fscanf(fid,'%i',1);
9
10	for k=1:n,
11	a=fscanf(fid,'%g',3);
12	T(k)=a(1);
13	S(k)=a(2);
14	D(k)=a(3);
15	poly3(k).t=a(1);
16	poly3(k).s=a(2);
17	poly3(k).dens=a(3);
18	end
19	for k=1:n,
20	a=fscanf(fid,'%g',9);
21	P(k,:)=a';
22	poly3(k).coeffs=a';
23	end
24	P(:,10)=T';
25	P(:,11)=S';
26	P(:,12)=D';