pell.c