/*
 * @(#) smt1_5.c - Program for calculation of the hitting factor of
 * some bombs, variant 1 (Task 1_5 at the Special Military Training course).
 * (c) 1996 Ivan Maidanski <ivmai@chat.ru> http://ivmai.chat.ru
 * Freeware program source. All rights reserved.
 **
 * Language: ANSI C
 * Tested with: Borland C++ v3.1
 * Last modified: 1996-12-06 15:40:00 GMT+03:00
 */

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <time.h>

#define KSI_N 100

double norm_rnd_func(double mx, double sigmax)
{
 int i;
 long ksi=0;
 double z;
 for (i=0;i<KSI_N;i++)
  ksi+=rand();
 z=(2.0*ksi/RAND_MAX-KSI_N)*sqrt(3.0/KSI_N);
 return ((z-(3.0+(1.0-z*z)*z)/KSI_N/20)*sigmax+mx);
}

int main()
{
 int k,n,m,l,r;
 int l1,l2,i,j;
 float x,y,sigma_v,sigma_u,sigma_x,sigma_y;
 randomize();
 printf("Input N:\n");
 scanf("%d",&k);
 printf("Input L:\n");
 scanf("%d",&l);
 printf("Input R:\n");
 scanf("%d",&r);
 printf("Input n:\n");
 scanf("%d",&n);
 printf("Input m:\n");
 scanf("%d",&m);
 printf("Input Sigma_v:\n");
 scanf("%f",&sigma_v);
 printf("Input Sigma_u:\n");
 scanf("%f",&sigma_u);
 printf("Input Sigma_x:\n");
 scanf("%f",&sigma_x);
 printf("Input Sigma_y:\n");
 scanf("%f",&sigma_y);
 l1=0;
 for (i=0;i<k;i++)
 {
  l2=0;
  x=norm_rnd_func(0,sigma_x);
  y=norm_rnd_func(0,sigma_y);
  for (j=0;j<n;j++)
   if (fabs(norm_rnd_func(x,sigma_v))*2<l &&
       fabs(norm_rnd_func(y,sigma_u))<r && (++l2)==m)
    break;
  if (l2==m)
   l1++;
 }
 printf(" W=%f\n\n",(float)l1/k);
 return 0;
}