nxn matrix inversion in verilog code

Hi,

I have a project on writing a n*n matrix inversion by using verilog code. I hope can doing this by using QR decomposition. Can anyone give me some clue on:

1) how can making the code is suitable for n*n matrix by just changing the parameter ?

2) a module that can do floating point division.