/*
 * @(#) smt2_4.c - Program for calculation of effectiveness of the radar
 * station as a system of queuing by solving a differential equation
 * (Task 2_4 at the Special Military Training course).
 * (c) 1997 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: 1997-05-04 21:15:00 GMT+04:00
 */

#include <stdio.h>

#define MAX_N 10
#define EPS 0.000001

int main(void)
{
 int n,j;
 float lambda,mu,h;
 float p[MAX_N+1],ps[MAX_N+1];
 float pp[MAX_N+1];
 int i;
 float s;
 printf("n= ");
 scanf("%d",&n);
 printf("h= ");
 scanf("%f",&h);
 printf("Lambda= ");
 scanf("%f",&lambda);
 printf("ServiceTime= ");
 scanf("%f",&mu);
 mu=1/mu;
 printf("\n");
 ps[0]=1.0;
 for (i=n;i>0;i--)
  ps[0]=ps[0]*lambda/(mu*i)+1.0;
 ps[0]=1.0/ps[0];
 for (i=1;i<=n;i++)
  ps[i]=ps[i-1]*lambda/(mu*i);
 p[0]=1.0;
 for (i=1;i<=n;i++)
  p[i]=0.0;
 for (j=1;;j++)
 {
  for (i=0;i<=n;i++)
   pp[i]=p[i];
  p[0]+=(-pp[0]*lambda+pp[1]*mu)*h;
  for (i=1;i<n;i++)
   p[i]+=(pp[i-1]*lambda-pp[i]*(lambda+mu*i)+pp[i+1]*mu*(i+1))*h;
  p[n]+=(pp[n-1]*lambda-pp[n]*mu*n)*h;
  for (i=0,s=0.0;i<=n;i++)
   s+=(ps[i]-p[i])*(ps[i]-p[i]);
  if (s<(n+1)*EPS)
   break;
 }
 printf("t= %f\n",h*j);
 for (i=0;i<=n;i++)
  printf("P%d= %f\n",i,p[i]);
 printf("\n");
 return 0;
}