Temperature distribution in drag flow restart: with(detools): ode: =diff(T (y),y, y)=-2/ki de: ="T(y) k Forced (Dirichlet")boundary conditions >T f: =simplify(dsolve(fode, T(0)=0, T(1)=0), T(y)))i I Oy( Digits: =4: k: =1: 2: =1: eql: =rhs(T f): Natural (Cauchy") boundary conditions >T n:=simplify (dsolve((ode, T(0)=0), T(y)))i Tn=T(n bc n: =subs(y=l, diff(rhs (r n) y))=-subs(y=l, rhs(T n))i solve(subs(9=1, bc n),C1)i >c1:=3/4:eq2:=rhs(Tn): plot(eql, eq2), y=0.1, thickness=3)i 028 0.2 0.24 022 02 018 016 0.14 0.12 008 006 0.04 002 0 04 06 08