roz1

d:\MElL\Semestr iIMnforraatyka :i\Wykład\. remngowe (J. Rokicki i\Zadanie [l;-‘~pp

//Instrukcje preprocesora ♦include <st.dio.h>

•include <math.h>

//Deklaracje funkcji

void wektorF{double x,double yU,double :||>, void vrk4( double x0, double yO|j, double b, double yl() );

//Funkcja mam void mair.O I

double y{4],yl[4;,x,xk; double h-0.01;

oid i'fum 'double, double*, double*).

x -0; xk*l0;

yioi-i *

ylll=2

yUI-i

yr31-2

8{yI0])+abs(yl1])>20 .

nia obliczeń

vrk4<x,y,h,4,

yI01-y!!0]; ylll-yim; y[21=yll2J; y 131 • yi 131;

if<i%10=»0)I

PJ*y(3n;

1

<ł wektorF(double x,double y[],double r101 exp(-x)-y[l]*y[0]-y[l]; f{l)=yl2J;

f[2l—y[3]*yt3)-y[0l*_V!2if

f(31-y[01*y(01-y[31ty[llj

ilosc rownan

wnania rozniczkowoąo zwyc2janego: // dY/dx - FunIX,Y),    Y(x0)-Y0

// metoda Rungego-Kut.ty !V-tego rzędu.

łdcfinc MAXN 10 // maksymaln. // wykonuje Jeden krok <

// wektorowego

Par.

>eLry formalne: wartość startów

er.aleznej

i-eiementowa)

x,y,prawastri - nazwa funkcji obliczającej prawe strony obliczona wartość zmiennej zależnej w punkcie x0+h (tablica rs-elementowa!

yO[|, double i

void (»fun: (double, double*^-, double*).

double klIMMCNI,

: |MAXN), k3(KAXN],

[MAXN1;