El script no me corre
Publicado por Pablo (1 intervención) el 14/07/2018 03:54:18
Hola Les hago una consulta, estoy tratando de correr este script pero me tirar error de dimensión. Y nose como solucionarlo.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
function [X,T,u] =ej3b(a,L,b,N)
dx = 0.02; dt = 0.02;
for n=1:N
a(n) = (4*L*b/(pi*pi)*sin(n*pi/2)*1/(n*n));
end
X = 0:dx:L;
T = 0:dt:3;
u=zeros (length(T) , length(X) );
XX=repmat (X ,[length(T) 1] );
TT=repmat (T',[1 length(X)] );
for n=1:N
u = u + a(n).*cos((n*pi*a*TT)/L).* sin((n*pi*XX)/L);
end
fprintf('X=%f T=%f u[X,T]=%f\n',XX, TT, u) ;
surf(XX,TT,u);
end
Valora esta pregunta
0
