function []= plotfunc(y)
%  [[  plotfunc(y)  ]]
%
% 関数y(x)のグラフを描画します．
% 引数:
% y : 文字式で入力されたxの関数
% 
% Copyright (C) 2007 清水 信行, 曽我部 潔, 小金澤 鋼一 All Rights Reserved.

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

disp('図の描画範囲を入力してください')
disp('入力の後にどれかキーを押すと図がダイアログに表示されます')
xmin=input('x min= ')
xmax=input('x max= ')

xx=xmin: (xmax-xmin)/100: xmax;
z=zeros(1,101);
for k=1:101
    x = xx(k);
    yy= eval(y);
    z(k)=yy;
end
plot(xx,z)
grid
title('Plot of the equation to be solved')
xlabel('x')
ylabel('y(x)')