MATLAB: Trouble finding solution of unknown variable in difficult integration problem

Question

I am trying to solve for 'h' using the first equation above. My script produces a result, however it is unexpected. The solution of 'h' should increase with 's', yet this is not the case. I have attached the script below, any help would be much appreciated.

clc; clear all; close all;
%%%%%%%%%% Physical Parameters %%%%%%%%%%
r=0.1;  
b=0.15;
W=60;
Kc=0.99*1000;    
Kphi=1528.4*1000;   
n=1.1;   
c=1.04*1000;    
phi=28; 
A1=(Kc/b+Kphi);
k=0.6;
Beta=0.0872665;
kx=0.043*Beta+0.036;
%%%%%%%%%% Defining equations %%%%%%%%%%
s=1
    
syms h 
syms Pheta
PhetaF=acos(1-h/r);
PhetaR=acos(1-k*h/r);
jx=r*(PhetaF-Pheta-(1-s)*(sin(PhetaF)-sin(Pheta)));
Sig=A1*r^n*(cos(Pheta)-cos(PhetaF));
Taux=(c+Sig*tand(phi))*(1-exp(-jx/kx));
%%%%%%%%%% Solving for h %%%%%%%%%%
fun=Taux*sin(Pheta)+Sig*cos(Pheta);
eqn=W==r*b*vpaintegral(fun,Pheta,[PhetaR,PhetaF]);
sol=vpasolve(eqn,h,.1)

Answer 1

I used a different approach - see below. I could only see an increase in h, with increasing s, as a step change around s = 4.5.

s = 1:0.1:10;
h0 = 0.1;   % Use h0 = 0.15 to get a different set of results.
for i = 1:numel(s)
    h(i) = fzero(@(h) Zfun(h,s(i)), h0);
end
plot(s,h,'o'),grid
xlabel('s'),ylabel('h')
legend(['initial guess h0 = ' num2str(h0)])
disp(h)
function Z = Zfun(h,s)
        r=0.1;  
        n=1.1;  
        rn = r^n;
        b=0.15;
        W=60;
        Kc=0.99*1000;    
        Kphi=1528.4*1000;
        k = 0.6;
        c=1.04*1000;    
        tanphi=tand(28); 
        Beta=0.0872665;
        kx=0.043*Beta+0.036;
        sigA = (Kc/b+Kphi)*rn;
        thetar = acos(1-k*h/r);
        thetaf = acos(1-h/r);
        
        sig = @(theta) sigA*(cos(theta) - cos(thetaf));
        jx = @(theta) r*(thetaf - theta - (1-s)*(sin(thetaf) - sin(theta)));
        taux = @(theta) (c + sig(theta).*tanphi).*(1 - exp(-jx(theta)/kx));
        kern = @(theta) taux(theta).*sin(theta) + sig(theta).*cos(theta);
        Z = r*b*integral(kern,thetar,thetaf) - W;
        
end

This results in