Métodos numéricos.

...

bufer1=ec[2][2];

bufer2=(-1)*ec[2][2];

for(i=1;i

for(j=1;j

ec1[i][j]=ec[i][j];

}

}

for(i=1;i

ec1[1][i]=ec[1][i]*bufer1;

ec1[2][i]=ec[2][i]*bufer2;

}

suma_c1=ec1[1][1]+ec1[2][1];

suma_c2=ec1[1][2]+ec1[2][2];

suma_c3=ec1[1][3]+ec1[2][3];

A=suma_c3/suma_c1;

/* Haciendo lo mismo pero ahora eliminando a A */

bufer1=ec[2][1];

bufer2=(-1)*ec[1][1];

for(i=1;i

for(j=1;j

ec1[i][j]=ec[i][j];

}

}

for(i=1;i

ec1[1][i]=ec[1][i]*bufer1;

ec1[2][i]=ec[2][i]*bufer2;

}

suma_c1=ec1[1][1]+ec1[2][1];

suma_c2=ec1[1][2]+ec1[2][2];

suma_c3=ec1[1][3]+ec1[2][3];

B=suma_c3/suma_c2;

printf("\n\n%0.2fA+ %dB=%0.2f",suma_x,n,suma_y);

printf("\n%0.2fA+%0.2fB=%0.2f",suma_x2,suma_x,suma_xy);

/* Sustituyendo los valores de A y B en la forma de la ecuacion

la linea recta tendremos */

printf("\n\n");

printf("La ecuacion de la recta que mas se ajusta a los puntos es:");

printf("\n\n");

printf("_______________________________________________");

printf("\n");

printf("y=%fx+%f",A,B);

printf("\n");

printf("_______________________________________________");

printf("\n");

getch();

}

[pic 4]

---------------------------------------------------------------

Runge Kutta

#include

#include

float k1(), k2();

float L,R,V,Ti,H;

float Tf,K1,K2;

float i,j,Iant,Inva;

float k1(float Iant)

{

float x;

x= H * (-(R/L)*(Iant) + (V/L));

return(x);

}

float k2(float K1)

{

float x;

x= H * (-(R/L)*(Iant+K1) + (V/L));

return(x);

}

int main()

{

int cont, cont1, cont2;

cont=cont1=cont2=0;

i=j=Iant=Inva=0.0;

L=R=V=Ti=H=0.0;

clrscr();

gotoxy(35,1);printf("METODO DE RUNGE-KUTTA\n");

for(i=1; i

gotoxy(1,5);printf("Teclee el valor de la autoinductancia.............:\n");

gotoxy(52,5);scanf("%f",&L);

gotoxy(1,6);printf("Teclee el valor de la Resistencia.................:");

gotoxy(52,6);scanf("%f",&R);

gotoxy(1,7);printf("Teclee el valor del Voltaje.......................:");

gotoxy(52,7);scanf("%f",&V);

gotoxy(1,8);printf("Teclee el valor del tiempo inicial................:");

gotoxy(52,8);scanf("%f",&Ti);

gotoxy(1,9);printf("Teclee el valor del tiempo final..................:");

gotoxy(52,9);scanf("%f",&Tf);

gotoxy(1,10);printf("Teclee el valor del paso o incrementos del tiempo:");

gotoxy(52,10);scanf("%f",&H);

printf("Paso Tiempo T Corriente \n");

printf("....................................\n");

for(i=Ti; i

{

cont=cont+1;

cont1=cont1+1;

K1=k1(Iant);

K2=k2(K1);

Inva=Iant + (K1+K2)/2;

j=i+0.1;

printf("%d