Código Fuente en Java para la resolución de Métodos Numéricos

por

Hace mucho en una clase de programación de JAVA me pidieron realizar un programa que resolviera incógnitas por diverso métodos numéricos entre los cuales se encuentran:

  • Método de Bisección
  • Método de NEWTON-RAPHSON
  • Método de GAUSS
  • Método de GAUSS-SEIDEL PARA 3 ECUACIONES
  • Método de GAUSS-JORDAN
  • Método de LU PARA 3 ECUACIONES
  • Método de Interpolación de NEWTON
  • Método de Interpolación de LAGRANGE

Y pues eh decidido compartir el código de este proyecto, la verdad el código no es del todo mío… Fue facilitado por un compañero, porque para ser honestos a mi no me acaba de agradar la programación orientada a objetos, prefiero la estructurada, pero en gustos se rompen generos…

NOTA:el codigo fuente lo puedes descargar a continuación:

>>> DESCARGAR Codigo Fuente <<<

El código fuente se muestra a continuación:

  1. import java.lang.Math;
  2. import java.io.*;
  3. //Clase ProyectoJava
  4. public class ProyectoJava{
  6. //Metodo Constructor
  7. public static void main(String args[]){
  8. ProyectoJava Proyecto = new ProyectoJava();
  9. Proyecto.menu(); //solo se manda a llamar a menu desde aqui, ya que menu llama a los demas metodos posteriromente
  10. }
  13. //(1)Metodo Biseccion
  14. public void MetodoBiseccion(){
  15. double a;
  16. double b;
  17. double tol;
  18. System.out.println("\t\t\t\"METODO DE BISECCION\"");
  19. System.out.println("Extremo Izquierdo: ");
  20. a=lee();
  21. System.out.println("Extremo Derecho: ");
  22. b=lee();
  23. System.out.println("Tolerancia: ");
  24. tol=lee();
  25. double c;
  26. do{
  27. c=(a+b)/2.0;
  28. if(((c*c-5)*(a*a-5))<0){
  29. b=c;
  30. }
  31. else{
  32. a=c;
  33. }
  34. }while(Math.abs(a-b)>tol);
  35. System.out.println("La raiz es: "+c);
  36. }
  39. //(2)Metodo NewRaphson
  40. public void MetodoNewtonRaphson(){
  41. double a;
  42. double tol;
  43. double b;
  44. double c;
  45. System.out.println("\t\t\t\"METODO DE NEWTON-RAPHSON\"");
  46. System.out.println("Primera Aproximacion: ");
  47. a=lee();
  48. System.out.println("Tolerancia: ");
  49. tol=lee();
  50. do{
  51. b=a-(a*a-a-2)/(2*a-1);
  52. c=Math.abs(a-b);
  53. a=b;
  54. }while(c>tol);
  55. System.out.println("La raiz es: "+b);
  56. }
  59. //(3) metodo gauss
  60. public void MetodoGauss(){
  61. int i,j,k,n;
  62. double a[][] = new double [5][6];
  63. double cte,x1,x2,x3;
  64. double x[] = new double [5];
  65. System.out.println("\t\t\t\"METODO DE GAUSS\"");
  66. System.out.println("Dame el numero de incognitas (de hasta 5) ");
  67. n=leeint();
  68. System.out.println("Ingrese coeficientes");
  69. for(i=0;i<n;i++)
  70. {
  71. System.out.println("Renglon "+(i+1));
  72. for(j=0;j<=n;j++)
  73. {
  74. System.out.println(" Ingrese a "+(i+1)+" "+(j+1));
  75. a[i][j]=lee();
  76. }
  77. }
  78. for(i=0;i<n-1;i++){
  79. for(j=i+1;j<=n;j++){
  80. cte=(-a[j][i])/(a[i][i]);
  81. for(k=i;k<n+1;k++){
  82. a[j][k]=((a[i][k])*cte)+a[j][k];
  83. }
  84. }
  85. }
  86. x3=a[n-1][n]/a[n-1][n-1];
  87. x2=(a[n-2][n]-x3*a[n-2][n-1])/a[n-2][n-2];
  88. x1=(a[n-3][n]-x2*a[n-3][n-2]-x3*a[n-3][n-1])/a[n-3][n-3];
  89. System.out.println("x0= "+x1+" \nx1= " +x2+" \nx2= " +x3);
  90. }
  94. //(4)gauss seidel
  95. public void MetodoGaussSeidel(){
  96. double x0,x1,x2,tol,e;
  97. int i,j;
  98. double a[][]=new double [3][4];
  99. System.out.println("\t\t\t\"METODO DE GAUSS-SEIDEL 3 ECUACIONES\"");
  100. System.out.println("Ingrese tolerancia");
  101. tol=lee();
  102. System.out.println("Ingrese coeficientes");
  103. for(i=0;i<3;i++)
  104. {
  105. System.out.println("Renglon "+(i+1));
  106. for(j=0;j<=3;j++)
  107. {
  108. System.out.println(" Ingrese a "+(i+1)+" "+(j+1));
  109. a[i][j]=lee();
  110. }
  111. }
  112. x1=0.0;
  113. x2=0.0;
  114. do{
  115. e=x1;
  116. x0=(a[0][3]-x1*a[0][1]-x2*a[0][2])/a[0][0];
  117. x1=(a[1][3]-x0*a[1][0]-x2*a[1][2])/a[1][1];
  118. x2=(a[2][3]-x0*a[2][0]-x1*a[2][1])/a[2][2];
  119. }while(Math.abs(e-x1)>tol);
  120. System.out.println("x0= "+x0+" \nx1= " +x1+" \nx2= " +x2);
  121. }
  124. //(5) gauss jordan
  125. public void MetodoGaussJordan(){
  126. int i,j,k,n;
  127. double a[][] = new double [5][6];
  128. double cte;
  129. double x[] = new double [5];
  130. System.out.println("\t\t\t\"METODO DE GAUSS-JORDAN\"");
  131. System.out.println("No. de incognitas (maximo 5): ");
  132. n=leeint();
  133. System.out.println("Dame los coeficientes: ");
  134. for(i=0;i<n;i++)
  135. {
  136. System.out.println("Renglon "+(i+1));
  137. for(j=0;j<=n;j++)
  138. {
  139. System.out.println(" Ingrese a "+(i+1)+" "+(j+1));
  140. a[i][j]=lee();
  141. }
  142. }
  143. for(i=0;i<n-1;i++){
  144. for(j=i+1;j<=n;j++){
  145. cte=(-a[j][i])/(a[i][i]);
  146. for(k=i;k<n+1;k++){
  147. a[j][k]=((a[i][k])*cte)+a[j][k];
  148. }
  149. }
  150. }
  151. for(i=0;i<n;i++)
  152. {
  153. System.out.println("Renglon "+(i+1));
  154. for(j=0;j<=n;j++)
  155. {
  156. System.out.println(" a "+(i+1)+" "+(j+1)+" =" +a[i][j]);
  157. }
  158. }
  159. }
  162. //(6) LU
  163. public void MetodoLu(){
  164. int i,j;
  165. double Y1,Y2,Y3,X1,X2,X3;
  166. double A[][] = new double [3][3];
  167. double L[][] = new double [3][3];
  168. double U[][] = new double [3][3];
  169. double B[] = new double [3];
  170. System.out.println("\t\t\t\"METODO DE LU 3 ECUACIONES\"");
  171. System.out.println("Ingrese coeficientes d");
  172. for(i=0;i<3;i++)
  173. {
  174. System.out.println("Renglon "+(i+1));
  175. for(j=0;j<3;j++)
  176. {
  177. System.out.println(" Ingrese a "+(i+1)+" "+(j+1));
  178. A[i][j]=lee();
  179. }
  180. }
  181. System.out.println("Ingrese los terminos independientes");
  182. for(i=0;i<3;i++){
  183. System.out.println("Termino "+(i+1));
  184. B[i]=lee();
  185. }
  186. for(i=0;i<3;i++){
  187. L[i][0]=A[i][0];
  188. U[i][i]=1;
  189. }
  190. U[0][1]=A[0][1]/L[0][0];
  191. L[1][1]=A[1][1]-L[1][0]*U[0][1];
  192. L[2][1]=A[2][1]-L[2][0]*U[0][1];
  193. U[0][2]=A[0][2]/L[0][0];
  194. U[1][2]=(A[1][2]-L[1][0]*U[0][2])/L[1][1];
  195. L[2][2]=A[2][2]-L[2][0]*U[0][2]-L[2][1]*U[1][2];
  196. Y1=B[0]/L[0][0];
  197. Y2=(B[1]-L[1][0]*Y1)/L[1][1];
  198. Y3=(B[2]-L[2][0]*Y1-L[2][1]*Y2)/L[2][2];
  199. X3=Y3;
  200. X2=Y2-U[1][2]*X3;
  201. X1=Y1-U[0][1]*X2-U[0][2]*X3;
  202. System.out.println("El resultado es:\nx0= "+X1+" \nx1= " +X2+" \nx2= " +X3);
  203. }
  206. //(7) Interpolacion Newton
  207. public void MetodoInterNewton(){
  208. double a[][] = new double [5][2];
  209. double x,y,fx1x0,fx2x1,fx3x2,fx4x3,fx2x1x0,fx3x2x1,fx4x3x2,fx3x2x1x0,fx4x3x2x1,fx4x3x2x1x0;
  210. int i;
  211. System.out.println("\t\t\t\"INTERPOLACION DE DIFERENCIAS DE NEWTON P/ 5 PTOS.\"");
  212. System.out.println("Valor a interpolar: ");
  213. x=lee();
  214. System.out.println("Dame los 5 pares de puntos");
  215. for(i=0;i<5;i++){
  216. System.out.println("Dame x "+i);
  217. a[i][0]=lee();
  218. System.out.println("Dame f(x) "+i);
  219. a[i][1]=lee();
  220. }
  221. fx1x0=(a[1][1]-a[0][1])/(a[1][0]-a[0][0]);
  222. fx2x1=(a[2][1]-a[1][1])/(a[2][0]-a[1][0]);
  223. fx3x2=(a[3][1]-a[2][1])/(a[3][0]-a[2][0]);
  224. fx4x3=(a[4][1]-a[3][1])/(a[4][0]-a[3][0]);
  225. fx2x1x0=(fx2x1-fx1x0)/(a[2][0]-a[0][0]);
  226. fx3x2x1=(fx3x2-fx2x1)/(a[3][0]-a[1][0]);
  227. fx4x3x2=(fx4x3-fx3x2)/(a[4][0]-a[2][0]);
  228. fx3x2x1x0=(fx3x2x1-fx2x1x0)/(a[3][0]-a[0][0]);
  229. fx4x3x2x1=(fx4x3x2-fx3x2x1)/(a[3][0]-a[0][0]);
  230. fx4x3x2x1x0=(fx4x3x2x1-fx3x2x1x0)/(a[4][0]-a[0][0]);
  231. y=a[0][1]+fx1x0*(x-a[0][0])+fx2x1x0*(x-a[0][0])*(x-a[1][0])+fx3x2x1x0*(x-a[0][0])*(x-a[1][0])*(x-a[2][0])+fx4x3x2x1x0*(x-a[0][0])*(x-a[1][0])*(x-a[2][0])*(x-a[3][0]);
  232. System.out.println("f(x) en ese punto es: "+y);
  233. }
  236. //(8) Interpolacion Lagrange
  237. public void MetodoInterLagrange(){
  238. int n,i;
  239. double x,y;
  240. double a[][] = new double [4][2];
  241. System.out.println("\t\t\t\"METODO DE INTERPOLACION LAGRANGE 2-4 PTOS.\"");
  242. do{
  243. System.out.println("Dame el numero de puntos");
  244. n=leeint();
  245. }while(n<2 || n>4);
  246. System.out.println("Dame los pares de puntos");
  247. for(i=0;i<n;i++){
  248. System.out.println("Dame x "+i);
  249. a[i][0]=lee();
  250. System.out.println("Dame f(x) "+i);
  251. a[i][1]=lee();
  252. }
  253. System.out.println("Dame el valor a interpolar ");
  254. x=lee();
  255. switch(n){
  256. case 2:
  257. y=(((((x-a[1][0])*a[0][1])/(a[0][0]-a[1][0]))+(((x-a[0][0])*a[1][1])/(a[1][0]-a[0][0]))));
  258. System.out.println("f(x) en ese punto es: "+y);
  259. break;
  260. case 3:
  261. y=((x-a[1][0])*(x-a[2][0])*a[0][1])/((a[0][0]-a[1][0])*(a[0][0]-a[2][0]))+((x-a[0][0])*(x-a[2][0])*a[1][1])/((a[1][0]-a[0][0])*(a[1][0]-a[2][0]))+((x-a[1][0])*((x-a[0][0])*a[2][1])/((a[2][0]-a[0][0])*(a[2][0]-a[1][0])));
  262. System.out.println("f(x) en ese punto es: "+y);
  263. break;
  264. case 4:
  265. y=((x-a[1][0])*(x-a[3][0])*(x-a[2][0])*a[0][1])/((a[0][0]-a[1][0])*(a[0][0]-a[2][0])*(a[0][0]-a[3][0]))+((x-a[0][0])*(x-a[2][0])*(x-a[3][0])*a[1][1])/((a[1][0]-a[0][0])*(a[1][0]-a[2][0])*(a[1][0]-a[3][0]))+((x-a[0][0])*(x-a[1][0])*((x-a[3][0])*a[2][1])/((a[2][0]-a[0][0])*(a[2][0]-a[1][0])*(a[2][0]-a[3][0]))+((x-a[0][0])*(x-a[1][0])*((x-a[2][0])*a[
  266. 3][1])/((a[3][0]-a[0][0])*(a[3][0]-a[1][0])*(a[3][0]-a[2][0]))));
  267. System.out.println("f(x) en ese punto es: "+y);
  268. break;
  269. default:
  270. System.out.println("INVALIDO");
  271. break;
  272. }
  273. }
  276. //para leer desde teclado
  277. public double lee(){
  278. double num;
  279. try{
  280. InputStreamReader isr = new InputStreamReader (System.in);
  281. BufferedReader br = new BufferedReader(isr);
  282. String sdato;
  283. sdato = br.readLine();
  284. num = Double.parseDouble(sdato);
  285. }
  286. catch(IOException ioe){
  287. num=0.0;
  288. }
  289. return num;
  290. }
  292. //para leer un entero
  293. public int leeint(){
  294. int num;
  295. try{
  296. InputStreamReader isr = new InputStreamReader (System.in);
  297. BufferedReader br = new BufferedReader(isr);
  298. String sdato;
  299. sdato = br.readLine();
  300. num = Integer.parseInt(sdato);
  301. }
  302. catch(IOException ioe){
  303. num=0;
  304. }
  305. return num;
  306. }
  309. //para salir del programa
  310. public int Fuera(){
  311. int sal;
  312. System.out.println("\n\n\nSI DESEAS OTRO METODO PRESIONA [1]");
  313. sal=leeint();
  314. return sal;
  315. }
  320. //despliega menu
  321. public void menu(){
  322. int a;
  323. int p;
  324. do{
  325. do{
  326. System.out.println("\n\n\t\t\tMETODOS NUMERICOS\n\n");
  327. System.out.println("\t1.-Biseccion\n\t2.-Newton-Raphson\n\t3.-Gauss\n\t4.-Gauss-Seidel\n\t5.-Gauss-Jordan\n\t6.-LU\n\t7.-Interpolacion Newton\n\t8.-Interpolacion Lagrange");
  328. System.out.println("\n\nEscoja el numero del metodo que desea usar:");
  329. a=leeint();
  330. }while(a<1 || a>8);
  331. switch(a){
  332. case 1:
  333. MetodoBiseccion(); //manda a llamara a cada uno de los metodos
  334. p=Fuera();
  335. break;
  336. case 2:
  337. MetodoNewtonRaphson();
  338. p=Fuera();
  339. break;
  340. case 3:
  341. MetodoGauss();
  342. p=Fuera();
  343. break;
  344. case 4:
  345. MetodoGaussSeidel();
  346. p=Fuera();
  347. break;
  348. case 5:
  349. MetodoGaussJordan();
  350. p=Fuera();
  351. break;
  352. case 6:
  353. MetodoLu();
  354. p=Fuera();
  355. break;
  356. case 7:
  357. MetodoInterNewton();
  358. p=Fuera();
  359. break;
  360. case 8:
  361. MetodoInterLagrange();
  362. p=Fuera();
  363. break;
  364. default:
  365. System.out.println("Opcion incorrecta");
  366. p=1;
  367. break;
  368. }
  369. }while(p==1);
  370. }
  372. }
