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

//vypíšeš při čtení parametry amplituda, frekvence, doba trvání. gcc -o ukol09 ukol09.c -lm   ./ukol09 2.5 200 500
//program vyplivne vzorky sinusovky

int main(int x, char *y[]){

  int z;
  float amplituda, frekvence, dobatrvani, cas, vystup, vzorky, sinth, costh, yted, yzpoz, ydvakratzpoz, DSP;	//definujeme proměnné typu float

    amplituda = atof(y[1]);
    frekvence = atof(y[2]);
    dobatrvani = atof(y[3]);
    vzorky = dobatrvani/0.125;
    if(frekvence <= 0 || amplituda <= 0 || dobatrvani <= 0){
	printf("Žádný zápor!\n");	//když je podmínka splněna, program pokračuje dál, když však není splněna, program skončí
	return 1;
    }

    sinth = sin(2 * M_PI * (frekvence/8000));	//výpočet sin theta
    costh = cos(2 * M_PI * (frekvence/8000));	//výpočet cos theta
    

    yted = sinth * amplituda;	//počáteční hodnota sin theta
    yzpoz = 0;
    ydvakratzpoz = 0;

    for(z=0;z<vzorky;z++){	//cyklus vypíše výsledný čas, sinus a sinus pomocí DSP algoritmu (vzorky)
	vystup = amplituda * sin(2 * M_PI * frekvence * 0.000125 * z);
	cas=0.000125 * z;
	DSP = yted;
	ydvakratzpoz = yzpoz;
	yzpoz = yted;
	yted = 2 * costh * yzpoz - ydvakratzpoz;

	printf("%f %f %f\n",cas,vystup,DSP);
    }
}