On Tue, 27 Mar 2012, Allin Cottrell wrote:
 you can do it in gretl:
 matrix A = {10,55,60,30,80,43,100}'
 matrix A39 = 39*ones(rows(A), cols(A))
 matrix A61 = 61*ones(rows(A), cols(A))
 matrix Atest = (A .> A39) .* (A .< A61)
 print Atest
 scalar n_40_60 = sumc(Atest) 
Sorry, you don't need A39 and A61; this will do it:
  scalar n_40_60 = sumc((A .> 39) .* (A .< 61))
And new in CVS, the boolean && operator works as it ought to:
  scalar n_40_60 = sumc(A .> 39 && A .< 61)
Allin Cottrell