v = 1:16
v1 = 1:4
v2 = 5:8
v3 = 9:12
v4 = 13:16
A = [v1; v2; v3; v4]
B = [v1', v2', v3', v4']
C = A - B
D = C'
E = D;
E(3, :) = 1
E(4, :) = [];
E(:, 4) = []
f = @(x) (x - 1).^2;
g = @(x) x.*sin(x);
F = f(E)
G = g(E)
figure
subplot(1, 2, 1)
fplot(f)
grid on
subplot(1, 2, 2)
fplot(g)
grid on