Metodo de newton raphson

package newton;

/**

*

* @author Alejandro

*/

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 {

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 for(j=i+1;j<=n;j++){

cte=(-a[j][i])/(a[i][i]);

for(k=i;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);

}

...