clear;
mac = [5.1 5.2 5.3; 4.1 4.2 4.3; 0 0 1];
x = [-1:0.1:1];
y = x.*x.*x;
fun = [x;y;ones(size(x))];

fun = mac*fun;

x = fun(1,:);
y = fun(2,:);
plot(x,y);