double COMPUTE_PI(char *seq,unsigned long seq_length, int charge_increment);