bakugu = [,];
yy = [];
lowerbound =  1;
gap = 10000;
for n = lowerbound:lowerbound+gap
    bb = 0;
    Ndigits = dec2base(n,10) - '0';
    f = factor(n) ;
    t = [];
    for r = 1:length(f)
        Rdigits = dec2base(f(r),10) - '0';
        s = length(Rdigits);
        for q = 1:s
            t = [t,Rdigits(q)];
        end
    end
    if length(Ndigits) == length(t) && length(f) ~= 1
        lo = log10(n);
        Y = ceil(log10(n))+isinteger(lo);
        lg = length(Ndigits);
        for w = 1:lg
            for m = 1:lg
                if Ndigits(w) == t(m)
                    bb = bb+1;
                end
            end
        end
    else
        bb = 1;
    end
    if bb == 0
        validateattributes(t, {'numeric'}, {'integer', 'nonnegative', '<', 10});
        y = polyval(t, 10);
        bakugu = [bakugu;[n,y]];
    end
end