#! /usr/bin/dc

10k

3 v sy
[La Lb Lc Ld Le Lf Lg Lu Lt Li Lj sz sz sz sz sz sz sz sz sz sz sz q
] sQ 

#a=x1 b=y1 c=x2 d=y2
[Sd Sc Sb Sa
lc la ld lb - ly * + - 2 / la +
Ld lb Lc La - ly * - - 2 / Lb +] sN #returns x y

#a=x1 b=y1 c=x2 d=y2 e=xplus f=yplus g=depth
[Sg Sd Sc Sb Sa

lc la - 2 / Se
ld lb - 2 / Sf

la		#x1
lb		#y1
la le +		#x1+xplus
lb lf + 	#y1+yplus
lN x 
Su		#u=y3 
St 		#t=x3

la le +		#x1+xplus
lb lf +		#y1+yplus
lc		#x2
ld		#y2
lN x
Sj		#j=y4
Si		#i=x4


[l ]P lt n [ ]P lu n [ ]P li n [ ]P lj p sz

1 lg <Q

#1
la 		#x1
lb 		#y1
la le + 	#x1+xplus
lb lf + 	#y1+yplus
lg 1 - 		#depth-1
lK x

#3
lt		#x3
lu		#y3
li		#x4
lj		#y4
lg 1 -		#depth-1
lK x

#5
la le +		#x1+xplus
lb lf +		#y1+yplus
lc		#x2
ld		#y2
lg 1 -
lK x

lQ x
] sK

[#PLOT 2]p
[o]p
[e]p
0 0 1 0 8 lK x
[x]p
