Ktoś poprosił mnie pewien czas temu o napisanie własnej implementacji Cosinusa Hiperbolicznego (http://pl.wikipedia.org/wiki/Funkcje_hiperboliczne) w oparciu o szereg Taylora. Wersja na szybko jako proof of concept:
[cpp]
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
double silnia (double n){
if (n <= 1)
return 1;
else
return silnia(n – 1)*n;
}
double mycosh(double z){
double result =1;
int ilosc_wyr_ciagu=30;
int i;
for (i=1;i<ilosc_wyr_ciagu;i++){
result += pow(z,2*i)/silnia(2*i);
/*
//odkomentowac jesli chcemy zobaczyc jak maleje delta i czemu ilosc_wyr_ciagu=30
double delta;
printf ("iteracja:%d, delta:%.40f\n",i, ( delta=(pow(z,2*(i-1))/silnia(2*(i-1)))-(pow(z,2*i)/silnia(2*i)) )) ;
if(delta == 0)printf("ZERO");
*/
}
return result;
}
int main(int argc, char *argv[])
{
double arg,c1,c2,blad;
arg= 3.113408;
c1 = cosh(arg);
c2 = mycosh(arg);
blad = abs(c2-c1)/c2;
printf( " z=%f; cosh(z)=%f; mycosh(z)=%f; blad=%.90f \n",arg,c1,c2,blad);
system("PAUSE");
return 0;
}
[/cpp]
