Gradient based and least squares based iterative algorithms for matrix equations AXB + CXTD = F