#include <stdio.h>
#include <math.h>
#include<stdlib.h>

double precyzja[] ={0.01,0.0001,0.000001};
FILE *fp;


int liczba_wyrazow(double x, double precyzja){
    int i = 1;
    double suma = x;
    while(1){
        suma *=  x*x /(2*i*(2*i+1));
        if(suma < precyzja) break;

        i++;
    }
    return i - 1;
}

double silnia(int s){
    if(s==0 || s==1) return 1;
    else return(s*silnia(s-1));
}

double potega(double x, int n ){
    double val = 1.0;
    int i;
    for(i=1; i<=n ;i++){
        val *= x;
    }
    return val;
}

double sinus(double x, int n){
    int i;
    double val = x;

    for(i=1; i<=n; i++){
        double p = potega(-1.0,i);
        double px = potega(x,2*i+1);
        double s=silnia(2*i+1);

        val += (p * px) / s;
    }
}

double sinus_tyl(double x, int n){
    int i;
    double val = x;

    for(i=n; i>=1; i--){
        double p = potega(-1.0,i);
        double px = potega(x,2*i+1);
        double s=silnia(2*i+1);

        val += (p * px) / s;
    }
    return val;
}

double sinus_przod_z_zapamietaniem(double x, int n){
    int i;
    double val = x;
    double wyraz = x;

    for(i=1; i<=n; i++){

        wyraz *=  x*x /(2*i*(2*i+1));

        if(i % 2 == 0)    val += wyraz;
        else val -= wyraz;


    }
    return val;
}

double sinus_tyl_z_zapamietaniem(double x, int n){
    int i;
    double val = x;
    double wyraz = potega(x,2*n + 1)/silnia(2*n+1);

    for(i=n; i>=1; i--){


        if(i % 2 == 0)    val += wyraz;
        else val -= wyraz;

        wyraz /=  x*x /(2*i*(2*i+1));
    }
    return val;
}

void print(double x, double val, double sin,double prec){
    fprintf(fp,"x: %lf wynik: %.20f błąd: %lf, precyzja:%.20lf\n",x,val,fabs(sin-val),prec);
}

void wypisz_wyniki(double x) {
    int i;
    int n;
    double val;
    fprintf(fp, "x: %lf wynik: %.20f\n", x, sin(x));
    for (i = 0; i < 3; i++) {
        n = liczba_wyrazow(x, precyzja[i]);
        print(x,sinus(x,n),sin(x),precyzja[i]);
        print(x,sinus_tyl(x,n),sin(x),precyzja[i]);
        print(x,sinus_przod_z_zapamietaniem(x,n),sin(x),precyzja[i]);
        print(x,sinus_tyl_z_zapamietaniem(x,n),sin(x),precyzja[i]);
    }
}


int main(){

    fp = fopen("wyniki.txt","w");
    if (fp == NULL) {
        fprintf(fp, "Nie mogę otworzyć pliku test.txt do zapisu!\n");
        return 1;
    }
    wypisz_wyniki(3);
    return 0;
}