t=linspace(0,50,N);