شبیه سازی معادلات توزیع دما در راکتور هسته ای در MATLAB
global n delta rspan options r T
n = 100;
delta = 1.e-3;
rspan=linspace(delta,1,n+1);
options = odeset('RelTol',1e-12,'AbsTol',1e-12,'NormControl','on','MaxOrder',5);
T0=newton_mfile(@cp6_BV,2,1.e-10,0);
r=[0;r];
T=[[T0,0];T]...