matl

 a=(1+sqrt(t))/2;
b=(1-sqrt(t))/2;
c=1/sqrt(5);
x=1:20;
y=1:20;
x=1/c*(a.^y-b.^y);