x=rand(3000,1);

X=[ones(size(x)) ,x,  x.^2,  x.^5];

a=5;

b=3.4;

c=-3;

d=20;

func=[d a b c]*X';

fitglm(X,fun);

Estimate          SE                  tStat                 pvalue