帮忙matlab求两条曲线交点程序，不知问题出在哪里。

帮忙看一下两条曲线交点的问题。菜鸟+新手弄不明白。 % temperature diffusion 温度扩散曲线 z = 0:300; %depth in mbsf 深度 t1 = 500; % time b = 0.054; %thermal gradient DT = 2; %temperature change k = 3E-7; %thermal diffusivity To = 2.2; %initial temperature T1 = To+b.*z+DT.*erfc(z./2.*sqrt(k.*t1)); % 温度随深度扩散曲线 hold on; %plot temperature diffusion line plot(T1,z,'g'); %green %Xlabel('Temperature (C)'); %Ylabel('Depth (mbsf)'); % hydrate stability Line 稳定方程曲线 % 3.35% salt water 100% pure methane P = 0.101325+1030.*9.81.*(750+z).*0.000001; %净水压力 LP = log10(P); TB = 1./(0.00379-0.000283.*LP)-273.15; %稳定方程曲线 % plot hydrate equilibrium plot(TB,z,'k--'); %black dash %到这里都没问题 % junction point T z 求两条曲线交点，下面这行出错。 result=solve('T = To+b.*z+DT.*erfc(z./2.*sqrt(k.*t1))','T = 1./(0.00379-0.000283.*log10(0.101325+1030.*9.81.*(750+z).*0.000001))-273.15'); x = double(result,T); Y = double(result,z);