البته برنامش دارم ، اگه کسی بتونه توضیح بده ممنون میشم.
fourth order runge-kutta method for solving the following problem
//y'=f(x,y), y(x0)=y0
#include <iostrem.h>
#include <conio.h>
#include <math.h>
double f(double x,double y){return x+y}
void main(){
double x,y,x0,y0,h,m,k1,k2,k3,k4,z;
cout<<"enter x0,y0,h,m\n";
cin>>x0>>y0>>h>>m;
x=x0;
y=y0;
cout<<"\nx\ty";
cout<<"\n-------"<<'\n'<<x<<'t'<<y;
for(int i=1;i<=m;i++){
k1=h*f(x,y)
x=x+0.5*h;
z=y+0.5*k1;
k2=h*f(x,z);
z=y+0.5*h;
k3=h*f(x,z);
x=x+0.5*h;
z=y+k3;
k4=h*f(x,z);
y=y+(k1+2*k2+2*k3+k4)/6;
cout<<'\n'<<x<<'\t'<<y;
{
{