summaryrefslogtreecommitdiff
path: root/1752/CH3/EX3.10/exa3_10.sce
blob: 814df94bab8c8a0fef7235af10153d3a19bec07b (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
//Exa 3.10
clc;
clear;
close;
//given data
L=0.12;// in meter
t=.15*10^-2;// thickness in m
K=55.5;// in W/mK
h=23.5;// in W/mK
T_L=357;// in K
T_0=313;// in K

// Formula m=sqrt(h*rho/(K*A)) and rho=%pi*d and A=%pi*d*t, putting value of rho and A
m=sqrt(h/(K*t));
mL=m*L;
mL=floor(mL);
// Formula (T_L-T_infinite)/(T_0-T_infinite)= 1/cosh(m*L)
T_infinite=(T_L-T_0/cosh(mL))/(1-1/cosh(mL));
T_infinite=ceil(T_infinite);
measurement_error=T_infinite-T_L;
disp("Measurement Error is : "+string(measurement_error)+" K")

// Note: In the book, Unit of answer is wrong