display("test")
qlr(2.5,[1,2,3,-1,-2,-3],[1,4,9,1,4,9])

display("equation")

f=@(x)7.960-0.1537*x+0.001076*x^2
display("point 1")

qlr(55,[50,50,50,70,70,70,80,80,80,90,90,90,100,100,100],[3.3,2.8,2.9,2.3,2.6,2.1,2.5,2.9,2.4,3.0,3.1,2.8,3.3,3.5,3])

f(55)
display("pint 2")
qlr(65,[50,50,50,70,70,70,80,80,80,90,90,90,100,100,100],[3.3,2.8,2.9,2.3,2.6,2.1,2.5,2.9,2.4,3.0,3.1,2.8,3.3,3.5,3])

f(65)

display("point 3")
qlr(75,[50,50,50,70,70,70,80,80,80,90,90,90,100,100,100],[3.3,2.8,2.9,2.3,2.6,2.1,2.5,2.9,2.4,3.0,3.1,2.8,3.3,3.5,3])

f(75)

display("point 4")
qlr(85,[50,50,50,70,70,70,80,80,80,90,90,90,100,100,100],[3.3,2.8,2.9,2.3,2.6,2.1,2.5,2.9,2.4,3.0,3.1,2.8,3.3,3.5,3])


display("point 5")
qlr(95,[50,50,50,70,70,70,80,80,80,90,90,90,100,100,100],[3.3,2.8,2.9,2.3,2.6,2.1,2.5,2.9,2.4,3.0,3.1,2.8,3.3,3.5,3])

display("point 6)
qlr(100,[50,50,50,70,70,70,80,80,80,90,90,90,100,100,100],[3.3,2.8,2.9,2.3,2.6,2.1,2.5,2.9,2.4,3.0,3.1,2.8,3.3,3.5,3])

f(100)
