Section 10.2 Heat Conduction in a Cylinder with an Electrical Heat Source
Following the same steps as BS&L used.
>  restart;
>  A:=r-> 2*Pi*r*L;  area of conduction 
>  Q:=r-> A(r)*qr(r);  Heat conducted across cylinder of radius r in the direction of increasing r: qr is the flux of energy.
>  Q(r);  the rate of thermal energy into cylindrical shell at r: eq. 10.2-2 
>  Q(r + dr);  rate of thermal energy out across cyllindrical surface at r+dr: eq. 10.2-3
>  V:=r->A(r)*dr;  Volume of shell
>  Epro:=(r)->V(r)*Se;  Thermal production rate by electrical current: eq. 10.2-4
>  Epro(r);
>  LHS:=limit((Q(r+dr)-Q(r))/dr,dr=0);  Left hand side of eq. 10.2-5
>  de:=simplify(r*(LHS=Epro(r)/dr)/A(r)); energy balance, out - in = energy production by electrical dissipation giving 10.2-6
>  sol:=dsolve({de,qr(0)=finite},qr(r));  Solving 10.2-6 with BC 10.2-8 to get 10.2-9
>  assign(sol);  Giving qr(r) a value
>  qr:=unapply(qr(r),r);  Making qr a function of r
>  deFL:=qr(r)=-k*D(T)(r);  Fourier's Law
>  sol:=dsolve({deFL,T(R)=T0},T(r)); solving the de with the bc 10.2-12 to get 10.2-13
>  assign(sol);T:=unapply(T(r),r);
>  dT:=simplify(T(r)-T0);
>  dToverCon:=simplify((dT/(Se*R^2/(4*k))));dTmax:=simplify(T(0)-T0);dTavg:=simplify(int((T(r)-T0)*r,r=0...R)/int(r,r=0...R)); QatR:=Q(R);  Our answers agree with eqs. 10.2-13, 14, 15 and 16

>