% Anime.m
% 
% 解析結果のアニメーション表示
% 
% Copyright (C) 2007 清水 信行, 曽我部 潔 All Rights Reserved.

% =================================================================================
% プログラムの使用条件
% 1. 本プログラムおよびデータを商用で，頒布してはならない．
% 2. 本プログラムおよびデータを著作者の許可なく営利目的で利用してはならない．
% 3. 本プログラムおよびデータの改変は，営利目的でない限り自由である．
% 4. 本プログラムおよびデータを使用することにより生じた損害などについては，著者，
%    プロクラム著作者，(社)日本機械学会，および(株)コロナ社は一切の責任を負わない．
% 5. 本プログラムの著作権は，本書編著者ならびにプログラム著作者に帰す．
% =================================================================================

close; % Close the figure window

for i=1:(Te-Ts)/dt+1
	
	plot([0,0],[0,0])
	axis square
	axis([xmin,xmax,ymin,ymax])
	hold on
	
	t=A_data(i,1);
	fprintf('t= %f\n',t)
	
	for j=1:nae
	  k=2+4*(j-1);
	  xL=A_data(i,k); yL=A_data(i,k+1);
		xR=A_data(i,k+2); yR=A_data(i,k+3);
	  plot([xL,xR],[yL,yR],'linewidth',0.5);
	end
	
	drawnow
	hold off
	set(1,'doublebuffer','on')
end
