How I plot qSol(r) wrt r and TSol(r) wrt r ?

`%heat conduction with an electrical heat source in copper wire`

clc; clear;close all;%parameters

R = 0.2; % radius in m

L = 10; % length of wire in m

ke = 60000000; % electrical conductivity of copper in S/m

k = 401; % thermal conductivity of copper in W/m.K

I = 10; % electrical current in Amp

Se = (I*I)/ke; % rate of heat production per unit volume

T0 = 20 % Temp in Celcius

%Energy Balance in cylindrical shell of thickness del(r) and length L

% Doing the energy balance and taking del(r) common and limiting it to 0, a

% linear differential equation is formed.

% SOLVING DIFFERENTIAL EQUATION

syms q(r) %q is the heat flux

ode = diff(r*q,r) == Se*rcond = q(0)==0 %B.C at r = 0,q~=inf gives a expression which states that at r = 0 , q = 0

qSol(r)=dsolve(ode,cond)%Temperature Profile

syms T(r)ode = diff(T,r) ==(-Se*r/(2*k))cond = T(R)==T0;TSol(r)=dsolve(ode,cond)%Plot Flux profile

## Best Answer