METODOS NUMERICOS

phantom
04 de Abril del 2005
Hola, como estan todos, necesito un poco de ayuda, alguien me podria pasar los metodos de punto fijo, newton rhapson y secante. por favor.

Los necesito. Muchas Gracias.

Que esten bien

angelucha84
04 de Abril del 2005
#include<stdio.h>
#include<conio.h>
#include<math.h>

float funcion(float *p);
float derivada(float *p);
void main()
{
float po,p,TOL,fp,fpo,dpo,a,b,fa,fb,ban=0;
int No,i;//ban=0;
clrscr();
gotoxy(8,3);printf("Programa para encontrar una solucion de f(x)=0 dada la funcion:");
gotoxy(30,5);printf("x^3+4x^2-10 = 0");
gotoxy(15,8);printf("Dame aproximacion inicial: ");scanf("%f",&po);
gotoxy(15,9);printf("Dame el margen de error: ");scanf("%f",&TOL);
gotoxy(15,10);printf("Dame el numero maximo de iteraciones: ");scanf("%d",&No);
i=1;
while(i<=No)
{
fpo=funcion(&po);
dpo=derivada(&po);

p=po-fpo/dpo;
if(fabs(p-po)<TOL)
{
ban=1;break;
}
else
i++;
po=p;
}
if(ban==1)
{
gotoxy(25,13);printf("x= %f para hacer f(x)=0",p);
gotoxy(20,15);printf("Procedimiento completado satisfactoriamente");
}
else
{
gotoxy(20,15);printf("El metodo fracaso despues de %d iteraciones",i-1);
gotoxy(20,16);printf("Procedimiento completado sin exito");
}
getch();
}

float funcion(float *po)
{ //Mi funcion es: xï3+4x'2-10 =0
float fpo;
fpo=pow(*po,3)+4*pow(*po,2)-10;
return fpo;
}

float derivada(float *po)
{ //Mi funcion es: xï3-2xï2-5=0
float dpo;
dpo=3*(pow(*po,2))+(8*(*po));
return dpo;
}

angelucha84
04 de Abril del 2005
#include<stdio.h>
#include<conio.h>
#include<math.h>

float funcion(float *po);
void main()
{
float TOL,po,p,TI,EXP,fpo,ap;
int No,i,ban=0;
clrscr();
gotoxy(30,1);printf("METODO DE PUNTO FIJO");
gotoxy(8,3);printf("Programa para encontrar una solucion de f(x)=0 dada la funcion:");
gotoxy(30,5);printf("1/2(10-x^3)^1/2 = 0");
gotoxy(15,8);printf("Dame aproximacion inicial: ");scanf("%f",&po);
gotoxy(15,9);printf("Dame el margen de error: ");scanf("%f",&TOL);
gotoxy(15,10);printf("Dame el numero maximo de iteraciones: ");scanf("%d",&No);
for(i=1;i<=No;i++)
{
fpo=funcion(&po);
p=fpo;
if(fabs(p-po)<TOL)
{
ban=1;break;
}
else
po=p;
}
if(ban==1)
{
gotoxy(25,13);printf("x= %f para que f(x)=0",p);
gotoxy(20,15);printf("Procedimiento completado satisfactoriamente");
}
else
{
gotoxy(20,15);printf("Procedimiento completado insatisfactoriamente");
}
getch();
}

float funcion(float *po)
{ //Mi funcion es: 1/2(10-x^3)^1/2=0
int cant,j,num,i;
float fpo;
fpo=0.5*pow(10-pow(*po,3),0.5);
return fpo;
}

angelucha84
04 de Abril del 2005
#include<stdio.h>
#include<conio.h>
#include<math.h>

float funcion(float *po);
void main()
{
float p0,p1,TOL,q0,q1,p,fp;
int No,i,ban=0;
clrscr();
printf("Dame aproximacion inicial po: ");scanf("%f",&p0);
printf("Dame aproximacion inicial pi: ");scanf("%f",&p1);
printf("Dame el margen de error: ");scanf("%f",&TOL);
printf("Dame el numero maximo de iteraciones: ");scanf("%d",&No);
i=2;
q0=funcion(&p0);
q1=funcion(&p1);

while(i<=No)
{
p=(p1-(q1*(p1-p0)))/(q1-q0);
if(fabs(p-p1)<TOL)
{
ban=1;break;
}
i++;
p0=p1;q0=q1;p1=p;
fp=funcion(&p);
q1=fp;
}
if(ban==1)
{printf("np: %fn",p);printf("Procedimiento completado satisfactoriamente");}
else
printf("nEl metodo fracaso despues de %d iteracionesnProcedimiento completado insatisfactoriamente",i-2);
getch();
}

float funcion(float *p)
{ //Mi funcion es: 3x+2=0
int cant,j,num,i;
float fp;
fp=(3*(*p))+2;
return fp;
}

angelucha84
04 de Abril del 2005
Ojalá k te sirvan

Unas por otras


Angélica [email protected]