Дан массив P(m), m<=20. Определить количество положительных элементов K и квадрат их суммы S2

P = [1 5 -6 8 -9 7 5 8 6 3 24 -8 0 1 -5 8 12 -89];
x = find(P >= 0);
K = length(x); 
S2 = bsxfun(@power, sum(P(x)), 2);
variables = {'Исходный массив'};
sx = ['' sprintf('%s', variables{:}) sprintf('\n')]; 
sx = [sx sprintf('  %2.2f \n', P)];
variables = {'Количество элементов'};
sr = ['' sprintf('%s', variables{:}) sprintf('\n')]; 
sr = [sr sprintf('  %2d \n', length(P))];
variables = {'Количество положительных элементов'};
sw = ['' sprintf('%s', variables{:}) sprintf('\n')]; 
sw = [sw sprintf('  %2d \n', K)];
variables = {'Квадрат суммы положительных элементов'};
sq = ['' sprintf('%s', variables{:}) sprintf('\n')]; 
sq = [sq sprintf('  %2d \n', S2)];
fid = fopen('forum.txt', 'wt'); 
fprintf(fid, '%s', sx); 
fprintf(fid, '%s', sr);
fprintf(fid, '%s', sw); 
fprintf(fid, '%s', sq); 
fclose(fid);

P =input('введите значения P:  ')

P =input('введите значения P(не больше 20ти значений):  ')

n = input('введите количество элементов массива:  ');
P = rand(n, 1);

P = input('введите значения P(не больше 20ти значений):  ');
disp('Количество элементов равно')
disp(length(P)) % отобразить количество элементов
if length(P) > 20 
    disp('Количество элементов больше 20')
    disp('Введите меньшее количество элементов!')
    break
end

P = input('введите значения P(не больше 20ти значений):  ');
disp('Количество элементов равно')
disp(length(P)) % отобразить количество элементов
if length(P) > 20 
    disp('Количество элементов больше 20')
    error('Введите меньшее количество элементов!')
end