What do you expect your values to look like? Are you expecting a vector something like
(35184372088832*pi^2 - 219978174876521)/(70368744177664*pi),
(281474976710656*pi^2 - 1766249117313443)/(562949953421312*pi)
Did you notice that your k does not change in the loop, and so k==n is either always true (if n happens to be 1, which it is not) or else always false (if n happens to be different from 1, which is the case.)
Did you notice that your statements in the loop do not depend upon S?
Did you notice that "plot on" is not set so you are plotting the same thing every time?
Did you notice that k is a scalar, but you are indexing it as if it were a row vector or 2D array ?