Numerical Mathematics - A Collection of Solved Problems

NUMERIČKA INTEGRACIJA 299kb=0h=gg-dga=(fun(dg)+fun(gg))/2pom=h*ado 50 k=1,15x=dg+h/210 a=a+fun(x)x=x+hif(x.lt.gg) go to 10t(k)=h/2*ab=1if(k.eq.1) go to 20do 15 m=1,k-1i=k-mb=4*b15 t(i)=(b*t(i+1)-t(i))/(b-1)20 b=4*bvint=(b*t(1)-pom)/(b-1)if(dabs(vint-pom).le.eps) returnpom=vint50 h=h/2kb=1endLista u potprogramu ima sledeće značenje:dg – donja granica integrala;gg – gornja granica integrala;fun – ime funkcijskog potprograma kojim se definiše podintegralna funkcija;eps – zahtevana tačnost izračunavanja;vint – vrednost integrala sa tačnošću eps, ukoliko je kb = 0;kb – kontrolni broj (kb = 0 integral je korektno izračunat; kb = 1 tačnostizračunavanja integrala nije postignuta sa 15 predvid¯enih koraka, tj.sa brojem podsegmenata 2 15 ).U naredbicommon navedena je zajednička promenljiva kojom se definiše parametaru podintegralnoj funkciji. U konkretnom slučaju, podintegralnu funkciju zaeliptički integral L(c) definišemo na sledeći način: