Okay, this is my attempt to get the required results that are in the first post:

Code:double sin(float angle) { double guess; for (guess = angle; !(fabs(guess * guess - angle) < EPSILON); guess = ((x * guess * guess * -1) / ((guess * 2) +1)!) { } return guess; }