Nusselt number and skin friction for bvp4c

Hi,

I need to include the nusselt number and skin friction in the program.

How should I include these in program. Please help me in this regard

Program is given below

close all

clc

phi = 0.01;

Betaf = 207;

Betas = 17;

Beta = 0.5;

kof = 0.613;

kos = 400;

m = -1;

bk= 0.5;

bv = 0.5;

GRT = 10;

GRC = 10;

Alpha = 0.5;

Rhof = 997.1;

Rhos = 8933;

P = 0.5;

n = 0.3;

Br = 0.01;

A1 = (1-phi).^2.5;

A2 = 1/(1 + 1/Beta);

A3 = (1-phi)+phi.*((Rhos.*Betas)./(Rhof.*Betaf));

G1 = A2.*A3;

G3 = A1.*A2;

G4 = (kos + 2*kof - 2*phi.*(kof - kos))/(kos + 2*kof + phi.*(kof - kos));

G5 = (1 + 1/Beta)./A1;

G6 = G5./G4;

dydx=@(x,y)[y(4);

y(5);

y(6);

bv.*y(4).*y(5)-(1+bv.*y(2)).*(GRT.*G1.*y(2)+GRC.*G1.*y(3)-G3.*P);

bk.*y(5).^2-Br.*G6.*y(4).^2-(bk-bv).*y(2).*Br.*G6.*y(4).^2+bk.*bv.*Br.*G6.*y(2).^2*y(4).^2;

Alpha.*y(3)];

BC = @(ya,yb)[ya(1);yb(1);ya(2)-(1+m);yb(2)-1.0;ya(3)-(1+n);yb(3)-1.0];

yinit = [0.01;0.01;0.01;0.01;0.01;0.01];

solint = bvpinit(linspace(-1,1,11),yinit);

S1 = bvp4c(dydx,BC,solint);

Answer 1

[Torsten]

Nu1 = -(1+bk*(1+m))*S1.y(5,1)
Nu2 = -(1+bk)*S1.y(5,end)
Tau1 = (1-bv*(1+m))*S1.y(4,1)
Tau2 = (1-bv)*S1.y(4,end)

[Syed Mohiuddin]

Thank you very much for the answer