Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- ```{r}
- #C1
- parabArea = function(N){
- counter = 0;
- a = 2;
- for (i in 1:N){
- u = runif(1, -sqrt(a), sqrt(a))
- v = runif(1, -sqrt(a), sqrt(a))
- w = runif(1, -a, a)
- realW = u*u + v*v;
- if (w > realW){
- counter = counter + 1;
- }
- }
- areaCuboid = (2*sqrt(a))**2 * 2*a
- areaParab = areaCuboid*counter/N;
- realArea = pi*a*a/2;
- absoluteError = abs(areaParab-realArea);
- relativeError = abs(areaParab-realArea) / abs(realArea);
- print("\n");
- print(absoluteError);
- print(relativeError);
- return (areaParab);
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement