function nfa = calNFA_Bars(N, r, m,n)
%% nfa = calNFA_Bars(N, r, m,n)
% For input output description follow GG.m script and comments
%
% Written by Boshra Rajaei 3/11/16
% b.rajaei@sadjad.ac.ir
%
% Copyright (c) 2016-2017 boshra rajaei <b.rajaei@sadjad.ac.ir>
%
% This program is a free software: you can redistribute it and/or modify
% it under the terms of the GNU Affero General Public License as
% published by the Free Software Foundation, either version 3 of the
% License, or (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU Affero General Public License for more details.
%
% You should have received a copy of the GNU Affero General Public License
% along with this program. If not, see <http://www.gnu.org/licenses/>.
%
%%
t = 5;
t = degtorad(t);
c2 = 3;
Ntest = c2*N*(N-1)/2;
PH0 = (t*r^2/(m*n))^2;
nfa = Ntest * PH0;