%
% File dcf.m
% Differential equation for dc.m
%
function xdot=dcf(t,x)
global r k j
om=x(1,1);
cur=(40-k*om)/r;
xdot(1,1)=k*cur/j;
