% function mg1 for Mackey-Glass function [mgx] = mg1(x,xd) mgx = -0.1*x + .2*xd/(1+xd^10);