import java.lang.Math;
import java.io.*;
//Clase ProyectoJava
public class  ProyectoJava{
  
  //Metodo Constructor
  public static void main(String args[]){
  ProyectoJava Proyecto = new ProyectoJava();
  Proyecto.menu();	//solo se manda a llamar a menu desde aqui, ya que menu llama a los demas metodos posteriromente
  }
    
  
    //(1)Metodo Biseccion	
  public void MetodoBiseccion(){
    double a;
    double b;
    double tol;
    System.out.println("\t\t\t\"METODO DE BISECCION\"");
    System.out.println("Extremo Izquierdo: ");
    a=lee();
    System.out.println("Extremo Derecho: ");
    b=lee();
    System.out.println("Tolerancia: ");
    tol=lee();
    double c;
    do{
    c=(a+b)/2.0;
    if(((c*c-5)*(a*a-5))<0){
      b=c;
    }
    else{
      a=c;
    }
      }while(Math.abs(a-b)>tol);
      System.out.println("La raiz es: "+c);			
  }


  //(2)Metodo NewRaphson
  public void MetodoNewtonRaphson(){
    double a;
    double tol;
    double b;
    double c; 
    System.out.println("\t\t\t\"METODO DE NEWTON-RAPHSON\"");
    System.out.println("Primera Aproximacion: ");
    a=lee();
    System.out.println("Tolerancia: ");
    tol=lee();
    do{
      b=a-(a*a-a-2)/(2*a-1);
      c=Math.abs(a-b);
      a=b;	
      	}while(c>tol);
      	System.out.println("La raiz es: "+b);
    }


