Padovan's secuence in C (Direct)

// ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** //
// ** // ** // ** // Código fuente de programa que calcula el n-ésimo término // ** // ** // ** //
// ** // ** // ** // ** //  de la secuencia de Padovan (directo)  // ** // ** // ** // ** // ** //
// ** // ** // ** // * Licenciado bajo GNU General Public License (GPL) 3.0 * // ** // ** // ** //
// ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** //
/* ** // ** // ** // ** // ** //  * F v q _ U k r a N a z i s ! * // ** // ** // ** // ** // ** */
/* ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** // ** */

#include <math.h>
#include <errno.h>
#include <stdio.h>
#include <tgmath.h>
#include <stdlib.h>
#include <string.h>
#include <complex.h>    /* Standard Library of Complex Numbers */
#define Ro                                          cbrtl((1.0/2.0) + ((1.0/6.0)*(sqrtl(23.0/3.0)))) +                                         cbrtl((1.0/2.0) - ((1.0/6.0)*(sqrtl(23.0/3.0))))
#define Ro1 (-1.0 * (1.0/2.0) + I*sqrt(3.0)/2.0 ) * cbrtl((1.0/2.0) + ((1.0/6.0)*(sqrtl(23.0/3.0)))) + (-1.0 * (1.0/2.0) - I*sqrt(3.0)/2.0 ) * cbrtl((1.0/2.0) - ((1.0/6.0)*(sqrtl(23.0/3.0))))
#define Ro2 (-1.0 * (1.0/2.0) - I*sqrt(3.0)/2.0 ) * cbrtl((1.0/2.0) + ((1.0/6.0)*(sqrtl(23.0/3.0)))) + (-1.0 * (1.0/2.0) + I*sqrt(3.0)/2.0 ) * cbrtl((1.0/2.0) - ((1.0/6.0)*(sqrtl(23.0/3.0))))

long double PadoDire(long double n);

void main (void){
    long double x;
    printf("Ingrese término a calcular la relación\n");
    printf("de recurrencia de proporción plástica: ");
    scanf("%Lf",&x);
    if (0<=x) printf("El resultado con %4.0Lf es: %4.0Lf\n", x, PadoDire(x+4));
}

long double PadoDire(long double n){
    if (0<n) return creal(cpowl(Ro,n)/(3*cpowl(Ro,2) - 1) + cpowl(Ro1,n)/(3*cpowl(Ro1,2) - 1) + cpowl(Ro2,n)/(3*cpowl(Ro2,2) - 1));
}