Anonymous user
Trigonometric functions: Difference between revisions
m
→{{header|REXX}}: optimized a few statements here and there, mostly used multiplication instead of division (* 0.5 instead of /2) and the like.
(Added zkl) |
m (→{{header|REXX}}: optimized a few statements here and there, mostly used multiplication instead of division (* 0.5 instead of /2) and the like.) |
||
Line 1,992:
say 'Using' showdigs 'decimal digits precision.'; say
do j=-180 to +180 by 15
stuff = right(j,4) 'degrees, rads='show( d2r(j)),
' sin='show(sinD(j)),
Line 2,001:
end /*j*/
say; do k=-1 to +1 by 1/2
say right(k,4) 'radians, degs='show( r2d(k)),
' Acos='show(Acos(k)),
Line 2,018:
cos: procedure; arg x; x=r2r(x); a=abs(x); numeric fuzz min(9,digits()-9)
if a=pi() then return -1; if a=pi()
if a=pi()/3 then return .5; if a=2*pi()/3 then return -.5
return .sinCos(1,1,-1)
Line 2,032:
parse value format(x,2,1,,0) 'E0' with g 'E' _ .; g=g*.5'E'_%2; m.1=p
do j=2 while p>9; m.j=p; p=p%2+1; end /*j*/
do k=j+5 to 0 by -1; numeric digits m.k; g=
numeric digits m.0; return (g/1)i
e: return,
Line 2,042:
exp: procedure; arg x; ix=x%1; if abs(x-ix)>.5 then ix=ix+sign(x); x=x-ix
z=1; _=1; w=z; do j=1; _=_*x/j; z=(z+_)/1; if z==w then leave; w=z; end
if z\==0 then z=
pi: return, /*a bit of overkill, but hey !! */
Line 2,061:
d2r: return r2r(arg(1)*pi()/180) /*convert degrees ──► radians. */
r2d: return d2d((arg(1)*180/pi())) /*convert radians ──► degrees. */
r2r: return arg(1)//(
show: return left(left('',arg(1)>=0)format(arg(1),,showdigs)/1,showdigs)
tellErr: say; say '*** error! ***'; say; say arg(1); say; exit 13
Line 2,067:
AsinErr: call tellErr 'Asin(x), X must be in the range of -1 ──► +1, X='||x
AcosErr: call tellErr 'Acos(x), X must be in the range of -1 ──► +1, X='||x
sqrtErr: call tellErr "sqrt(x), X can't be negative, X="||x
/* ┌───────────────────────────────────────────────────────────────┐
│ Not included here are (among others): │
| |||