pdbeta1_8;
eta8slice = eta1(:,40);
pdbeta1_7;
eta7slice = eta1(:,40);
pdbeta1_6;
eta6slice = eta1(:,40);
pdbeta1_5;
eta5slice = eta1(:,40);
pdbeta1_4;
eta4slice = eta1(:,40);
pdbeta1_3;
eta3slice = eta1(:,40);
pdbeta1_2;
eta2slice = eta1(:,40);
pdbeta1_1;
eta1slice = eta1(:,40);

figure
hold on
plot(x, eta1slice, 'r-.')
plot(x, eta2slice, 'y-.')
plot(x, eta3slice, 'g')
plot(x, eta4slice, 'c')
plot(x, eta5slice, 'b')
plot(x, eta6slice, 'm')
plot(x, eta7slice, 'r')
plot(x, eta8slice, 'y')
hold off