  //(3) metodo gauss
  public void MetodoGauss(){
    int i,j,k,n;
    double a[][] = new double [5][6];
    double cte,x1,x2,x3;
    double x[] = new double [5];
    System.out.println("\t\t\t\"METODO DE GAUSS\"");
    System.out.println("Dame el numero de incognitas (de hasta 5) ");
    n=leeint();
    System.out.println("Ingrese coeficientes");
    for(i=0;i<n;i++)
    {
      System.out.println("Renglon "+(i+1));
      for(j=0;j<=n;j++)
      {
        System.out.println(" Ingrese a "+(i+1)+" "+(j+1));
        a[i][j]=lee();
      }
    }
    for(i=0;i<n-1;i++){
      for(j=i+1;j<=n;j++){
       cte=(-a[j][i])/(a[i][i]);
          for(k=i;k<n+1;k++){
       a[j][k]=((a[i][k])*cte)+a[j][k];
       }
       }
    }
   	x3=a[n-1][n]/a[n-1][n-1];
   	x2=(a[n-2][n]-x3*a[n-2][n-1])/a[n-2][n-2];
   	x1=(a[n-3][n]-x2*a[n-3][n-2]-x3*a[n-3][n-1])/a[n-3][n-3];
    System.out.println("x0= "+x1+" \nx1= " +x2+" \nx2= " +x3);
    }

  
  
