double s; a,b; double W(double x){ double nl=0,nr=3.4,m,t; int i; for(i=0;i<50;i++){ m=(nl+nr)/2; t=m*exp(m); t