clear all
f=@(x) x*sin(x)-1
fplot(f,[0 12])
fzero(f,2)
fzero(f,4)
fzero(f,6)
fzero(f,9)
disp('Press any key to continue')
pause
clc
Ts=-15; Ti=20;
T=0; alpha=1.38e-7;
t=60*24*3600;
f=@(x) (T-Ts)/(Ti-Ts)-erf(x/2/sqrt(alpha*t))
fplot(f,[0 2])
fzero(f,1)
pause
clc
L=4.2
f=@(beta) 1+cosh(beta*L)*cos(beta*L) 
fplot(f,[0 2])
r=fzero(f,0.5);
omega=r^2*sqrt(21000/0.53)
r=fzero(f,1.5);
omega=r^2*sqrt(21000/0.53)
r=fzero(f,2);
omega=r^2*sqrt(21000/0.53)