clear all
clc
x0=4;
x1=4;
d=dec2bin(hex2dec('56AB0A'))
e=dec2bin(hex2dec('72A4FB'))
x0=x0+((x0*x0)|5)
x1=x1+((x1*x1)|13)
z=x0
for r=0:24
z=(circshift(z,1))+((z+d)>1)
y=x1&z
end
for r=0:24
y=(circshift(y,1))+(circshift(y,-1))+y+e
z=z&y
end