sin(x) = x - x3/3!+x5/5!-x7/7!….
// float sin(float theta) { int i;float sign,x,a; //
x=theta/M_PI;if (x < 0.0) {sign=-1;x=-x;} else sign=+1;
i=(int)(x+0.5);a=x-i;if ((i-i/2*2) != 0) sign=-sign;
return sign*sin2(a*M_PI);
} // float sin2(float x) { int i;float y,r; //
r=x;y=x;
for (i=0;i < 10;i++) {y*=-x*x;r+=1.0/factrl(1+2*(i+1))*y;}
return r;
} // long factrl(int n) { int i;long la;la=1;for (i=2;i ⇐ n;i++) la*=i;return la; } //