pdtheta8;
u8slice = u1(:,40);
pdtheta7;
u7slice = u1(:,40);
pdtheta6;
u6slice = u1(:,40);
pdtheta5;
u5slice = u1(:,40);
pdtheta4;
u4slice = u1(:,40);
pdtheta3;
u3slice = u1(:,40);
pdtheta2;
u2slice = u1(:,40);
pdtheta1;
u1slice = u1(:,40);

figure
hold on
plot(x, u1slice, 'r-.')
plot(x, u2slice, 'y-.')
plot(x, u3slice, 'g')
plot(x, u4slice, 'c')
plot(x, u5slice, 'b')
plot(x, u6slice, 'm')
plot(x, u8slice, 'y')
plot(x, u7slice, 'r')
hold off
