function [e] = errors(d)
e = d - model(d);