“…specific heat of fluid and nanoparticles 1); dy(2,1) = y(3,1); dy(3,1) = (y(2,1)^2-y(1,1)*y(3,1)-K*y(5,1) + (K1)*y(2,1)-lmda*(y(6,1)-Nr*y(8,1)-Nc*y(10,1))/ (1 + K)); dy(4,1) = y(5,1); dy(5,1) = (y(2,1)*y(4,1) + K*(2*y(4,1) + y(3,1))-y(1,1)*y(5,1))/(1 + 0.5*K); dy(6,1) = y(7,1); dy(8,1) = y(9,1); dy(7,1) = -((eps + 4*Rd*(theta-1)*(1 + (theta-1)*y(6,1))^2)*y(7,1)^2 + Pr*y(1,1)*y(7,1) + Pr* (Nb*y(7,1)*y(9,1) + Nt*y(7,1)^2))/(1 + eps*y(6,1) + 1.3333*Rd*(1 + (theta-1)*y(6,1))^3); dy(9,1) = -(Nt./Nb)*dy(7,1)-Pr*Le*y(1,1)*y(9,1) + Pr*Le*sigma*y(8,1)*((1 + dlta*y(6,1))^n)* (exp(-E./(1 + dlta*y(6,1)))); dy(10,1) = y(11,1); dy(11,1) = -Lb*(y(1,1)*y(11,1)) + Pe*(y(9,1)*y(11,1) + dy(9,1)*(y(10,1) + dta1)); end Boundary conditions: global m gmaNbNt S res = [ya(1)-S, ya(2)-1, ya(4) + m*ya(3), ya(7) + gma*(1-ya (6)), Nb*ya(9) + Nt*ya (7),ya (10) (6);Nr = vec (7);Rd = vec (8);theta = vec(9);Le = vec(10); sigma = vec(11);dlta = vec(12);eps = vec(13);n = vec (14);E = vec(15);m = vec(16);gma = vec (17);lmda = vec(18);Lb = vec (19);…”