load data.mat
y = f1(data, b);
save result.mat -7 y