select, flag=twiss, clear=true;
select, flag=twiss, range=#S/#E, column=name,s,x,y;

TWISS, DELTAP=0.0;
write, table=twiss, file=twiss.prt;

plot, haxis=s, vaxis=betx,bety, hmin= StartPlot, hmax= EndPlot, vmin=0.0, vmax=150.0, colour=100;
plot, haxis=s, vaxis=dx,dy, hmin= StartPlot, hmax= EndPlot, vmin=-1.0, vmax=6.0, colour=100;
plot, haxis=s, vaxis=mux,muy, hmin= StartPlot, hmax= EndPlot, vmin=0.0, vmax=30.0, colour=100;
plot, haxis=s, vaxis=x,y, hmin= StartPlot, hmax= EndPlot, vmin=-0.010, vmax=0.040, colour=100;

return;