long double から減算を行うたびに、間違った答えが返されます
long double Runge(long double a, long double b, int m)
{
long double R;
R=(b-a)/((2^m)-1);
cout<<setprecision(16)<<R<<endl;
return(R);
}
たとえば、a=0.53733 b=0.53736 m=2 の場合、R は -2.805572983821669e-006 に等しくなります
なぜこれが起こるのですか?