//	Übung 5 - Aufgabe A/B  (dynamische Arrays, Funtionen)
//	g++ -Wno-deprecated u5AB.cc vectors.cc
//	icc -ansi -w2 -Wall u5AB.cc vectors.cc

//		45min von u4A  nach u5AB	[GH]

#include <iostream.h>
#include "vectors.hh"

int main()
{
 double *u, *w, *z;		// Pointer fuer dynamische Arrays
    int n,			// Laenge von u und w
        i;
 double u_norm, w_norm, uw_prod,// Normen und inneres Produkt
        alpha, beta;		// skalare Groesse

//			Datenberechnung und Speicherallokierung
 cout << " Eingabe n = "; cin >> n;

//			Speicher fuer u, w, z
 u = new double [n];
 w = new double [n];
 z = new double [2*n];

//			Berechnung u,w
 for ( i=0; i<n; i++ )
   {
     u[i] = fu( i );
     w[i] = fw( n, i );
   }

//		Berechnung der Normen und des Produktes
 uw_prod = InnerProd( n, u, w );
 u_norm  = EuklidNorm( n, u );
 w_norm  = EuklidNorm( n, w );

//		Berechnung alpha und beta
 alpha = uw_prod/u_norm;
 beta  = uw_prod/w_norm;

//		Berechnung von z
 Calc_z( n, z, alpha, beta, u, w );

//		Datenausgabe
 Print( 2*n, z );

 cout << " alpha = " << alpha << "    beta = " << beta << endl << endl;

 cout << " || u || = " << u_norm << "    || w || = " << w_norm
      << "    < u,w > = " << uw_prod << endl;

//		Speicherfreigabe
 delete [] z;
 delete [] w;
 delete [] u;

}