Values on plot do not match manually inputted values for an anonymous function.

Question

0 个投票

This is the equation I am using. The two unknown variables are "t" and "T". My code is:

clear all;
clc;
syms T Tsur Ti r A V rho c SB Se t;
Tsur = 300;       % Kelvin, temp of chamber
Ti = 800;         % Kelvin, initial temp of ball
r = 0.006;        % Meters, radius of ball
A = 4*pi*r^2;     % Meters^2, surface area of ball
V = (4/3)*pi*r^3; % Meters^3, volume of ball
rho = 8000;       % kg/m^3, density of ball
c = 480;          % J/kg-K, specific heat of ball
SB = 5.6700e-08;     % W/m^2-K^4, Stefan-Boltzmann constant
Se = 0.1;         % Surface emmisivity
e = 1e-8;
m = 10000;
t = @(T) ((rho*V*c)/(4*Se*A*SB*(Tsur^3)))*((log(abs((Tsur+T)/(Tsur-T))) - (log(abs((Tsur+Ti)/(Tsur-Ti)))) + 2*(atan(T/Tsur)-atan(Ti/Tsur))));
T = 400:800;
plot(T,t(T))