Heated-Rod Equation Recursive Question
조회 수: 2 (최근 30일)
이전 댓글 표시
investigating the heat rod function.
Length = 50 cm uniform intiial temperature u0 = 100dC and thermal diffusivity k = 0.15. I was given the following code in the text, and it is fine until I put x into the sine term.
it says my error is in line 8. I know it is the x.
How can I get the correct input arguements? Also how would I get the graphs of u(x, 30) abd u(25,t)?
댓글 수: 0
채택된 답변
Shadaab Siddiqie
2020년 12월 9일
It looks like you want to create function in a file and run it. Refere Create function in files for more information. Here is the code
function u = expand(x,t)
k = .15;
L = 50;
u0 = 100;
S = 0;
N= 50;
for n = 1:2:2*N+1
S = S + 1/n*exp((-n^2)*(pi^2)*k*t/(L^2))*sin((n*pi*x)/L);
end
u=(4*u0*S)/pi;
end
To plot you can run
a1 = expand(1:30,30);
a2 = expand(25,1:30);
plot(a1);
plot(a2);
from the command line or from another file.
댓글 수: 0
추가 답변 (0개)
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!