%《MATLAB及在電子信息課程中的應用》第四章第四節調用的子程序
% 范德堡方方程的子程序
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
% 用ode函數解微分方程所需的被積函數
function  yprime = vdpl (x,y)
global r			% r值由主程序通過全局變量傳送
yprime = [ y(2);  r*(1 - y(1).^2).*y(2) - y(1) ];	% 兩行單列向量