penta function

C code for solving pentadiagonal linear equations