I'm highly suspicious of
rp=[0 4 2 -1 0; 0 1 -4 -3 0;0 0 0 0 0] %x-y values
If you want to do translation you don't need (x, y): you need (x, y, w), where w is the reciprocal scaling factor. So all your points are at infinity, and translating them isn't going to do much.
rp=[0 4 2 -1 0; 0 1 -4 -3 0;1 1 1 1 1] %(x,y,w) values
should get you somewhere.