Игла Бюффона на Си

[Bilymo]

Здравствуйте! Будьте любезны пожалуйста помочь. Я пишу программу, по определению числа пи, при условии x<=l*sin(θ), где 0<θ<pi, 0<x<l - эта проекция угла θ, а l - единичная длина иглы

#include <stdio.h>
#include <omp.h>
#include<conio.h>
#include<stdlib.h>
#include<time.h>
#include<math.h>
#include<iostream>
 
int main()
{
    int a,l,n,m;
double p,val,pi,x,fi,eps;
double mp;
x=1;
l=1;
srand(time(NULL));
pi = 4 * ((double)atan(1.0));
for(int i = 0; i < 10; ++i)
{
mp = (double)(rand()%22)/7;
printf("%f \n",mp);
}
eps=0.01;a=1;l=1;
printf("dsds \n");
do
{
fi=(double)(rand()%22)/7;
n=n+1;
if (x<=l*sin(fi))
{
m=m+1;
p=(2*n/m);
}
printf("steps %d pi: %e \n",m,p);
}
while(fabs(p-pi)<eps);
 
 
getche();
}

Диагноз: неверный результат, и есть ли библиотечная переменная может быть, в которой затабулировано исходное число пи? Ни где этого не нашёл

[BDA]

Число Пи есть в math.h - M_PI. Насчет неверного подсчета не помогу (лень разбираться в теории). Не забывайте оформлять код кнопочкой решетка #.

[Bilymo]

Сообщение от BDA

Число Пи есть в math.h - M_PI. Насчет неверного подсчета не помогу (лень разбираться в теории). Не забывайте оформлять код кнопочкой решетка #.

Проблема ещё заключается в том, что оператору rand()%M_PI нельзя вроде задать нецелый тип... а тут никакой заумной теории нету, условие должно быть правильно..., как я портировал с паскаля. Но вот проблема с rand....
Кстати, если rand()%M_PI - необъявленный индентификатор

[Bilymo]

Сообщение от Bilymo

Проблема ещё заключается в том, что оператору rand()%M_PI нельзя вроде задать нецелый тип... а тут никакой заумной теории нету, условие должно быть правильно..., как я портировал с паскаля. Но вот проблема с rand....
Кстати, если rand()%M_PI - необъявленный индентификатор

решил вообщем с http://stackoverflow.com/questions/1...ffons-needle-c