Help me to reduce its accuracy?

I am solving entropy generation and heat transfer problem;

whose code is:

function Needle
a=0.01;
r=2.0;
Pr=3.0;
Ec=0.3;
Nr=10;
solinit= bvpinit(linspace(a,10,2500),[1 0 0 0 1]);
sol = bvp4c(@Needleode,@Needlebc,solinit);
eta = sol.x;
f = sol.y;

disp(f(3,1))
Cf=8*a^(1/2)*f(3); 
disp(Cf)

nusselt=-2*((a)^(1/2))*f(5)*(1+(4/(3*Nr))*(r)^3);
disp(nusselt)
plot(eta,f(4,:))
hold on

function dfdeta = Needleode(eta,f)
r=2.0;
Pr=3.0;
Nr=10;
Ec=0.3;
u=f(4)*(r-1)+1;
v=eta*(1+(4/(3*Nr))*(u)^3);
ff1= -(1/(2*eta))*(f(1)*f(3)+2*f(3));
t1=f(5)+(1/2)*(Pr*f(1)*f(5));
t2=4*eta*Ec*Pr*f(3)^(2);
t3=(2/(3*Nr))*(u^3)*f(5);
t4=(4/Nr)*(u^2)*eta*(f(5)^2)*(r-1);

ff2=-(1/v)*(t1+t2+t3+t4);
dfdeta=[f(2);f(3);ff1;f(5);ff2];


function res = Needlebc(f0,finf)
a=0.01;e=0.3;
res = [f0(1)-(a*e/2)
    f0(2)-e/2
    finf(2)-(1-e)/2
    f0(4)-1
    finf(4)-0
    
    ];

I am getting this error:

Warning: Unable to meet the tolerance without using more than 2000 mesh points.

The last mesh of 2500 points and the solution are available in the output argument.

The maximum residual is 0.0724122, while requested accuracy is 0.001.

> In bvp4c (line 323)

In Needle (line 8)

3.3796

2.7037

4.1376

I have to reduce value of Cf from 2.7037 to 2.6621 and value of nusselt nmbr from 4.1376 to 4.0623to obtain my desired results. suggest me suitable mesh points and initial guess.