area of 10m2 .overa 0.05 − −t 0.000018t^2

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);
Error in dTdt (line 2)

clc
clear
[t,T]=ode45(@dtdT,[0 3600],(40));
plot(t,T,'-')
xlabel('Time (s)')
ylabel('Temperature (oC)')

Error in timetemp (line 3)
[t,T]=ode45(@dtdT,[0 3600],(40));

Answer 1

[Chunru]

% 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