- #1

MathematicalPhysicist

Gold Member

- 4,677

- 362

can anyone spot my error in my code, i appreciate:

Code:

```
#include <stdio.h>
#include <math.h>
main()
{
long double initial,area,Segments,dx;
initial=-2.5;
area=0;
printf("insert the number of sgements to subdivide the interval \n");
scanf("%f", &Segments);
if(Segments!=0) dx=4.5/Segments;
else printf("error \n");
while(area<=84.515625 && initial<2)
{
area+=0.5*dx*((initial+dx)*(initial+dx)*10-(initial+dx)*(initial+dx)*(initial+dx)+initial*initial*10-initial*initial*initial);
if((fabs(area-84.515625))<0.001){
printf("the numerical integration yields %f with accuracy of 0.001 from analytical result \n",area);
break;
}
initial+=dx;
}
}
```