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)