function hcurve(pl,pr,level)


A = [0,-1;1,0];
B = A;
if level == 0
   pmid = (pl+pr)/2;
   pup = pmid + A*(pmid-pl);