I have a simple programg=1 ;
b= 1;
s= 1;
R= 2;
syms r;
func = @(r) 1 -(b/r)^2 - g^-2*((2/15)*(s/R)^9*(1/(r-1)^9 - 1/(r+1)^9 -9/8*r*(1/(r-1)^8- 1/(r+1)^8))-(s/R)^3*(1/(r-1)^3 - 1/(r+1)^3- 3/2*r*(1/(r-1)^2-1/(r+1)^2)));
minrootgbsR = fsolve(func,1.1)
k=1;
for i=1:length(minrootgbsR);
if(isreal(minrootgbsR(i))==1);
vec(k)=minrootgbsR(i);
k=k+1;
end
endvec;
max(vec)
clear all
The result of max(vec) must be 1.352104649412185 , and matlab give me1.3576
Why ? any ideas ? i try to use the multiprecision toolbox but dont work