function plotSolutionVideo(sol,y,t)
N=length(sol(:,1,1));
for i = 1:N
    figure(1)  
    u1 = sol(i,:,1);
    u2 = sol(i,:,2);
    u4 = sol(i,:,4);
    u5 = sol(i,:,5);
    u9 = sol(i,:,9);
    u10 = sol(i,:,10);
    u11 = sol(i,:,11);
    u12 = sol(i,:,12);
    u13 = sol(i,:,13);
    u14 = sol(i,:,14);
    K = length(sol(i,:,13));
    eS=sqrt((1+exp(-u14)).^(-1));
    eS2=sqrt(u10.*u10+u11.*u11+u12.*u12);
    sum = 0;
    for a = 1:K
        sum = sum + 0.74*(1+exp(-u13(a)))^(-1);
    end
    %imshow(processo(:,:,1,i))
      %plot(0.74.*(1+exp(-u13)).^(-1),y);
      plot(u4,y);
      xlabel(sum);
      %axis([0 1 0 1])
      %axis([-0.001 0.001 0 1])
      F(i) = getframe(gcf) ;
      drawnow
end
  % create the video writer with 1 fps
  writerObj = VideoWriter('ac100_t1000.avi');
  writerObj.FrameRate = 10;
  % set the seconds per image
% open the video writer
open(writerObj);
% write the frames to the video
for i=1:length(F)
    % convert the image to a frame
    frame = F(i) ;    
    writeVideo(writerObj, frame);
end
% close the writer object
close(writerObj);
end

