Coding ODE metode runge kutta orde 4 dengan bahasa c

Bismillahirrahmanirrahim

Alhamdulillah, washalatu wasalamu 'ala rasulillah, waba'du

Berikut ini saya berikan contoh program untuk menghitung ODE metode range kutta orde 4 dengan bahasa c.

======================================================================================================================

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

int main()

{

int i=0;

double xnew=0, ynew=0, X=5, xi, yi, h=0.1;

double k1, k2, k3, k4;

xi=0, yi=4;

printf("i = %i, y(x=%.5f) = %.5f \n",i,xi,yi);

i=1;

while(xi<(X-h))

{

k1=2*xi-3;

k2=2*(xi+0.5*h)-3;

k3=2*(xi+0.5*h)-3;

k4=2*(xi+h)-3;

ynew=yi+(k1+2*k2+2*k3+k4)*(h/6);

xnew=xi+h;

printf("i = %i, y(x=%.5f) = %.5f \n",i,xnew,ynew);

xi=xnew;

yi=ynew;

i++;

}

return 0;

}

======================================================================================================================

Selamat mencoba!