% Octave用   2005/09/19　小川智浩　作成．
% 懸垂線を記述する．

clear
format long e

%懸垂線の記述--------------------------------------

%図形に関するパラメータ
disp('catenary')
a = 1

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

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

for t = (-pi):dt:(pi)
	XX = [XX t];
	YY = [YY a*cosh(t/a)];
end

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

%図形描画
	axis([min(XX)-Waku_Xmin max(XX)+Waku_Xmax min(YY)-Waku_Ymin max(YY)+Waku_Ymax]);
	
	plot(XX,YY,'-;catenary;')
	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')

%　epsファイルに出力
        gset terminal postscript eps;
        gset output "graph_catenary.eps";

	plot(XX,YY,'-;catenary;')
%　X軸の追記
	hold('on')
	plot([(min(XX)-Waku_Xmin*0.99) (max(XX)+Waku_Xmax*0.99)],[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

