utils/matlab/ini_poly3.m

b30af0ec84 Alis*0001 function [poly3] = ini_poly3()
                0002 % P=INI_POLY3
                0003 %
                0004 % Reads the file 'POLY3.COEFFS' and returns coefficients in P
                0005
                0006 fid=fopen('POLY3.COEFFS','r');
                0007
                0008 n=fscanf(fid,'%i',1);
                0009
                0010 for k=1:n,
                0011  a=fscanf(fid,'%g',3);
                0012  T(k)=a(1);
                0013  S(k)=a(2);
                0014  D(k)=a(3);
                0015  poly3(k).t=a(1);
                0016  poly3(k).s=a(2);
                0017  poly3(k).dens=a(3);
                0018 end
                0019 for k=1:n,
                0020  a=fscanf(fid,'%g',9);
                0021  P(k,:)=a';
                0022  poly3(k).coeffs=a';
                0023 end
                0024 P(:,10)=T';
                0025 P(:,11)=S';
                0026 P(:,12)=D';