# $CVSHeader: pal8e/dat/stat/plot-Hrmsd,v 1.2 2005/02/15 01:10:24 indraneel Rel $

## Run using GNUPlot

set data style lines
set mxtics 5
show mxtics
set mytics 5
show mytics
#set key below box

set xrange [7:90]
#set yrange [0:0.2]
set xlabel "Helix length"
set ylabel "RMSD of Helix C{/Symbol a} about central axis"
set title "Helix RMSD vs. length using our algorithm"

f(x) = m*x + c

m = 0.0118
c = 2.1756

set terminal x11
plot "Hrmsd-orig" smooth unique lt 2 linewidth 4 title \
		"Average RMSD of Unbroken Helix", \
	'Hrmsd-brok' smooth unique lt 3 linewidth 4 title \
		"Average RMSD of Broken Helix", \
	f(x) lt 1 linewidth 4 title "Approximated RMSD of linear helix"

set terminal postscript eps enhanced "Times-Roman" 22
set output "Hrmsd.eps"
replot


