How do I get Julia to treat all elements of matrices as "floats" (like C type float)

[una-dinosauria]

You can pass the type at the beginning of the constructor in zeros

...
A = zeros(Float32,i,j)	      # Matrix initialization 
X =  zeros(Float32,i,j)
Y =  zeros(Float32,i,j)
...