phi=[1:1:54]; for j=1 end