mgequation.m

一个时滞系统的工具箱

function output = mgequation(t, x, tt, xt, initfun) 
tmp = ydelay(-17, t, tt, xt, initfun);

output  = -0.1*x + 0.2*tmp/(1 + tmp^10);