a5083c82b50a0140

d:\MEiL\Semestr I1\Informatyka I7\Wykład\... treningowa :i (J. Rokicki) Uadan i e [2j.cpp 2

tldefino MAXN 10 // maksymalna ilosc rownar:

// wykonuje jeden krok całkowania

II wektorowego równania rozniczkowego zwyczjanego:

H dY/dx - Fun(x,Y), Y(xOi-YO // metoda Rungego-Kutty lv-r.ego rzędu.

// Parametry formalne:

// x0 - wartość startowa zm. niezależnej

// yO - wartość startowa zm. zalozncj {tablica n-eiementowa)

// n    liczba równań

// funlx,y,prawastr> - nazwa funkcji obliczającej prawe strony // yl - obliczona wartość zmiennej zależnej w punkcie xO+h //    (tablica n-eiementowa)

void vrk4 ( double x0, double yO 11, double h, mt n, void (‘funt :<ioui>le, doubie*, double*-, double ylI) )

double ki [MAXN], k2 [HAXN], k3 JHAXN |, k-1 |MAXN]; double ytmp[MAXH];

fur. ( x0, yO, kl>; for { i»0; Kn; ł+i)

I

kl[i] * h;

ytmpti 1 yO[i) ♦ kl[i)/2.0;

fun( x0+h/2.0, ytmp, k2); for ( i«=0; i<n; •-'•i)

k2(i| *= h;

ytmp(i) - yOliJ * Ic2 (11/2.0;

fun{ x0*h/2.0, ytmp, k3}; for ( i=0; i<n; ++i)

k3(l] *= h;

ytmp[ij = yO(il *■ k3[i];

fun( xO<-h, ytmp, k4;; for { i=0; i<n; • k4(ij — h;

; i-0; i<n; **i)

yim ** yom + (kl (i] + 2. *k21 i ] * 2.*k3[il * k4[l|l/6.;