/* Tehtava 1 */

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

// Alkioiden lkm taulukossa (tulostusteknisista syista vain 20 alkiota)
#define MAX 20

int main(void) {

int a[MAX], i, summa=0, pos=0;

  /* alustetaan satunnaislukugeneraattori */
  srand(time(NULL));

  /* alustetaan taulukko satunnaisluvuilla väliltä -50 ... 49 */
  for(i=0; i<MAX; i++)
    a[i] = (rand() % 100) - 50;

  /* Tulostetaan taulukko */
  printf("Alkuperainen taulukko:\n");
  for(i=0; i<MAX; i++)
    printf("%d ", a[i]);
  printf("\n");

  /* kohta c: lasketaan kuinka moni taulukon alkioista on positiivinen */
  for(i=0; i<MAX; i++) 
    if ( 0<a[i] )
      pos++ ;

  printf("Positiivisia alkiota %d kpl.\n", pos);

  /* kohta b: lasketaan yhteen taulukon alkiot a[i], missä i=0, 5, 10, ... */
  for(i=0; i<MAX; i+=5)
    summa += a[i] ;

  printf("Alkioiden summa on %d.\n", summa);

  /* kohta a: korvataan kukin alkio neliollaan */
  for(i=0; i<MAX; i++)
    a[i] *= a[i];

  /* Tulostetaan taulukko */
  printf("Nelioity taulukko:\n");
  for(i=0; i<MAX; i++)
    printf("%d ", a[i]);
  printf("\n");

  return(0);

}





/* lasketaan eksponenttifunktio Taylorin sarjasta */

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

int main(void) {

int i, n ;

/*
  Mita kay laskun tarkkuudelle jos nama ovat float tyyppisia muuttujia?
*/
double x, summa, termi, raja ;

  printf("Anna x, termien maara ja pienin laskettava termi:\n") ;
  if ( scanf("%lf %d %lf", &x, &n, &raja) != 3 ) {
    printf("Senkin tomppeli.\n") ;
    printf("Anna luvut muodossa reaaliluku, kokonaisluku ja reaaliluku!\n") ;
    exit(1);
  }

  /*jos virheellinen termien lkm, lopetetaan */
  if (n <= 0) {
    printf("Termien lukumaaran oltava > 0\n") ;
    exit(2);
  }

  /* Alustetaan ensimmainen termi ja lasketaan Taylorin sarjan summa */
  summa = termi = 1.0;

  i=1;
  while ( (i<n) && (raja<termi) ) {
    summa += termi *= x/i++;
    /*
      yllaolevassa lasketaan seuraavasti:
      a) lasketaan apu = x/i
      b) kasvatetaan i:ta yhdella
      c) lasketaan termi = termi * apu
      d) lasketaan summa = summa + termi
      sijoittaminen suoritetaan siis oikealta vasemmalle
    */
  }

  printf("exp(%lf)=%.20lf kun summassa on %d termia.\n", x, summa, i);

  return(0);

}

/*

Esimerkkiajo:

tke@carmine:~/demo4$ ./4_2 
Anna x, termien maara ja pienin laskettava termi:
1 100 1e-50
exp(1.000000)=2.71828182845904553488 kun summassa on 43 termia.

*/


/*
  Tehtava 3.
  ratkaisulla x = (1-x)^(1/3) ohjelma toimii kun alkuarvo on valilta ]0,1[
  ratkaisulla x = 1 - x^3 iterointi ei vaikuta onnistuvan millaan alkuarvoilla
  ylla siis ^ on potenssiin korotus joten ^(1/3) on kuutiojuuri
*/

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

#define MAX_LOOP 10000

int main(void) {

int n=0, c;  double x0, x1;

  printf("\n");
  printf("Funktiot ovat:\n");
  printf(" 1) x = (1-x)^(1/3)\n");
  printf(" 2) x = 1-x^3\n");

  printf("Kumpaa funktiota kokeillaan, 1 vai 2 ?\n") ;
  if ( (scanf("%d", &c)!=1) || ( c!=1 && c!=2 ) ) {
    printf("Valitse 1 tai 2.\n");
    exit(1);
  }

  printf("Anna alkuarvo:\n");
  scanf("%lf", &x0);

  /*
    Seuraavassa kaytetaan XXX ? YYY : ZZZ -notaatiota.
    Siis jos XXX on tosi niin YYY muuten ZZZ
    Eli jos c==1 niin x1 = cbrt(1 - x0) muuten x1 = 1 - x0*x0*x0;
  */
  x1 = c==1 ? cbrt(1 - x0) : 1 - x0*x0*x0;

  while ( (fabs(x1-x0)>1.0e-5) && (n<MAX_LOOP) ) {
    x0 = x1 ;
    x1 = c==1 ? cbrt(1 - x0) : 1 - x0*x0*x0;
    n++;
  }

  if ( c==1 )
    printf("x=%lf, (1-x)^(1/3) = %lf, n=%d\n", x1, cbrt(1 - x1), n);
  else
    printf("x=%lf, 1-x^3 = %lf, n=%d\n", x1, 1 - x1*x1*x1, n);

  if ( n == MAX_LOOP )
    printf("Ratkaisua ei ilmeisestikaan loytynyt.\n");

  return(0);

}