function dxdt = complexbloch1(t,x)

NN=150;
w=2;

for n=1:NN,
    V(n)=0.5*(cos(n/2))^2+0.01*n;
end

% for n=NN/2-1:NN/2+1,
%     V(n)=1;
% end


dxdt(1)=(w*i*(x(2)-2*x(1))-i*V(1)*x(1)-x(1))';
for I=2:NN-1,
dxdt(I)=(w*i*(x(I+1)+x(I-1)-2*x(I))-i*V(I)*x(I))';
end
dxdt(NN)=(w*i*(x(NN-1)-2*x(NN))-i*V(NN)*x(NN)-x(NN))';

dxdt = dxdt';