% Harmonic series ... s3hn.m
n=input('Give number of terms n= ');

a=1:n;
a=1 ./ a;

disp(['Harmonic series of order ' num2str(n) ' is equal to ' num2str(sum(a)) ]);
disp(['ln(n)+ Gamma is ' num2str((log(n)+0.577)) ]);
