pdtheta19;
u19slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta18;
u18slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta17;
u17slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta16;
u16slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta15;
u15slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta14;
u14slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta13;
u13slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta12;
u12slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta11;
u11slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta10;
u10slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta9;
u9slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta8;
u8slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta7;
u7slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta6;
u6slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta5;
u5slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta4;
u4slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta3;
u3slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta2;
u2slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1

pdtheta1;
u1slice = u1(:,40);
clear a1 c1 eta1 f1 pi1 rmac1 sig1 v1


%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

utot = [u1slice u2slice u3slice u4slice u5slice u6slice u7slice u8slice u9slice u10slice u11slice u12slice u13slice u14slice u15slice u16slice u17slice u18slice u19slice]
mesh(t, x, utot)
