function mainfunction()

a = 1:4;
b = a * a';
disp('input array:')
disp(a)
disp('scalar product a * a''');
disp(b)
endfunction