The running results in matlab and R language are inconsistent

See more:

MatLab

R

matlab

Yr=Y;
Us{1}=Yr*C;
U=Us{1};
Fxy=Y-Us{1}*C
R
Yr <- Y
Us <- list()
Us[[1]] <- Yr %*% C
U <- Us[[1]] 
Fxy <- Y-(U[[1]] %*% t(C))

What I have tried:

The U, C, and Y values in the two pieces of code are the same except for the number of digits left after the decimal point. In matlab, four digits are kept, and in R, the calculation is accurate. Y is a 246 2 matrix, C is a 2 1 vector, and U is a 246-row list with dimensions that can be multiplied. Why does it work in matlab and R? Different results