(* file geomanal/geomplot.math  -- Make a plot rmsd vs residue     *)

Off[General::spell1];
$DefaultFont={"Times-Roman",13};

                                     (* Change the name of the input file. *)
 (*===>*)    data=ReadList["geomplot.list",{Number,Number}];

     column1 = Transpose[data][[1]]; 
     column2 = Transpose[data][[2]]; 
     max=Max[column2];
     number=Dimensions[column1][[1]]; 
     ListPlot[column2,  
      AxesLabel->{"Seq. Res. No.","Geometry"}, 
      Ticks->{Range[0, number , 5],Range[0.0, max ,1.]},
      PlotLabel->"GeomPlot", AxesOrigin->{0,0},PlotJoined->True, 
      PlotRange->{{0, number +2 }, {0., max +2 }} ];