  //(4)gauss seidel
  public void MetodoGaussSeidel(){
    double x0,x1,x2,tol,e;
    int i,j;
    double a[][]=new double [3][4];
    System.out.println("\t\t\t\"METODO DE GAUSS-SEIDEL 3 ECUACIONES\"");
    System.out.println("Ingrese tolerancia");
    tol=lee();
    System.out.println("Ingrese coeficientes");
    for(i=0;i<3;i++)
    {
      System.out.println("Renglon "+(i+1));
      for(j=0;j<=3;j++)
      {
        System.out.println(" Ingrese a "+(i+1)+" "+(j+1));
        a[i][j]=lee();
      }
    }
    x1=0.0;
    x2=0.0;
    do{
      e=x1;
      x0=(a[0][3]-x1*a[0][1]-x2*a[0][2])/a[0][0];
      x1=(a[1][3]-x0*a[1][0]-x2*a[1][2])/a[1][1];
      x2=(a[2][3]-x0*a[2][0]-x1*a[2][1])/a[2][2];
    }while(Math.abs(e-x1)>tol);
    System.out.println("x0= "+x0+" \nx1= " +x1+" \nx2= " +x2);
    }
  
  
  //(5) gauss jordan
  public void MetodoGaussJordan(){
    int i,j,k,n;
    double a[][] = new double [5][6];
    double cte;
    double x[] = new double [5];
    System.out.println("\t\t\t\"METODO DE GAUSS-JORDAN\"");
    System.out.println("No. de incognitas (maximo 5):   ");
    n=leeint();
    System.out.println("Dame los coeficientes: ");
    for(i=0;i<n;i++)
    {
      System.out.println("Renglon "+(i+1));
      for(j=0;j<=n;j++)
      {
        System.out.println(" Ingrese a "+(i+1)+" "+(j+1));
        a[i][j]=lee();
      }
    }
    for(i=0;i<n-1;i++){
      for(j=i+1;j<=n;j++){
       cte=(-a[j][i])/(a[i][i]);
          for(k=i;k<n+1;k++){
       a[j][k]=((a[i][k])*cte)+a[j][k];
       }
       }
    }
    for(i=0;i<n;i++)
    {
      System.out.println("Renglon "+(i+1));
      for(j=0;j<=n;j++)
      {
        System.out.println(" a "+(i+1)+" "+(j+1)+" =" +a[i][j]);
      }
    }
    }

  
  //(6) LU
  public void MetodoLu(){
    int i,j;
    double Y1,Y2,Y3,X1,X2,X3;
    double A[][] = new double [3][3];
    double L[][] = new double [3][3];
    double U[][] = new double [3][3];
    double B[] = new double [3];
    System.out.println("\t\t\t\"METODO DE LU    3 ECUACIONES\"");
    System.out.println("Ingrese coeficientes d");
    for(i=0;i<3;i++)
    {
      System.out.println("Renglon "+(i+1));
      for(j=0;j<3;j++)
      {
        System.out.println(" Ingrese a "+(i+1)+" "+(j+1));
        A[i][j]=lee();
      }
    }
    System.out.println("Ingrese los terminos independientes");
    for(i=0;i<3;i++){
      System.out.println("Termino  "+(i+1));
      B[i]=lee();
      }
    for(i=0;i<3;i++){
   			L[i][0]=A[i][0];
   			U[i][i]=1;
      }
    U[0][1]=A[0][1]/L[0][0];
    L[1][1]=A[1][1]-L[1][0]*U[0][1];
    L[2][1]=A[2][1]-L[2][0]*U[0][1];
    U[0][2]=A[0][2]/L[0][0];
    U[1][2]=(A[1][2]-L[1][0]*U[0][2])/L[1][1];
    L[2][2]=A[2][2]-L[2][0]*U[0][2]-L[2][1]*U[1][2];
    Y1=B[0]/L[0][0];
    Y2=(B[1]-L[1][0]*Y1)/L[1][1];
    Y3=(B[2]-L[2][0]*Y1-L[2][1]*Y2)/L[2][2];
    X3=Y3;
    X2=Y2-U[1][2]*X3;
    X1=Y1-U[0][1]*X2-U[0][2]*X3;
    System.out.println("El resultado es:\nx0= "+X1+" \nx1= " +X2+" \nx2= " +X3);
    }

  
  //(7) Interpolacion Newton
  public void MetodoInterNewton(){
    double a[][] = new double [5][2];
    double x,y,fx1x0,fx2x1,fx3x2,fx4x3,fx2x1x0,fx3x2x1,fx4x3x2,fx3x2x1x0,fx4x3x2x1,fx4x3x2x1x0;
    int i;
    System.out.println("\t\t\t\"INTERPOLACION DE DIFERENCIAS DE NEWTON P/ 5 PTOS.\"");
    System.out.println("Valor a interpolar: ");
    x=lee();
    System.out.println("Dame los 5 pares de puntos");
    for(i=0;i<5;i++){
        System.out.println("Dame x "+i);
        a[i][0]=lee();
        System.out.println("Dame f(x) "+i);
        a[i][1]=lee();
        }		
    fx1x0=(a[1][1]-a[0][1])/(a[1][0]-a[0][0]);
    fx2x1=(a[2][1]-a[1][1])/(a[2][0]-a[1][0]);
    fx3x2=(a[3][1]-a[2][1])/(a[3][0]-a[2][0]);
    fx4x3=(a[4][1]-a[3][1])/(a[4][0]-a[3][0]);
    fx2x1x0=(fx2x1-fx1x0)/(a[2][0]-a[0][0]);
    fx3x2x1=(fx3x2-fx2x1)/(a[3][0]-a[1][0]);
    fx4x3x2=(fx4x3-fx3x2)/(a[4][0]-a[2][0]);
    fx3x2x1x0=(fx3x2x1-fx2x1x0)/(a[3][0]-a[0][0]);
    fx4x3x2x1=(fx4x3x2-fx3x2x1)/(a[3][0]-a[0][0]);
    fx4x3x2x1x0=(fx4x3x2x1-fx3x2x1x0)/(a[4][0]-a[0][0]);
    y=a[0][1]+fx1x0*(x-a[0][0])+fx2x1x0*(x-a[0][0])*(x-a[1][0])+fx3x2x1x0*(x-a[0][0])*(x-a[1][0])*(x-a[2][0])+fx4x3x2x1x0*(x-a[0][0])*(x-a[1][0])*(x-a[2][0])*(x-a[3][0]);
          System.out.println("f(x) en ese punto es: "+y);
    }


