W dniu 2012-11-15 11:23, AK pisze:
>
>> Oczywiście wszystko się zmienia, jeżeli mnożenie przez wagi jest pod
>> pętlą (a tak jest np. w implementacji Simpsona z Wiki -
>
> Chyba tylko jakis _kompletny nieuk_ mnozylby nie wiadomo po co pod petla
> :)).

Ojtam, ja nie tylko mnożyłem pod pętlą, ale i budowałem
wielkie wektory tych wag - do zastosowań bojowych zupełnie
bez sensu;)

Roman W w swoim rozwiązaniu wstawił mnożenie pod pętlą.
Jak całkujesz skomplikowaną funkcję zawierajacą
takie kwiatki: http://mathworld.wolfram.com/Polylogarithm.html
to dodatkowe mnożenie niewiele zmienia.


Dziwi mnie ta gałąź. Sławek się burzy, a wjego kodzie nie
było mnożenia pod pętlą:

double integrate_simpson(double x[], double y[], int n)
{
int i;
double odd = 0., even = 0.;

for(i = 1 ; i <= n ; i+=2) odd += y[i];
for(i = 2 ; i <= n ; i+=2) even += y[i];

return (2.0*odd-y[1]-y[n] + 4.0*even)*(x[3]-x[1])/6.0;
}


Jak ktoś nie wierzy kompilatorowi, można wprost mu kazać:

//n parzyste!
double even=-0.5*y[0];
double odd = 0.0;

for(i = 1 ; i <= n ; i+=2)
{
odd += y[i];
even += y[i-1]
}
return (2.0*even + 4.0*odd + y[n])*stale

Albo od razu
if (n%2==0) return (2.0*even + 4.0*odd + y[n])*stale;
else return (2.0*even + 4.0*odd + a*y[n] + b*y[n-1] + c*[n-2] )*stale;

gdzie a,b,c to współczynniki, o których pisałem wczoraj.


pzdr
bartekltg