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	enderton	1.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';