% NewmarkMain.m
% 
% ニューマーク法のメイン部分
% 呼び出される関数はNewmark.m
% 
% Copyright (C) 2007 清水 信行, 曽我部 潔, 藤川 猛 All Rights Reserved.

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

clear all;hold off
    M=1;C=0;KL=1000;            %定数設定
    eq={'500*D1^3' '1500*D1^2'} %非線形復元力設定
%時間積分     
tic; sol=Newmark(M,C,KL,eq,0,[1,0],0.5,0.01);toc;
%結果の表示
subplot(2,1,1);
X=sol(1,:);Y=sol(2,:);plot(X,Y);xlabel('time step');ylabel('Disp'); grid;%plot1
subplot(2,1,2);
X=sol(2,:);Y=sol(4,:);plot(X,Y);xlabel(' Disp.');ylabel('Fint');grid;    %plot2
