sp.m

%script sp.m  somme des puissances pimes

disp('Somme des puissances pimes des n premiers entiers :');
n = input('entrez n : ');
p = input('entrez p : ');
x = 1:n;
snp = sum(x.^p)