			#!/usr/bin/python

			import sys

			#format output = (r- = (a, b), r+ = (a, b), tmedio, terror)
			def retta(yi, yim, dyi, dyim, xi, xim):

			rho = yi - yim
			eps = dyi + dyim
			dx = xi - xim

			out = [[0, 0], [0, 0], 0, 0]

			out[0][0] = (rho - eps) / dx
			out[0][1] = yi - dyi

			out[1][0] = (rho + eps) / dx
			out[1][1] = yi + dyi

			t = 0

			for i in range(2):
			t -= out[i][1] / out[i][0]

			out[2] = t / 2 + xi

			err = 0

			err -= out[1][1] / out[1][0]
			err += out[0][1] / out[0][0]

			out[3] = err

			return out


			def main(argv):

			if len(argv) < 6:
			print "More arguments bitch"
			print "Format: <y_i> <y_{i-1}> <err y_i> <err y_{i-1}> <x_i> <x_{i-1}>"
			return

			args = []

			for i in range(6):
			args.append(float(argv[i]))

			out = retta(args[0], args[1], args[2], args[3], args[4], args[5])

			print "r-: ", out[0][0], "(x -", args[4], ") + ", out[0][1]
			print "r+: ", out[1][0], "(x -", args[4], ") + ", out[1][1]

			print "h critique:", out[2], "+-", out[3]


			if __name__ == "__main__":
			main(sys.argv[1:])

rSTICAZZI/labo/H8/retta.py
166e407a59a5master

retta.py
View Options