MATLAB: Use built-in bvp4c to solve thermal problem

Question

I am writing some code while using the built in function bvp4c but I keep obtaining errors that I am not sure how to fix. Here is the code I have, any help would be greatly appreciated!

function res = bcfun(Ta,Tb)
BCa = 473; BCb = 273;
res = [Ta(1)- BCa
    Tb(1) - BCb];
end
function dydx = odefun(x,yw)
% dydx = odefun(x,T,yw)
hc = 40;
P = 0.016;
eps = 0.4;
k = 240;
Ac = 1.6E-5;
Seg = 5.67E-8; 
Ts = 293;
kAc = k*Ac;
A1 = hc*P/kAc;
A2 = eps*Seg*P/kAc;
dydx = A1*(T-Ts)+ A2*(T^4 - Ts^4);
end
function Problem3
solinit = bvpinit(linspace(0,0.1,20),[0.2 0.2]);
sol = bvp4c(@odefun, @bcfun, solinit);
plot(sol.x, sol.y(1,:),'r')
xlabel('x'); ylabel('y');
end

Answer 1

Problem3   %function call
function res = bcfun(Ta,Tb)
BCa = 473; BCb = 273;
res = [Ta(1)- BCa
    Tb(1) - BCb];
end
function dydx = odefun(x,Y)
% dydx = odefun(x,T,yw)
hc = 40;
P = 0.016;
eps = 0.4;
k = 240;
Ac = 1.6E-5;
Seg = 5.67E-8; 
Ts = 293;
kAc = k*Ac;
A1 = hc*P/kAc;
A2 = eps*Seg*P/kAc;
dydx = [Y(2);Y(1).^4.*9.450000000000002e-8+Y(1).*(5.0e2./3.0)-4.952980313402783e4];
end
function Problem3
solinit = bvpinit(linspace(0,0.1,20),[0.2 0.2]);
sol = bvp4c(@odefun, @bcfun, solinit);
plot(sol.x, sol.y(1,:),'r')
xlabel('x'); ylabel('y');
end