Evaluating ordinary differential equations by RK-4 or classical method both order first and second

//To compute ordinary differential equations by RK-4 or classical method for order first

#include

#include

#include

void main()

{

cout<<“\t first order R-K 4th Classical Method\t”<

double h,n,m1,i,m2,m3,m4,x[100],y[100];

cout<<“Number of intervals “;

cin>>n;

cout<

cout<<“Enter the value of h”<<<“h = “;

cin>>h;

cout<<“Enter the value of x0″<<<“x0 = “;

cin>>x[0];

cout<<“Enter the value of y0″<<<“y0 = “;

cin>>y[0];

cout<

for(i=0;i<=n;i++)

{

x[i+1]=x[i]+h;

m1=h*.5*(1+x[i])*(y[i]*y[i]);

m2=h*.5*(1+(x[i]+.5*h))*(y[i]+m1*.5*h)*(y[i]+m1*.5*h);

m3=h*.5*(1+(x[i]+.5*h))*(y[i]+m2*.5*h)*(y[i]+m2*.5*h);

m4=h*.5*(1+(x[i]+h))*(y[i]+m3*h)*(y[i]+m3*h);

y[i+1]=y[i]+(m1+2.0*m2+2.0*m3+m4)/6;

cout<<“x”<<<“= “<<

cout<<“y”<<<“= “<<

cout<

}

getch();

}

// To compute ordinary differential equations by RK-4 or classical method for order second

#include

#include

#include

void main()

{

cout<<“\t 2nd order order R-K 4th Order Method\t”<

double h,n,m1,i,m2,m3,m4,p1,p2,p3,p4,x[100],y[100],z[100];

cout<<“Number of intervals “;

cin>>n;

cout<

cout<<“Enter the value of h”<<<“h = “;

cin>>h;

cout<<“Enter the value of x0″<<<“x0 = “;

cin>>x[0];

cout<<“Enter the value of y0″<<<“y0 = “;

cin>>y[0];

cout<<“Enter the value of z0″<<<“z0 = “;

cin>>z[0];

cout<

for(i=0;i<=n;i++)

{

x[i+1]=x[i]+h;

m1=h*z[i];

p1=h*(6*x[i]+3*y[i]-2*z[i]);

m2=h*(z[i]+.5*p1);

p2=h*(6*(x[i]+.5*h)+3*(y[i]+m1*.5)-2*(z[i]+.5*p1)) ;

m3=h*(z[i]+.5*p2);

p3=h*(6*(x[i]+.5*h)+3*(y[i]+m2*.5)-2*(z[i]+.5*p2)) ;

m4=h*(z[i]+.5*p3);

p4=h*(6*(x[i]+.5*h)+3*(y[i]+m3*.5)-2*(z[i]+.5*p3)) ;

y[i+1]=y[i]+(m1+2.0*m2+2.0*m3+m4)/6;

z[i+1]=z[i]+(p1+2.0*p2+2.0*p3+p4)/6;

cout<<“x”<<<“= “<<

cout<<“y”<<<“= “<<

cout<<“z”<<<“= “<<

cout<

}

getch();

}

To evaluate ordinary differential equations by RK-4 (Classical) Method in C++ Programming Language for Numerical Methods for Engineering Students