double mysin(double x) { double a=x,x1=x,b; int n=1,f=1,s=1; while(myabs(b)>1e-5) { b=b*(n+1)*(n+2); x1=x*x*x1; s=-s; b=x1/f*s; a=a+b; n+=2; } return a; } double myabs(double x) { (x>=0)?x:-x; return (x); }