  //(8) Interpolacion Lagrange
  public void MetodoInterLagrange(){
    int n,i;
    double x,y;
    double a[][] = new double [4][2];
    System.out.println("\t\t\t\"METODO DE INTERPOLACION LAGRANGE  2-4 PTOS.\"");
    do{
    System.out.println("Dame el numero de puntos");
    n=leeint();	
    }while(n<2 || n>4);
    System.out.println("Dame los  pares de puntos");
    for(i=0;i<n;i++){
        System.out.println("Dame x "+i);
        a[i][0]=lee();
        System.out.println("Dame f(x) "+i);
        a[i][1]=lee();
        }
    System.out.println("Dame el valor a interpolar ");
    x=lee();
    switch(n){
      case 2:
      y=(((((x-a[1][0])*a[0][1])/(a[0][0]-a[1][0]))+(((x-a[0][0])*a[1][1])/(a[1][0]-a[0][0]))));
      System.out.println("f(x) en ese punto es: "+y);
      break;
      case 3:
      y=((x-a[1][0])*(x-a[2][0])*a[0][1])/((a[0][0]-a[1][0])*(a[0][0]-a[2][0]))+((x-a[0][0])*(x-a[2][0])*a[1][1])/((a[1][0]-a[0][0])*(a[1][0]-a[2][0]))+((x-a[1][0])*((x-a[0][0])*a[2][1])/((a[2][0]-a[0][0])*(a[2][0]-a[1][0])));
      System.out.println("f(x) en ese punto es: "+y);
      break;
      case 4:
      y=((x-a[1][0])*(x-a[3][0])*(x-a[2][0])*a[0][1])/((a[0][0]-a[1][0])*(a[0][0]-a[2][0])*(a[0][0]-a[3][0]))+((x-a[0][0])*(x-a[2][0])*(x-a[3][0])*a[1][1])/((a[1][0]-a[0][0])*(a[1][0]-a[2][0])*(a[1][0]-a[3][0]))+((x-a[0][0])*(x-a[1][0])*((x-a[3][0])*a[2][1])/((a[2][0]-a[0][0])*(a[2][0]-a[1][0])*(a[2][0]-a[3][0]))+((x-a[0][0])*(x-a[1][0])*((x-a[2][0])*a[
3][1])/((a[3][0]-a[0][0])*(a[3][0]-a[1][0])*(a[3][0]-a[2][0]))));
      System.out.println("f(x) en ese punto es: "+y);
      break;
      default:
      System.out.println("INVALIDO");
      break;	
    }
  }

  
  //para leer desde teclado
  public double lee(){
    double num;
    try{
      InputStreamReader isr = new InputStreamReader (System.in);
      BufferedReader br = new BufferedReader(isr);
      String sdato;
      sdato = br.readLine();
      num = Double.parseDouble(sdato);
    }
    catch(IOException ioe){
      num=0.0;
    }
    return num;
    }

