#macro suavizar(P, N, i)
    // suaviza a juntaentre o arco i e o arco j =  mod(i+1, N)
    #local j = mod(i+1, N);
    #local vv = (P[j][1] - P[i][2])/2;
    #declare P[i][2] = P[i][3] - vv;
    #declare P[j][1] = P[j][0] + vv;
#end