%
function v=uN(x,t,N)

v=zeros(size(x));

for n=1:N
    cn=4*(2*sin(1/2*n*pi)-sin(n*pi))/n^2/pi^2;
    v=v+cn*cos(n*pi*t/10)*sin(n*pi*x/10);
end