Метод Рунге-Кутта 4-го порядка

#include <stdio.h>
#include <iostream>
#include <math.h>
using namespace std;
 
double f(double x, double y)
{
      return(2*y/x+x*x*pow(2.7,x));               
}
 
int main(void){
   int n,i;
   double x,y,h,end,k1,k2,k3,k4;
   
   x=-1;
   
   y=-2.33
   
   end=2
   
   h=0.3;   
 
   n=(end-x)/h;
   
   for(i=1;i<=n;i++){
      k1=h*f(x, y);
      k2=h*f(x+h/4, y+k1/4);
      k3=h*f(x+0.5*h, y+0.5*k2);
      k4=h*f(x+h, y+k1-2*k2+2*k3);
      x+=h;
      y+=(k1+4*k3+k4)/6;
      cout<<"x="<<x<<" "<<"y="<<y<<endl;
   }
 
}

 k2=h*f(x+h/2, y+(h/2)/k1);
      k3=h*f(x+0.5*h, y+0.5*k2*h);

#include <stdio.h>
#include <iostream>
#include <math.h>
using namespace std;
 
double f(double x, double y)
{
          return(2*y/x+x*x*pow(2.7,x));               
}
 
int main(void){
   int n,i;
   double x,y,h,end,k1,k2,k3,k4;
   
   x=1.1;
   
   y=0.341;
   
   end=2;
   
   h=0.3;   
 
   //n=(end-x)/h;
   
   for(i=1;i<=4;i++){
      k1=h*f(x, y);
      k2=h*f(x+h/4, y+k1/4);
      k3=h*f(x+0.5*h, y+0.5*k2);
      k4=h*f(x+h, y+k1-2*k2+2*k3);
      x+=h;
      y+=(k1+4*k3+k4)/6;
      cout<<"x="<<x<<" "<<"y="<<y<<endl;
      cout<<x*x*(pow(2.7,x)-2.7)<<endl;
   }
 system("PAUSE");
}