m_2=trapz(fr,Y_m2);%moment m2
T_02=sqrt(m_0/m_2);
Ym_m1=(fr.^(-1)).*S;
m_10=trapz(fr,Ym_m1);
T_E=m_10/m_0;
neta=((m_0*m_2)/m_1^2)-1; %largeur du spectre netacarré
Y_m4=(fr.^(4)).*S;
m_4=trapz(fr,Y_m4);
epsilone=1-(m_2^2/(m_0*m_4));
% on suppose qu'on est en eau profonde
Hs=Hm;
T_m=T_02;
T_ptheo=T_m/0.87;
Ts=0.90*T_ptheo;
ti1=sprintf('U=%d \t, ', U);
ti2 = sprintf('F=%d \t', F);
ti3=('Spectre de JONSWAP , ')
title([ti1,ti2,ti3]);
Tp=1/f_pic;
fprintf('la periode de pic \n Tp=%d ', Tp);
fprintf('la periode theorique \n T_ptheo=%d ', T_ptheo);
fprintf('The default values: \n f_p=%d \n alpha=%d \n gamma=%d \n m_0=%d \n
m_1=%d ', f_pic, alpha, gamma, m_0, m_1)
fprintf('\n m_2=%d, \n m_10=%d \n neta=%d \n epsilone=%d \n T_01=%d \n
T_02=%d
\n
T_E=%d
\n
m_4=%d
\n
T_ptheo',m_2,
m_10,neta,epsilone,T_01,T_02,T_E,m_4,T_ptheo)
end