From 25da8b87e3bba5366731e9c9b9194e911f0052ee Mon Sep 17 00:00:00 2001
From: arcan1s <esalexeev@gmail.com>
Date: Tue, 16 Jul 2013 18:00:30 +0400
Subject: [PATCH] Added lib 'summary_stat.c'

---
 stat_new/stat           | Bin 9431 -> 8005 bytes
 stat_new/stat_print.c   |  16 +++++++--
 stat_new/summary_stat.c |  70 ++++++++++++++++++++++++++++++++++++++++
 3 files changed, 84 insertions(+), 2 deletions(-)
 create mode 100644 stat_new/summary_stat.c

diff --git a/stat_new/stat b/stat_new/stat
index f9733725ebbf9939b22a4a8d79ba9dfbc59eac0e..f12f41dee628393a25aafc2ce91cccec72926c31 100755
GIT binary patch
delta 3512
zcmccadDL!#gwO{T2L=#;kPHkRd=s^nIsV}Si7_z1Xa)uaumB?i0|Og~#SJ1D7#O${
zAPgw&!35z3OcrF65QHj{fJ#bmPqt!=6I#OwQL}~<V&@vp$*UL_iMad_a=6pZWc0`7
z>As#9i~G&2`X5eCV6vEef=Pl=eey;oaYlJCX*pStS$whsvjn3Ah(Ec2S%@_`C%-s#
z@)TyJ$-9}A85t%&Wmc~DfLO!uf(PP8kh@u+Vjzlx7b4Eiz`y|V7f1|*1(3v{-T?_S
zNFa$bgGC^O0+KkA`SlDMND{0NwP4Z!NgU)Tupk2igGcil4v%iuRt^OQkJbYvtS^GW
zf=4(G!&E@%|E5J83Jm{MOE?r5_~jiK{;Puc86e4*5B~rE|6etWLxG__1LUcf7r+8Z
z9108>AkVx!0Om)5_#jWb+yLeWf%qWDzgz(3dx7{MH@}?l|Nno(x1OCB<2*Vac{D%y
z5D?;d{J;hd2L_Mh2OYrFad{30h8OMs|NsAzHUZ?I9saxt45pzTx*og<3?9ut_&hqF
zmNUg3-T`tV$ZKC>_k*PVKj4>d0Ez7Y@xbaQYq4p%f_VQQLFF23H?k`*Fq8_0dUQU0
zQTPA<|JY-!^Vk&_7(X8dS?1Af8^Nl;z)&IwH(n%kzX12-8Ei839?fq8JUVMnbe10J
z3_ak{>AS<H({;mf*B>BTkGuW>QLk0HUH^E@@aPq31#vyPT{n0zUg&0MJy62d?fR$F
z^+&VoAI5IiKh~}XYOZ*6yPnv00_5+`<1ae@{r~UL>H5K=+x3G-=fPdV3=9nWK*7HM
z6evil>Ot~Qxg@yU1CU(&zyJTC4tV*Lfq?<L8$R=E9r(;I==$a}zn1F}5bMQfegP0S
z&h-U|`=s0TPj3Lobs!(`Ao<`#Hv@nD`a^INx?PX(uXjBF<D;4LpxgCN^8rSPY^3W0
zklA-W^9zLD;g7oXnP0&70e{p<kkwE4qYhwos{qI)1&pAe=)BnZX`eHv14F$>=kXn&
zWcI@8@BjbsVEn<~!~g&Pe~)h09o+&T(xdqROE&}K1&H;6u5Uns-3&191rQga1||fG
z2XK%Z!h<z|(}Cdy*WdsDLC(DLnLp0;O*FsA0sg2HAb|^D0Wkjn$c_6ZaXK)7oCnGk
zFP8iPJLiK(=kXU6fBydmC#92zK}oK&_K!zr=?Ab#?BO&$e)$$qGGUl(&1=H+^!MaK
z-bkn5P+woiph$(_5XTVD;1Ew|D^U4aqoAOWr=Xw!qWz$>g0ljQRe)gSDNO#rD^{<I
zA8;`!sG4U&u%(_+nu4mOo+*^%Vqn0p8B>;v!N=3j)iKCc)m+asjmtUG#o5Od%3)xX
zX12C~7WP+u{r?a0v&--Q|3erU7)<{B|G$QTfuZ={|Nmzg7#ODf`~Uw90|Nug|Ns9@
zCfD*QGrjmVc@dvpy~nTr{|i8Y#V6p#C*j4<UCz<KU@v8@WvrqE5|&|LU~pkzU}*UL
z|Gzyb9vt}u+L)Yq*_f@Ec-T1@K=LUJ3=H%C{{Jro(g2h91<6A6f*iKs@Bja)AbBS~
zfqo`eK8ZeNS3ZT_$>IFEoFE6D`1}7q)8vbc!dx3zA$86K*2!D~YEBW*%8wVMnSp^}
z3kO7663T~FZ_-dc2b7kA(j45Aiv;u-H6||<kY@~;d{DqZ09x6Ad;?N6VKSqjyyyvL
zh`FGu1|+%x#GiatN_?`9pa|oH$%%sUj5U)R1q+zC6DD6|6xEytb=Wc}y$wnqgVNWa
z^s~u~GQyJ|FtTh;5L&>*xMDJ+sJ@^E$aDq<1`SZoV_;wind~U4FIWQ=t$~X2O>Pv`
z=bQi)oiKT)s5uwdT?`Bi8<;2m5;YTqIr{`u%Y?~}V)|ULK#O3RTqR~E*Z_4(1E~67
zU|=Yid{9gu-mLfl6+OZ{nMquq6B@o-I3`<(n+rgB2cW_FW^$vrJ{MB6V3)X=;0tIg
z;RUqSATpUz!hjRnqF~^jY$ah1<8g3Lu0rqxxF_#I@Fch=|B^80gvQPbp2=2{=CHW>
zzypdC1_qDGjgtCY7SP(+f_L&RNi!}`T?CT$;GO(S(vHc4ce0g~9a97E<SHpUP7X1U
z6a&MG$p@v(IX8gPD+2?=)X9v}+Zi<`Uz8SQW@o6Gd{J7QaRQjf$*=)Tb1|F%)7%UZ
zlNn`1nRyr*z%(zz2M}G)2TzF*^C2WN13$w7XqZ7nAS5$`00S)4AR-WwnL&^N7JLv9
z2+7PKglM=yguo;-gD?XJC}iuwTm-?)Ai^L3l0XQ7S<DQg46w!>SO7{eGl;=kT~Hw~
z#mpei@BtbHU;!w>%pk!43wDqYABX^9W(G+H58la5vgPR3fSRLB47>~qlLh6(>jR+K
z8l(q=K}|o9>;yFNIB-)BWG1BP2db(;>T|*3P~U?Dpm+)c0|Uq$Sg`>WW>^4L4-Z99
zVGohoKlz}nzVbP+v>?L|E=abA$T2Y70gHpfa`I+b)%q7;2_Z<JgV_+`Cs-Y3*fN3L
z3%3@cn1O*olo8Z|#SB|*Mg|5U23`hGYYn6rRw!9Bg4^?xFEa`nhJ)Ft;SW-q#K^!P
zz#zc@OaCA-5YC*uP+nZJ7Hl$RGHnLSE`sVsj#H4{^&nZ+BOn_=rE&dbMg|5!CJ6=y
zXqE%10pSNAaiqi!EkhXIgKgzy05v;7>cC|_0|Uc<uzivY0chs1FhTqe>*Ij*fG{@`
z$ZwKb7@iY{Y5+B8K?=c@4g&*&%;ZKzaYp^g4;95HD=4x|Zc!AG3}=Ek9zB4PL8_E$
zp(Z_mra*A@#=yYP4;JTUh+vt_q;!s}G&eV~s4~8|BysW$rNYV4${v#sE7wd8R;iOH
UM#z<97MEluGsI`+WtK1i01<KTF8}}l

literal 9431
zcmb<-^>JfjWMqH=CI&kO5bpr116T+`GB7L<26MrL1A_$v2ZIBH9D_6i8v_FaD+2=q
zOq~Oi1)~pu3}awmfYBUa6Brnn85kH?Kx`%u!N9=41T_Il%Rq%;G}IWd+d%d~rBQ4K
ziOt~$kqitlnt?$9ECA9Eaw~%Xn9IN*kN{yo=>(`dGN3d}ABYRmCjr$b0aXv9FMteY
zU|@jJF#m!42EsBR0~i<>WMKY+(k>9&8DKO>El4QfX-Nu*?ZOA*F)+CBLBa(_+dzaF
zU^GZANGR}WNealFAT}`==1%>b%p@}t{hSotoXot^3f&3|GhH(iz2ba5BXInHybbcA
zyI&~SSde!>_JBeaWVZ+dBrZYn0_PrWT-$RXC%@_DMlD%^1Lwa!c?r@B!(ao5B-|Jn
z7=*B>6~`eih(p{2hqw+7ac>4#%%fT^heJI#4slf+;tV;NNy(t#HDgH2FG$T}NGm8U
zNoGh(&dD!MWynp;$;nS<NGmVOEJ<ZZD=5m$D@kLBj|XXrFD^+eDv8fc%*<m*D@sje
zh>uUt&Ci1h$1}M5cse=98|fM8nS!zt0}~W7z-f?=7#Tp20iF#?r7}4|DxyG=3=9k}
z1R?nX6vrS86=vWNf{3#-Fo5$QOk4m-9F#U-;u1*Wpu7$fS3nX6r5~8M29h{({x(1o
zM|O_|k~qjsFmoJ`#9?V1B<F!7&IuBLVvpuG93I`QcX<>TJX#Nwu)fG<U|{fQKEiPr
zqTs*jB_0Ka|EgDb6d3sB9T@(rg7_IA$(IlQ|NsAA^%Rc+Lk1{`zq|nE9|G|~LGkhc
zn7<3e2L-{)4PgEz5Fh0AmkYrBRUkgd%P%K@`HMh&kXK(efcdjPe2^Dk7J&JaKzxwb
zUMBqi|KISfXXnK@kIqLP%}+iAgm@l5uz}Zs!Q=Qr2QYP9lh=XaMfv~#|G%V70J(oh
zx=;c`sE2NXPy&NT^AA3c&Zp%}v4?knJOm1WFR}YUQvV<D%Qt|;c7S*=^#^}2dvrea
z=se_c@Rhv>;{}h4KT4Tm5AOmg^k}_Q@+H=z^B6?`BdGoc+n?MD3=E}$p&p%2U)25o
z|3CH^>oaZz2FA~aLGJNrw(a0lU|=YbgF8VabiaT=0)t2A)BhJ@UsnJB|DRu;;bk#|
z1{o9V(fp>uqqBB}M`!5*kIv8;oxW3iI$b+_I(-{Fx?LCSI|1^0=kXW&|NZ~((dl}@
zqdRnkN4M*ReHS1SV5?pf{`>#`xa%8`<6nn#yT0+5;n6G73gUWnyGnR4Ug&0MJy2qX
zRYsuO^-XU9BgoLsi=ChLvGX}Fcyu0z7%cJc|Nmy!H;mn`Z>(LX)GPt1^5{Ic3sj!%
z1ErMxLVOMkFANd(zpP?qU_iI10qO+)z6hMM9y2@+K4A9fJpSVFU$9?K9Cv*I^1{nW
z3<wjCf<!#JT~B!Q$}nNK25b-$$e=2)LCnWppMZ3|)`O`>7{YkLqnn{S^o2+AC69xT
zm^>IS90w_7Xm))9@-s>}LOfyZ`lNn)v+E0-(l7Y?{r~^}-|adD;gFAiz&<(v_0enV
zZr3UN>s{Y;9_n^|!FU118Ub)vp$HxDU<4`bWWsB2{VJ&OFF^?b>gkskz<~#{?FGoT
zCup|K0Nd95gOR^=5m@C5bek`9Gk}FK0__PX(0X~Mf&$E=JM@4D<Ed^TkQgi=LFuF0
z^$A!vEHNGc+XoiHFc2IIgbakl%oA(ZDf~Su|3Rq}6csPPe2)L1Y_JF9pyRF&AQU(W
zy$*#2$7@|;(m*!@qw52RMW|T`kpy<&kBSK}YheC{=xlz$$auio^+74;ao0N__q{g8
zu5cG9>@*K`Ug#BJ=yrX;c%k#q!B_H(7cTzjbiIRBi$Hhjo^IbepbYBR`C*?baxM<~
z1qv5%NW6?>1Z5^rdLTCSb%3Hw0A%n!Q)E;B`~;Z_&RQ?eBc&aTcm*dENTTRw04KFM
z3=9m-wNDsIeULH&deX#_5s;i>h3u5Vpa1{A_9ELXJ7lv|e}WPfp6t)xFZBQae^7XU
zst7JnOXK1XkRDJ%#GYPW@b`Csi%4)ByjJdZg=90dpaPXyX@rcbKLjrRI*-4w{sDHz
z0gyI*ur^ph&EHe?3oLg5B&Yo2|9^0C?$LatB0BaksN(6Y{o~PD`T;B)dzfF|g@M5_
z+_Ce&W0+^>H^)%V&Z~|gKE0|dm=qX7J$iLPbg+-+Q^VUH2mdmIS{g5UzyJU5(fYPT
zj$gh7l$05KI$!y8esk&k@AzNjjz{xhMi0x^rLRB@kmy*)7{^%0ILG+IX=!>e^`8%?
zO#qn$QUPkBFfxFe8Vn4HVAB{F+&q0;ZB>i8f<qibJcC0#ovpZP6ciNl6ciLdv>%A(
z(nSXh3aaLrV9dpUE~m@Ipa7DxRW;1CRn;}m<YG`Y&$I;*3_7ajnG6h!(#+Nt3=H7D
z;GJ*(|ATy1^Zoz-5C#T@oFD)HuVG+dc>3@E|0@g(41fOp|NjLf{{R1f8%72OQ2zns
zbjGS62F3~jMrj^)jtPwH0w8fv-%sKD|NmMb1yBqU2Qg9AGcc%t`dADM3<2N%|6c%7
zz$f6wC*j4<UCz<KU@v8@WvrqE5(cRQ^|^L@|Np-d)XH+?6KG>{=4EqjW@gT3<6#H2
zte|G)FfcH1{`>zQRBtoD)cHb`fXwIt8Tjx2|5Om|#3#_t<jN<}$Lz|d(97b=r_sae
z!DrCM=D}yt%<jwQki%!;$Y<clr{Tn>;KV24#3$gy#{u>~C|orD{r?Y%NswXm!wqbZ
z{$&D`hS>}%N<nNG%>&|qFmzmi7sO^@VA#S7(Ju+*!}`_IP(Cy`8Khu*X#W`$w;&;y
z{MUc~^FjOsX#W`0AO`UZm?7qaq6)-sfXahP5D-5NL@+Qgfcjw|rUZy!U|<l1(l83t
z00*%@K>Ihau)YBjWI&`tP(vFe1nOskXbl#K{V;I>C?Do-n0>H51&sgiKg2&Dpzi+<
z<-@}B1C)OaDi5-Y9vC)mfo_)>G+y1HbQqLQgVJSCx(!NCgVM{O^foAsZmzqtvz3B|
zyI-iLf{C7mo`HgqfuXU1sezG#MnO?(QEE<VVsWY_1DY~J1_r(2%G{E~BnG|Wk|GG5
z0b^z6m82FGFzDswm!#@BI(h1rB&I`o>3OAkNu`-NDY}^{40<45Mq+UWgI-EyUU6kE
zgf1z9$dsiP6=&w>p>X1h81#x#a}q%spsa$N5(d4@yvz~?y^_?55(d4r%)Cqny`ub-
z#F9h?z0{2Ow4%h^)cA~)A_NbjBR;XHD6tZv1IEruE&@9vF*h@rK`%YO1Wf3GZGo7W
zR9wuUmz<xQo0?aG9-JieLHQ7r=3wa^G0p}t3YHHcJk-8DNKBc5fdQ2EVfqo{SRh%<
zu@X?n79s;74bk-5K=ZE+)Id;q3{eRo(e;DM*>C^<=R=h<z{&^Mcpc3Bu<`<yFCnHe
zF!({s2AF<W`OyFs_yAQ1Dqle11k(@8$6?U%EvRw^Sb0<cEl^<P57Zcxu@R7190LOb
zDDT7c!}9wzsD4;|26I159gGGIM}W!_7#~)ie1Pf)jl+TLL3ck$41}wo<E$|Ku<}a)
z8qhF*!t6))Kd3GTDTV2Wm4_N2K?Vi}SosOk52_>3_4h%?i$LxJxrYH(PJlEpFu>~t
zupHE_FgA$pMANST)vo~653cJ$3ZWP@UI}8t=qYIWVdL5dp!&an5;sT#6vOmExeT+=
z^uy}AEzk-OR_?>pLxTmzWdOMkl&)amcK}*n9)Rj^hO2}yVBrVhF)*xv+7HtY8wZAs
z55vZRVd~NC--f2Y0Xm-C03BbJ0XYt)nSp@;-Tng%3=ANp(By_RZVu6et{$C#hJk?r
zG$aqx&jGFfIiLZJt{$BaN?#x+!t{eihCnt85UKwzv|S3*53B!R<N2_0dszB}*$*|I
z;TbgiVESR>qc1=M&RD_=rXNPXM$->#e|&)Ihm5L%d<4NTeJ~o-Rz~)J0CXHP0Nf}=
zjN?LdW5(Yvs2gGK2My<e)Fwdn!|VZxfiQad`44r!0Z0-F!}NpN{Yb*_@(<)@kT}Sj
zFuOr~7@h&MAEpp0fyM=!2TS`<C16ScRIxEIFnoYUoGb$agFe)+AO)~;4U|8z*)IU@
z4KOeqM>8DVevnzXFf+(?(E1En4Ko81e7*-&oS6Yu|DlSr!0R<saaIOceT6E{2Cs)u
z#n~BP^$V&v2Lr6$Ko#d?fYk@6;#>@{avoKjn*m;aBS|nb@G!v2Wh7w+W(HmcSb2*o
z&Ij)|po;S|z{*EdaRCNcxrZt)$N($PP{oB1<rI<xGlMV#to%U|W?*IzVStq@sN$jw
zu<`;`Tns)>f+{Y~0L$N~;t~w7e2XeB$pFibsN%5R5<(P|E?F2pg0dq*5X!<P4$9Zq
z#6i=zj9B~&YDX|J@G^XW)~`^b86a~>ybKCx{X-A%91|}CdcP6W-eF`w_zUE2kXgCl
zas)H{LG?4pURXaAq!xrB^(JOK?gzIKFv9_te?j#RBZDBr41P#I4rDe6-vYY_d%R-{
zH&8tQnvq0LKi|RTU`Z!nk7K%Dm=V;*#!M%gj3Co_89=i)AoszFB}=e4NCXvoGh(*A
zK`PS0;+W~Jh!Ip*NMNM5N*w0*g2gfO>141tFT)0CdH|21FfcGI0E>fMfQmQc5WfU<
z4|=}1k3;<zu)k2v2Jx7g7#IYZBp4!~^&&_Ngatt242bdsR<Fr3LBbg{JqA(-9=ie6
z9bj>g2r71F0)?9-Lk61p9$@vD<(ofLJ$k;5hN=h6$AN5sR(}kjehiEa5{KbpChYOs
z4mJm-3PH~X`IiYZewX7g=OEY|Eaf?<r>mFDkX%$!TvD2rrkBhRAMfQH6z}6191<VT
zkdj)Inx0u)l3EmBk{h3#lb@Gb%z#BGK0YNsK0PNtDKRHLr6j+oI6kqof+0CSw;(6A
zBsE3P%-GbD0jr{TQ2!dzTaHi7D=DgEfb_3ZN^^578RDV6Z7?q$<X9A?#U+U)@!)B5
zu*2Pg9DQBmUHx3*<DuSwu~BrUFvPn=`Z@Y~Iy1y0n+cwAXNY(A@pp3ciT8JN3w8~O
z4{>zzaRnI%o|I26E(JRnJSFeofIJBw?-}A7j|i_2UuRHgg(N2Bq%wed`jDA>hWJzv
z2f6|P%>+<TFo3<2nU@})n4ZG`UM_&93VnG1x-`gqsAbS41?VbJ781nArxfSMXC&sO
z<fJmhd-}&0m!xFo#g`VRrofgWfP4)KLFj@7G~-h9QW&82#m6TV7em>Ql?$Ln4`}MW
n%Ho5diOM-AvA8(37`g%iO<jD5uQMz%L17C^b|@hWjz0zfcpf|3

diff --git a/stat_new/stat_print.c b/stat_new/stat_print.c
index df6b145..9a32d37 100644
--- a/stat_new/stat_print.c
+++ b/stat_new/stat_print.c
@@ -2,7 +2,8 @@
  * 
  * Usage:
  *                  printing_agl (input_file, output_file, number_of_molecules, 
- * true_label_molecules, num_of_molecules_in_aglomerates, aglomerates, statistic)
+ *          true_label_molecules, num_of_molecules_in_aglomerates, aglomerates, 
+ *                                                statistic, type_of_aglomerate)
  */
 
 #include <stdio.h>
@@ -11,7 +12,7 @@
 
 int printing_agl (char *input, char *output, const int *connect, int num_mol, 
                   const int *true_label_mol, const int *num_mol_agl, 
-                  const int *agl, const int *stat)
+                  const int *agl, const int *stat, int *type_agl)
 /* input            - name of file with coordinates
  * output           - name of output file
  * connect          - connectivity graph for all molecules
@@ -20,6 +21,7 @@ int printing_agl (char *input, char *output, const int *connect, int num_mol,
  * num_mol_agl      - massive of numbers of molecule in aglomerates
  * agl              - massive of aglomerates
  * stat             - massive of statistics
+ * type_agl         - massive of numbers of aglomerate types
  */
 {
   int i, iso, j, k, p, type, *label_matrix, **matrix;
@@ -28,8 +30,12 @@ int printing_agl (char *input, char *output, const int *connect, int num_mol,
  * type             - number of cycle in aglomerates
  * label_matrix     - massive of indexes of molecule
  * matrix           - connectivity graph
+ * f_out            - output file
  */
   
+  type_agl[0] = 0;
+  type_agl[1] = 0;
+  
   f_out = fopen (output, "a");
   
 //   head
@@ -68,9 +74,15 @@ int printing_agl (char *input, char *output, const int *connect, int num_mol,
         for (k=0; k<num_mol_agl[i]; k++)
           p += matrix[i][j];
       if (p == (2*num_mol_agl[i]-2))
+      {
         type = 0;
+        type_agl[0]++;
+      }
       else
+      {
         type = (p - (2*num_mol_agl[i]-2)) / 2;
+        type_agl[1]++;
+      }
       
 //       printing class of aglomerate
       fprintf (f_out, "  %7i=%1i=%-7i\n", num_mol_agl[i], type, iso);
diff --git a/stat_new/summary_stat.c b/stat_new/summary_stat.c
new file mode 100644
index 0000000..473edaa
--- /dev/null
+++ b/stat_new/summary_stat.c
@@ -0,0 +1,70 @@
+/* Library for summary statistic
+ * Usage:
+ *            summary_statistic (filename, number_of_step, number_of_molecules, 
+ *                                        type_of_aglomerate, summary_statistic)
+ */
+
+#include <stdio.h>
+#include <stdlib.h>
+
+
+int summary_statistic (char *filename, int step, int num_mol, const int *type_agl, const int *stat_all)
+/* filename         - name of output file
+ * step             - number of steps
+ * num_mol          - number of molecules
+ * type_agl         - massive of numbers of aglomerate types
+ * stat_all         - massive of summary statistics
+ */
+{
+  float conc, p, pn, type[2], x, y;
+  int i, index;
+  FILE *f_out;
+/* conc             - concentrate of aglomerates
+ * p                - probability of aglomerates
+ * pn               - weight probability of aglomerates
+ * f_out            - output file
+ */
+  
+  index = 0;
+  for (i=0; i<num_mol; i++)
+    if (stat_all[i] != 0)
+      index = i;
+  
+//   head
+  f_out = fopen (filename, "a");
+  fprintf (f_out, "SUMMARY STATISTIC:\n");
+  fprintf (f_out, "|   n   |   N   |    C    |    p    |    pn    |\n------------------------------------------------\n");
+  
+  for (i=0; i<index; i++)
+  {
+//     calculating concentrates
+    x = stat_all[i];
+    y = step;
+    conc = x / y;
+    
+//     calculating probabilityes
+    x = (i + 1) * stat_all[i];
+    y = step * num_mol;
+    p = x / y;
+    pn = (i + 1) * p;
+    
+    fprintf (f_out, " %7i %7i %9.2f %9.5f %9.5f\n", i+1, stat_all[i], conc, p, pn);
+  }
+  
+//   types of aglomerates
+  x = type_agl[0] + type_agl[1];
+  type[0] = type_agl[0];
+  type[1] = type_agl[1];
+  fprintf (f_out, "------------------------------------------------\n");
+  fprintf (f_out, "LINEAR=%7.5f\nCYCLE=%7.5f\n", type[0]/x, type[1]/x);
+  
+  fclose (f_out);
+  
+  return 0;
+}
+
+
+int main (int argc, char *argv[])
+{
+  return 0;
+}
\ No newline at end of file