  //para  leer un entero
  public int leeint(){
    int num;
    try{
      InputStreamReader isr = new InputStreamReader (System.in);
      BufferedReader br = new BufferedReader(isr);
      String sdato;
      sdato = br.readLine();
      num = Integer.parseInt(sdato);
    }
    catch(IOException ioe){
      num=0;
    }
    return num;
    }
  
  
  //para salir del programa
  public int Fuera(){
    int sal; 
    System.out.println("\n\n\nSI DESEAS OTRO METODO PRESIONA [1]");
    sal=leeint();
    return sal;
    }




  //despliega menu
  public void menu(){
    int a;
    int p;
    do{
    do{
    System.out.println("\n\n\t\t\tMETODOS NUMERICOS\n\n");
    System.out.println("\t1.-Biseccion\n\t2.-Newton-Raphson\n\t3.-Gauss\n\t4.-Gauss-Seidel\n\t5.-Gauss-Jordan\n\t6.-LU\n\t7.-Interpolacion Newton\n\t8.-Interpolacion Lagrange");
    System.out.println("\n\nEscoja el numero del metodo que desea usar:");
    a=leeint();
    }while(a<1 || a>8);
    switch(a){
      case 1:
      MetodoBiseccion();  //manda a llamara a cada uno de los metodos
      p=Fuera();
      break;
      case 2:
      MetodoNewtonRaphson();
      p=Fuera();
      break;
      case 3:
      MetodoGauss();
      p=Fuera();
      break;
      case 4:
      MetodoGaussSeidel();
      p=Fuera();
      break;
      case 5:
      MetodoGaussJordan();
      p=Fuera();
      break;
      case 6:
      MetodoLu();
      p=Fuera();
      break;
      case 7:
      MetodoInterNewton();
      p=Fuera();
      break;
      case 8:
      MetodoInterLagrange();
      p=Fuera();
      break;
      default:
      System.out.println("Opcion incorrecta");
      p=1;
      break;
      }
    }while(p==1);
    }

}

 

 

