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

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

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

     column0 = Transpose[data][[1]]; 
     column1 = Transpose[data][[2]]; 
     column2 = Transpose[data][[3]];
     max=Max[column1,column2];
     number=Dimensions[column1][[1]]; 
     ListPlot[column1, 
       Axes->True, 
       AxesLabel->{"Seq. Res. No.","RMSD"}, 
       Ticks->{Range[0, number , 5],Range[0.0, max ,0.2]},
       PlotLabel->"backbone", AxesOrigin->{0,0},PlotJoined->True, 
       PlotRange->{{0, number +2 }, {0., max +0.1 }}];
     ListPlot[column2,  
       Axes->True,
       AxesLabel->{"Seq. Res. No.","RMSD"}, 
       Ticks->{Range[0, number , 5],Range[0.0, max ,0.2]},
       PlotLabel->"side chain", AxesOrigin->{0,0},PlotJoined->True, 
       PlotRange->{{0, number +2 }, {0., max +0.1 }}];