Matlab中的非线性函数？

Question

我有如下转发：

Tx：x'= a1x + a2y + a3xy + a4。

Tx：y'= b1x + b2y + b3xy + b4。

function y = transformationTest (I,a1,a2,a3,a4,b1,b2,b3,b4) 
    [M N] = size(I);
    for i = 1: M
        for j = 1: N
            x1 = a1*i + a2*j + a3*i*j + a4;
            y1 = b1*i + b2*j + b3*i*j + b4;
            y(round(x1+1),round(y1+1)) = I(i,j);
        end
    end
end

我的实施是错误的还是什么因为我得到大图像而且毫无意义。