f=@(x)x(~mod(sum(dec2base(x,10)'-48).^.5,1)) f([1,4,9,16,25,1111]) f([1431,2,0,22,999999999]) f([22228,4,113125,22345]) f([421337,99,123456789,1133557799])
Standard input is empty
f =
@(x) x (!mod (sum (dec2base (x, 10)' - 48) .^ .5, 1))
ans =
1 4 9 1111
ans =
1431 0 22 999999999
ans =
22228 4 22345
ans = [](1x0)