(* file nmr/average_plot.math  -- Make a plot rmsd vs residue   *)

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

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

     column1 = Transpose[data][[1]]; 
     column2 = Transpose[data][[2]]; 
     column3 = Transpose[data][[3]];
     max2=Max[column2];
     max3=Max[column3];
     number=Dimensions[column1][[1]]; 
     ListPlot[column2,
      Axes->True,  
      AxesLabel->{"Seq. Res. No.","rmsd / A"}, 
      Ticks->{Range[0, number , 5],Range[0.0, max2 ,0.2]},
      PlotLabel->"backbone", AxesOrigin->{0,0},PlotJoined->True, 
      PlotRange->{{0, number+ 0.2 }, {0., max2 +0.2 }}];
     ListPlot[column3, 
      Axes->True, 
      AxesLabel->{"Seq. Res. No.","rmsd / A"}, 
      Ticks->{Range[0, number , 5],Range[0.0, max3 ,0.5]},
      PlotLabel->"side chain", AxesOrigin->{0,0},PlotJoined->True, 
      PlotRange->{{0, number +0.2 }, {0., max3 +0.2 }}];