%《MATLAB在電子信息課程中的應用》第二章第六節演示程序q26
% M文件及程序調試
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

% 求fibonacci數的程序
R=input('輸入fibonacci數的最大范圍R= '),f=[1 1]; i=1;
while f(i)+f(i+1)<R f(i+2)=f(i)+f(i+1);i=i+1;
end
f,plot(f)

display('輸入fibonacci數的最大范圍R= '),f=[1 1]; i=1;
keyboard
while f(i)+f(i+1)<R f(i+2)=f(i)+f(i+1);i=i+1;
end
f,plot(f)

% 求素數的程序
disp('求素數(prime number)的程序'),pause
N=input('輸入素數的最大范圍N= '), x=2:N;
for u=2:sqrt(N)
   n=find(rem(x,u)==0 & x~=u);
   x(n)=[];
end,x

%子程序是不能直接執行的,必須由主程序調用,因此下面的子程序只能閱讀,不能執行
%function y = humps(x)
%HUMPS	A function used by QUADDEMO, ZERODEMO and FPLOTDEMO.
%	HUMPS(X) is a function with strong maxima near x = .3 and x = .9.
%	See QUADDEMO, ZERODEMO and FPLOTDEMO.
%	Copyright (c) 1984-94 by The MathWorks, Inc.

%y = 1 ./ ((x-.3).^2 + .01) + 1 ./ ((x-.9).^2 + .04) - 6;