//Programa RK4 para sistemas de ecuaciones#include #include #include using namespace std;// Definicion de parametros del mtod (ejemplo de 3 ecuaciones diferenciales ordinarias)const int N = 3; // Numero de ecuacionesconst double paso = 0.1; // tamaNo del pasoconst double t0 = 0.0; //tiempo inicialconst double tfinal = 10; //tiempo final// function declarationsdouble f(const double & x, const double y[], const int & i); void UnPasoDeRK4(const double & x, double y[], const double & h, const int & nNum);//definicione de funcionesdouble f(const double & x, const double y[], const int & i){ if (0 == i) return 10*(y[1]-y[0]); //dx/dt x=y[0] y=y[1] z=y[1] if (1 == i) return 28*y[0]-y[1]-y[0]*y[2]; //dy/dt if (2 == i) return -(8/3)*y[2]+y[0]*y[1];// dz/dt}void UnPasoDeRK4(const double & x, double y[], const double & h, const int & num){ double k1[num], k2[num], k3[num], k4[num]; // Pendientes double l1[num], l2[num], l3[num]; //truco int i; for (i=0; i