我正在MATLAB中为特定系统求解扩散方程，但是我得到的第二个解决方案nsub（r）存在问题。理想情况下，我希望它只是r的函数，具有两个未确定的常量，而表达式的第一部分确实是一些besselfunctions，它还包含r，'IgnoreAnalyticConstraints'，true）。有谁知道如何解决这个问题，出了什么问题？这是我的代码：
clear all
clc
%%%%%% constants %%%%%%%
r_nw=50*10^(-9);
pitch_nw=500*10^(-9);
f_sub=10^18;
f_nw=10^18;
L_nw=100*10^(-6);
%%%%%% diffusion on wire %%%%%
tau_nw=1000*10^(-9);
D_nw=1*10^-9;
lambda_nw=sqrt(D_nw*tau_nw);
%%%%%% diffusion on substrate %%%%%%%
tau_sub=100*10^(-9);
D_sub=1*10^(-9);
lambda_sub=sqrt(D_sub*tau_sub);
%%%%%% diffusion equation on wire %%%%%%%%
syms nw(z)
nw = dsolve(D_nw*diff(nw,z,2) == -f_nw + nw/tau_nw,nw(L_nw)==0);
nw(z)=simplify(nw);
dnw = diff(nw,z);
%%%%%% diffusion equation on substrate %%%%%%
syms nsub(r)
nsub = dsolve(diff(nsub,r,2)+1/r*diff(nsub,r) == -f_sub + nsub/tau_sub);
dnsub = diff(nsub,r);