Math Is Fun Forum

Input y (y>1)(0 for quit):112349123746128387461279838741628934876129341237463924
The solution of x^x=1.12349e+053 is:34.4979.
Input y (y>1)(0 for quit):23452432324957346598724365234987623945623475238476239847623984762394862
348962349
The solution of x^x=2.34524e+079 is:47.3708.
Input y (y>1)(0 for quit):1102374123084712304821930487213498071234098732094712340123478901732081
The solution of x^x=1.10237e+069 is:42.4204.
Input y (y>1)(0 for quit):12039487123804981392847382905651208936410293874123098471238490382714809
23874019238741809238472189568293047213094
The solution of x^x=1.20395e+111 is:61.9786.
Input y (y>1)(0 for quit):12349812634906785610238470238109471324098465182394712380947238409734863
2840923174819028347120398416581372410938471234896358903471823409183409178236581293047128390472138
94
The solution of x^x=1.23498e+169 is:87.1489.
Input y (y>1)(0 for quit):19837469812374598123712894716328941723649127834537891263419283746192378
4126937845128374619234691287345617858912364198287346981273649182374612389476374573846132947123946
1273461923784519326493812874162398476397851374612384912347612398745713984612938764238914761937846
193246918
The solution of x^x=1.98375e+273 is:129.405.
Input y (y>1)(0 for quit):91236549123461289374612398451723469123874612987346123984762357894321034
7507236401823748941287346172389475197813264123784587346871346129847364982734512379847612348912734
6178235412783461298351364781923746129837643294872367461392841763489726349183247612839746328941723
64918237461298374612938476123795793649182387416923746127893476127983461932
The solution of x^x=1.#INF is:-1.#IND.
Input y (y>1)(0 for quit):

#include<iostream>
#include<math.h>
using namespace std;
int const pr=10;
double const e=2.7182818284590452354;
double const l2=0.69314718055994530942;
double omega(int n,double x);
int main(){
st: cout<<"Input y (y>1)(0 for quit):";
	double y;
	cin>>y;
	if(y==0.) goto en;
	cout<<"The solution of x^x="<<y<<" is:"<<pow(e,omega(pr,log(y)))<<".\n";
	goto st;
en: cout<<"Press any key to continue..."<<endl;
	char c;
	cin>>c;
	return 0;
}

double omega(int n,double x){
	if(n==0){
		return 0;
	}
	else{
		double om=omega(n-1,x);
		double omp=pow(e,om);
		double result=om-(om*omp-x)/(omp*(om+1)-((om+2)*(om*omp-x))/(2*om+2));
		return result;
	}
}