Pascal/Turbo Pascal - Ayuda Metodo Newton Raphsom

Hola a todos, no se si me pueden ayudar... Yo ya tengo hechas las librerias donde estan los procedimientos basicos de matrices: suma, resta, producto, inversa, etc.. Pero necesito la rutina principal para resolver un sistema de ecuaciones no lineales por el metodo de Newton Raphson. Me dieron la formula pero no se si es esa: Pn = P - (F(p)/F'(p))....

O sino, me podrian ayudar para poder hacer un programa que realize operaciones matriciales especificando la operacion en la linea de ordenes (comando). como por ejemplo: A+ B, y solo introduciendo esto, me haga la suma de dos matrices.

Por su comprension y Ayuda,,, gracias....