Homework answers / question archive / Eigenvalues without pain: CHEBFUN code, Eigenfunctions can be derived using the CHEBFUN with MATLAB since it has an overloaded eig function

Eigenvalues without pain: CHEBFUN code, Eigenfunctions can be derived using the CHEBFUN with MATLAB since it has an overloaded eig function

Management

Share With

Eigenvalues without pain: CHEBFUN code, Eigenfunctions can be derived using the CHEBFUN with MATLAB since it has an overloaded eig function. The following code solves for the eigenfunctions of the Robin problem shown in the text in Section 11.3. The code can also readily evaluate the series coefficient. Thus the whole procedure of the method of separation of variables can be automated with a code similar to this. You will be able to solve or verify many of the problems in this chapter with this code and also will be able to use it in your research. xi = chebfun(’xi’,[0,1]); A = chebop(0,1); A.op = @(xi,u) diff(u,2) ; % A.lbc = ’neumann’; % A.rbc = ’Robin’; A.lbc = @(u) diff(u,1) % neumann Biot = 1.0 A.rbc = @(u) diff(u,1) + Biot* u B = chebop(0,1); B.op = @(xi,u) -u ; %% % Then we find the eigenvalues with eigs. [F,L] = eigs(A,B) omega = sqrt(diag(L)) % eigenvalues % series coefficient can be evaluated i = 1; for i = 1:6 center = F(0,i); N1 = center * int ( F(:,i) ); N2 = int ( F(:,i).* F(:,i) ); B1(i) = int ( F(:,i) )/center A1(i) = N1/N2 % Answer C1 = 1.119 end On running the code you should get the following results for Biot = 1. The eigenvalues are 0.8603, 3.4256, 6.4373, 9.5293, 12.6453, 15.7713, and 18.9024. The series coefficients for Biot = 1 can be calculated as 1.1191, −0.1517, 0.0466, −0.0217, 0.0124, and −0.0080. It is only a matter of summing the series for any chosen time to find the temperature profiles. You can write a small piece of code and generate time–temperature plots for any given Biot number.