Sinus

#include <stdio.h>

#include <math.h>

#define PI 3.141592653589793

double rest(double x)

{

double x_mod;

x_mod=x;

while (x_mod>=2*PI)

x_mod=x_mod-2*PI;

return x_mod;

}

double fact(unsigned int x)

{

if (x==0)

return 1;

else

return (double) x*fact(x-1);

}

double sinus(double x, unsigned char rank)

{

double val=0;

int i;

for (i=0;i<rank;i++)

{

if (i%2==0)

val=val+pow(x,2*i+1)/fact(2*i+1);

else

val=val-pow(x,2*i+1)/fact(2*i+1);

}

return val;

}

void gen_vect(float *x,float l_lim,float r_lim, int nrp)

{

double step;

int i;

step=(r_lim-l_lim)/nrp;

x[0]=l_lim;

for (i=0;i<nrp-1;i++)

x[i+1]=x[i]+step;

}

float get_max(float *x, int nr_elem)

{

float maxim=x[0];

int i;

for (i=1;i<nr_elem;i++)

if (maxim<x[i]) maxim=x[i];

return maxim;

}

float get_min(float *x, int nr_elem)

{

float minim=x[0];

int i;

for (i=1;i<nr_elem;i++)

if (minim>x[i]) minim=x[i];

return minim;

}

void swap(float *x, float *y)

{

float temp;

temp=*x;

*x=*y;

*y=temp;

}

void ord_up(float *x,int nr_elem)

{

int i,j;

for (i=0;i<nr_elem-1;i++)

for (j=i+1;j<nr_elem;j++)

if(x[i]>x[j]) swap(&x[i],&x[j]);

}