

m = input('Enter 1 for Dirichlet condition, 2 for Neumann ')

x = 0:.05 :10;

if m == 1

Reflect =  moviein(40);
   for j = 1:40
      j
      t = (j-1)*.25;
      u = bigg(x+t) - bigg(t-x);
      plot(x,u)
      axis([0 10 -2 2])
      Reflect(:,j) = getframe;
   end

else 

Reflect =  moviein(40);
   for j = 1:40
      j
      t = (j-1)*.25;
      u = bigg(x+t) + bigg(t-x);
      plot(x,u)
      axis([0 10 -2 2])
      Reflect(:, j) = getframe;
   end
end
     



