/*
 * @(#) smt2_3.c - Program for determining matrix of the battle effectiveness
 * (Task 2_3 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 20:20:00 GMT+04:00
 */

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

#define G 9.81

#define MAX_N 10
#define MAX_M 10

int main(void)
{
 int n,m;
 float zmax,r2;
 float q[MAX_M][MAX_N];
 float c[MAX_M],ng[MAX_M],v[MAX_M],bt[MAX_M];
 float p[MAX_N],bs[MAX_N];
 int i,j;
 float d,rmin,phi1s,phi2s,phi1t,phi2t;
 printf("n= ");
 scanf("%d",&n);
 printf("m= ");
 scanf("%d",&m);
 printf("z_max= ");
 scanf("%f",&zmax);
 printf("R2= ");
 scanf("%f",&r2);
 printf("\n");
 for (j=0;j<m;j++)
 {
  printf("C[%d]= ",j+1);
  scanf("%f",&c[j]);
  printf("Ng[%d]= ",j+1);
  scanf("%f",&ng[j]);
  printf("V[%d]= ",j+1);
  scanf("%f",&v[j]);
  printf("Betta_T[%d]= ",j+1);
  scanf("%f",&bt[j]);
 }
 printf("\n");
 for (i=0;i<n;i++)
 {
  printf("P[%d]= ",i+1);
  scanf("%f",&p[i]);
  printf("Betta_S[%d]= ",i+1);
  scanf("%f",&bs[i]);
 }
 printf("\n");
 for (j=0;j<m;j++)
 {
  rmin=v[j]*v[j]/(G*sqrt(ng[j]*ng[j]-1.0));
  phi1t=bt[j]-2.0*rmin/r2;
  phi2t=bt[j]+2.0*rmin/r2;
  d=4.0*rmin/r2;
  for (i=0;i<n;i++)
  {
   phi1s=bs[i]-zmax/r2;
   phi2s=bs[i]+zmax/r2;
   phi1s=phi1s>=phi1t ? phi1s : phi1t;
   phi2s=phi2s<=phi2t ? phi2s : phi2t;
   q[j][i]=phi2s-phi1s;
   q[j][i]=q[j][i]>0.0 ? q[j][i]/d : 0.0;
  }
 }
 for (j=0;j<m;j++)
  for (i=0;i<n;i++)
    q[j][i]*=c[j]*p[i];
 for (j=0;j<m;j++)
 {
  for (i=0;i<n;i++)
   printf("\t%3.3f",q[j][i]);
  printf("\n");
 }
 printf("\n");
 return 0;
}