p = plot(x^2, (-0.5, 1.4)) + line([(0,0), (1,1)], color='green')
p += line([(0.5, 0.5), (0.5, 0.5^2)], color='purple')
p += point(((0, 0), (0.5, 0.5), (0.5, 0.5^2), (1, 1)), color='red', pointsize=20)
p += text('(1|1)', (-0.05, 0.1), color='red')
p += text('B', (1.01, 1.1), color='red')
p += text('Test', (0.48, 0.57), color='red')
p += text('D', (0.53, 0.18), color='red')
p.show(xmin=-0.5, xmax=1.4, ymin=0, ymax=2,axes_labels=['$x$','$f(x)$'])