42 comentarios en “Código Fuente en Java para la resolución de Métodos Numéricos”

    • Miguel:

      Lo que realmente deberia darte ladilla es Codear TU MISMO el source, en definitva eso si te daria ladilla. Ahora me llamas Pajuo? Digamos que si hubieses puesto atencion veras que en la parte superior del Source hay una opccion que dice PLAIN TEXT donde si tu le das clic veras que te pasa todo a texto plano listo para copiar en tu editor favorio. Ahora si eso no te basta podrias descargar el source directamente en RAR (indicado con letras rojas).

      Saludos!

      Responder

    • @marco:

      Si mi memoria no me falla todos los lenguajes de programación incorporan las funciones básicas tales, como sin(x), cos(x), etc. Seria cuestión de que le dieras una leída a la documentación.

      Saludos.

      Responder

  16. Pues, solo observo que mencionas que no te agrada el POO (Abstracción, Herencia, Polimorfismo y Encapsulamiento) y que prefieres el paradigma estructurado. No entiendo a que te refieres en tu ejemplo, ya que Java como cualquier lenguaje con capacidad POO puede escribirse de las dos maneras. Tu aportación es totalmente ESTRUCTURADA, pero igual ayuda y sirve bastante, ya que lo que importa es el algoritmo de cada método, con eso es suficiente. GRACIAS x tu APORTE

    Responder

Responder a luis alberto Cancelar respuesta