Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- // calcular numericamente int_{0.1}^{1.0} (dx/x) = -ln(0.1)
- #include <stdio.h>
- #include <math.h>
- #define n 1000
- int main (void)
- {
- int i;
- double h,x,soma=0.0f; // esse f não dá problema
- h=(1.0-0.1)/n;
- for (i=0;i<n;i++)
- {
- x=0.1+h*i;
- soma+=h*(1/x);
- }
- printf("\nIntegral calculada: %f. Valor real: %f\n\n",soma, -log(0.1));
- return 0;
- }
- // o valor da integral foi estimado acima do seu valor real.
- // int_{x_{i}}^{x_{i+1}} f(x)dx = h*f_{i+1}
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement