由于给出的微分方程比较复杂（未知变量过多，如I、t、b、n、I、N、d），所以求解这个方程，必须把b、n、I、N、d的具体数值代入，其他还得给出其初始条件，这样才能解出I(t)的解析值。即I=dsolve('DI=。。。'，'I(0)=0')t=0: n %时间从0-nI=eval(I) %计算对应t的电流值plot(t,I) %绘制t-I的曲线图如按上述方法，仍解不出可以用ode（）函数求得其数值解。
