This final project discusses a new modification on Heun’s method based on a contraharmonic mean for solving an initial value problem with high efficiency. The process begins by deriving the formula for the method discussed followed by stability analysis, consistency analysis, local truncation error (LTE). From the numerical results presented it can be seen that this method is more stable …