Rick,
In your code:
delta_U = tnonct(x_c, DF, alpha);
is delta_U .

But in your picture , it is delta_R ,is that a problem ?