% The function name is "f" not "dtdT"
[t,T]=ode45(@f,[0 3600],(40));
plot(t,T,'-')
xlabel('Time (s)')
ylabel('Temperature (oC)')
function dTdt=f(t,T)
dTdt=(10000-0.05*t-0.000018*t^2)-(1010)*(T-(25-0.0003*t))+(0.20*4180)*(40-T)/(250*4180);
end
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!