//http://mathonweb.com/help_ebook/html/algorithms.htm #include #include using namespace std; constexpr double pi=3.14159265358979323846264338327950288419; constexpr double e = 2.7182818284590; double squareroot(double m); double power(double n,int p){ double product=1; for(unsigned i=0;itwopi){ n-=twopi; } int sign=1; if(n>pi){ n-=pi; sign=-1; } double s=n-power(n,3)/6+power(n,5)/120-power(n,7)/5040+power(n,9)/362880-power(n,11)/39916800+power(n,13)/6227020800-power(n,15)/1307674368000+power(n,17)/355687428096000-power(n,19)/121645100408832000; return s*sign; } double cosinus(double n){ double c = sinus(n+pi/2); return c; } double tangent(double n){ double t=sinus(n)/cosinus(n); return t; } double arctan(double n){ if(n>1) return (pi/2 - arctan(1/n)); if(n>(2-squareroot(3))) { return (pi/6+arctan((squareroot(3)*n-1)/(squareroot(3)+n))); } // cout<0.000001); return guess; } double cubicroot(double m){ double prevguess=m; double guess=prevguess; do{ prevguess=guess; guess=guess-((guess*guess*guess-m)/(m*guess*guess)); }while(abs(prevguess-guess)>0.000001); return guess; } template T abs(T n){ if(n<0) return n*-1; return n; } tuple scientificnotation(double n, bool out=false){ int p=0; if(n>=1){ while(n>=10){ n/=10; p++; } }else{ while(n<1){ n*=10; p--; } } tuple t; t= make_tuple(n,p); if(out) cout< t=scientificnotation(n); double m=get<0>(t); int p=get<1>(t); double x=squareroot(m); double y=(x-1)/(x+1); double lnx=2*(y+y*y*y/3+y*y*y*y*y/5+y*y*y*y*y*y*y/7+power(y,9)/9+power(y,11)/11+power(y,13)/13+power(y,15)/15); double l=2*lnx+p*2.30258509299; return l; } int round(double n){ if(n<0) return int(n-0.5); return int(n+0.5); } double exp(double x){ int n= round(x); double r=x-n; double en=n<0?power(1/e,abs(n)):power(e,n); double er=1+r+r*r/2+r*r*r/6+r*r*r*r/24+r*r*r*r*r/120+r*r*r*r*r*r/(120*6)+r*r*r*r*r*r*r/(120*6*7)+r*r*r*r*r*r*r*r/(120*6*7*8)+r*r*r*r*r*r*r*r*r/(120*6*7*8*9)+power(r,10)/(120*6*7*8*9*10); return en*er; } double pow(double x,double p){ return exp(p*ln(x)); } double arcsinus(double x){ return arctan(x/(squareroot(1-x*x))); } double arccosinus(double x){ return arctan(squareroot(1-x*x)/x); } int main(){ double deg=pi/180; cout<<"\n\nSINUS\n"<

My Math Library