Text-only Table of Contents (frame/ no frame)
(24) Matlab Subfunction and Function Handle Example Previous Top Next

  Subfunction and Function Handle Example


function ode_ex2
tstart =0;
tfinal =20;
trange=[tstart tfinal];
y0=[2;0];  % inital conds y(0)=2, dy/dt((0)=0
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
[t,y]=ode45(@vdpol,trange,y0,options);
plot(t,y(:,1),t,y(:,2),'r');
legend('y','dy/dt');

% define differential equation function
% solve van der Pol equation
% d^2y/dt^2 = mu*(1-y^2)*dy/dt-y

function dy=vdpol(t,y)
mu=2;
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=mu*(1-y(1)^2)*y(2)-y(1);

try:
>> ode_ex2


Previous Top Next


Sub_func_example.src  last modified Oct 17, 2007 Introduction Table of Contents
(frame/no frame)
Printable
(single file)
© Dartmouth College