#include <iostream>
#include <cstdlib>
#include <cmath>
#include <ctime>
using namespace std;

//Function name: Rand_Float
//This function will create a random mumber between 0 and 1.
//The input is a seed value
//The output is a number
double	Rand_Float(double a, double b)
{
	return ((double)rand()/RAND_MAX)*(b-a)+a;
}

int main()
{
	//Declare and Initialize

	int n(5000), k;
	double component_reliability(.95);
	double series_success(0), parallel_success(0);
	double num_1, num_2, num_3, num_4;
	double total, percent;
	unsigned int seed;

	// Get seed value from user
	cout<<"Please enter a seed value. ";
	cin>>seed;
	cout<<endl;
	srand (seed);


	//Determine the simulation reliability
	for (k=1; k<=n; k++)
	{
		num_1=Rand_Float(0, 1);
		num_2=Rand_Float(0, 1);
		num_3=Rand_Float(0, 1);
		num_4=Rand_Float(0, 1);

		if (((num_1<=component_reliability) &&
			(num_2<=component_reliability)))
		{
			series_success++;
		}
		else
		{
			if(((num_3<=component_reliability) ||
			(num_4<=component_reliability)))
		{
			parallel_success++;
		}
	}
	}

	//Determine the success as a percent
	total=series_success+parallel_success;
	percent=(total/n)*100;


		//Output results
	/*	cout<<"Simulation Riliability "<<endl;
		cout<<"Series: " <<series_success<<endl;
		cout<<"Parallel: " <<parallel_success<<endl;
		cout<<"Total: " <<total<<endl;                   */
		cout<<"Reliability: "<<percent<<" %"<<endl;


		//Exit Program
		return 0;

	}