Modos normales de vibración de una barra elástica

Raíces múltiples

Dada la energía E, las abscisas de los puntos de retorno, son las raíces de la ecuación

cosh(qL)·cos (qL)=-1

public class Funcion extends Ecuacion{
    double f(double x){
      double y=cosh(x)*Math.cos(x)+1;
      return y;
    }
    double cosh(double x){
      return((Math.exp(x)+Math.exp(-x))/2);
    }
}

Calcula las raíces

public class Aplicacion {
           public static void main(String[] args) {
		         double[] raices=new Funcion().
hallarRaices(1.0, 50, 1.0);
                   System.out.print("Raíces de la ecuación");
                   for(int i=0; i<raices.length; i++){
                         System.out.print("  "+raices[i]);
                   }
                   System.out.println("");
}

Integral definida

Método de Simpson

La amplitud y(x) de los distintos puntos x de la barra en el modo normal de vibración n es:

$y_{n} (x) = A {(sinh (q_{n} x) - sin (q_{n} x)) - \frac{sinh (q_{n} L) + sin (q_{n} L)}{cosh (q_{n} L) + cos (q_{n} L)} (cosh (q_{n} x) - cos (q_{n} x))}$

El valor de la constante de proporcionalidad A es la escala vertical. Para que todos los modos de vibración estén dibujados a la misma escala, se calcula A de modo que

$\int_{0}^{L} y_{n}^{2} (x) d x = cte$

public class Funcion extends Simpson{
    double q;
	  double c;
	public Funcion(double q, double c){
		this.q=q;
		this.c=c;
	}
	public double f(double x){
        return (fOnda(x)*fOnda(x));
    }
	double fOnda(double x){
       double y=sinh(q*x)-Math.sin(q*x)-c*(cosh(q*x)-Math.cos(q*x));
       return y;
    }
}

public class Aplicacion {
  public static void main(String[] args) {
      double q=raices[iModo]; //modo de vibración
	  double c=(sinh(q)+Math.sin(q))/(cosh(q)+Math.cos(q));
	  double resultado=new Funcion(q, c).integral(0.0, Math.PI/2, 10);
        System.out.println("integral "+resultado);
  }
}