<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Dear Kondaiah Samudrala<br>your problem is due to some negative frequency. So change ( tr2_ph = 1.0d-10 ) to some other value i.e -10 to -14 and rerun ph.x and then try QHA. Hope your problem will resolve. <br>Best of Luck<br><br>Muhammad Zafar<br>
PhD Scholar<br>
Department of Physics<br>
The Islamia University of Bahawalpur,Pakistan<br><br></td></tr></table>