practica4_4.cpp

Topics Practices: Programming and Numerical Methods Practice 1: Introduction to C Practice 2

/*PRACTICA 4 EJERCICIO 4
El coeficiente de fricci髇 para el flujo de una suspensi髇 de part韈ulas 
fibrosas se puede relacionar con el n鷐ero de Reynolds mediante la siguiente 
ecuaci髇 emp韗ica:

1/sqrt(f) = ln(Re(sqrt(f))/k + (14 - 5.6/k)

En esta ecuaci髇, f es el coeficiente de fricci髇, Re es el n鷐ero de Reynolds y 
k es una constante determinada por la concentraci髇 de la suspensi髇. Para una 
suspensi髇 con 0.08 % de concentraci髇, se tiene que k = 0.28. 緾u醠 es el valor
de f si Re = 3750?  */


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

float f(float n)
{
      float k,Re;
      k=0.28;
      Re=3750;
      return (1/k)*log(Re*sqrt(n)) + 14 - 5.6/k - 1/sqrt(n);
}
main()
{
      float a,b,k;
      printf("Aberigue coeficiente de friccion para el flujo de una suspension de particulas\n");
      printf("ingrese el extremo izquierdo del intervalo:  ");
      scanf("%f", &a);
      printf("ingrese el extremo derecho del intervalo:  ");
      scanf("%f", &b);
      printf("f(%f) = %f\n", a, f(a));
      printf("f(%f) = %f\n", b, f(b));
      if(f(a)<0 && f(b)>0)
      {
                while(b-a>=1e-5)
                {
                      k=(a+b)/2;
                      if(f(k)>0)
                         b=k;
                      else
                         a=k;
                 }
      printf("la raiz es %f \n", k);
      }
      else 
      {
      printf("no cumple la hipotesis!!\n");
      }
printf("PRESIONE CUALQUIER TECLA PARA TERMINAR!!!!!!");      
getch();      
}