% Octave用   2006/02/14　小川智浩　作成．
% インボリュート曲線を記述する．

clear
format long e

%インボリュート曲線の記述--------------------------------------

%図形に関するパラメータ
disp('involute')
r = 1.5

%描画に関するパラメータ
dt = pi/180;
Waku_Xmin = 2.0;
Waku_Xmax = 2.0;
Waku_Ymin = 2.5;
Waku_Ymax = 2.5;

%初期化
XX = [];
YY = [];

for t = 0:dt:(2*pi)
	XX = [XX r*(1+t^2)*sin(t-atan(t))];
	YY = [YY r*(1+t^2)*cos(t-atan(t))-r];
end

%-------------------------------------------------------

%図形描画
	axis([min(XX)-Waku_Xmin max(XX)+Waku_Xmax min(YY)-Waku_Ymin max(YY)+Waku_Ymax]);
	
	plot(XX,YY,'-;involute;')
	hold('on')
%　X軸の追記
        plot([(min(XX)-Waku_Xmin*0.99) (max(XX)+Waku_Xmax*0.99)],[0.0 0.0],'g-;;')

%　Y軸の追記
	plot([0 0],[(min(YY)-Waku_Ymin*0.99) (max(YY)+Waku_Ymax*0.99)],'g-;;')
	hold('off')

gset terminal windows

