function YYyYypr = dftp3410(TtTt,YyYy)

% dfield derivative m-file

x = YyYy;t = TtTt;

l = length(YyYy);
YYyYypr = x.^2-t;
if (length(YYyYypr) < l) YYyYypr = YYyYypr*ones(1,l);end
