%%
%% This is file `penddae.m',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% dae.dtx  (with options: `penddae')
%% 
%% (c) Tony Roberts, v1.0b, 16 May 2000, aroberts@usq.edu.au
%% Permission is hereby granted, free of charge, to any person to copy
%% this software and associated documentation files.  But under no
%% circumstances is it to be modified.  This copyright and permission
%% notice shall be included in all copies or substantial portions of the
%% Software.
%% 
%% Used by pendrun.m
function [f,k,m]=pendae(t,w,wd)
x=w(1); y=w(2); u=w(3); v=w(4); s=w(5);
xd=wd(1); yd=wd(2); ud=wd(3); vd=wd(4);
f=[-xd+u
   -yd+v
   -ud-x*s
   -vd-y*s+1
       x^2+y^2-1];
k=[0 0 1 0 0
   0 0 0 1 0
   -s 0 0 0 -x
   0 -s 0 0 -y
   2*x 2*y 0 0 0];
m=-diag([1,1,1,1,0]);
%% 
%%
%% End of file `penddae.m'.