#include <iostream>
#include <cmath>
#include <iomanip>

using namespace std;

#define sqr(x) ((x)*(x))

int main(){

  double Deltatau,tau;
  double vx_0[4], vy_0[4], vx_1[4], vy_1[4];
  double E;
  int N, n;

  cin >> Deltatau >> N;
  
  int k;

  //Entrer les conditions initiales
  //==================================
  for(n=0;n<4;n++){
    cin >> vx_0[n] >> vy_0[n];
  }
  //==================================


  //Afficher les conditions initiales 
  //=========================================================================== 

  tau = 0.0;

  cout << setprecision(9) << tau << " " << fabs(vx_0[0]*vy_0[1]-vx_0[1]*vy_0[0]) << endl; 
  
  for(k=0;k<4;k++){
    E = 0.5*( sqr(vx_0[k]) + sqr(vy_0[k]) );
    //    cout << setprecision(9) << tau << " " << vx_0[k] << " "<< vy_0[k]<< " "<< E;
  }
  //cout << endl;
  //===========================================================================

  //Evolution dans le temps 

  for(n=1;n<=N;n++){

    tau = n*Deltatau;

    for(k=0;k<4;k++){
      //Euler symplectique 
      //===========================================================
	vx_1[k] = vx_0[k] + Deltatau*vy_0[k];
	vy_1[k] = -Deltatau*vx_0[k] + (1.0-sqr(Deltatau))*vy_0[k];
      //=========================================================

	E = 0.5*( sqr(vx_1[k]) + sqr(vy_1[k]) );

	//	cout << setprecision(9) << tau << " " << vx_1[k] << " "<< vy_1[k]<< " "<< E;


      } 

    //    cout << endl;


    //    cout << setprecision(9) << tau << " " << fabs(vx_0[0]*vy_0[1]-vx_0[1]*vy_0[0]) << endl;

    //Mettre à jour les vitesses
    for(k=0;k<4;k++){
	vx_0[k] = vx_1[k]; vy_0[k] = vy_1[k];
    }

    cout << setprecision(9) << tau << " " << fabs(vx_0[0]*vy_0[1]-vx_0[1]*vy_0[0]) << endl;

  } //L'évolution arrive à sa fin

  return(0);
} //Fermer la fonction main