load penddat
t=penddat(:,1);
x1=penddat(:,2);
x2=penddat(:,3);
x3=penddat(:,4);
x4=penddat(:,5);
u=penddat(:,6);
clear penddat

plot(t, [x1 x2 x3 x4]);