Matlab - error al ejecutar el teta-metodo

1. function[u, tiempo]= theta_metodo(theta, delta_x, delta_t, T_f,u_0)
2. J=1/delta_x;
3. N=1/delta_t;
4. nu=delta_t*T_f/(delta_x*delta_x);
Estas son las primeras lineas de mi programa y al ejecutarlo me da el error:
??? Input argument "delta_x" is undefined.

Error in ==> theta_metodo at 2
J=1/delta_x;
Lo que no logro identificar el error. Me puede por favor ayudar???

Cómo lo estas ejecutando?. Podrias poner al foro la forma que los ejecutando.


Saludos.
JOSE JEREMÍAS CABALLERO
Asesorías en Matlab
programador en matlab
Servicios de programación matlab
[email protected]

http://matlabcaballero.blogspot.com

http://www.lawebdelprogramador.com/foros/Matlab/1371532-FORMA_DE_APRENDER_MATLAB.html

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

function[u, tiempo]= theta_metodo(theta, delta_x, delta_t, T_f,u_0)

if nargin <1

    theta=0.6

    delta_x=3;

    delta_t=0.8;

    T_f=3;

    u_0=12;

end

J=1/delta_x;

N=1/delta_t;

u=J/N;

nu=delta_t*T_f/(delta_x*delta_x);

tiempo=nu*delta_t;

Asi es como le estoy ejecutando:

function[u, tiempo]= theta_metodo(theta, delta_x, delta_t, T_f,u_0)
N=1/delta_t;
J=1/delta_x;
nu=delta_t*T_f/(delta_x*delta_x);


u=u_0(2:J); %Por ahora omitimos las condiciones de contorno de Dirichlet

e=ones(J-1,1);
ed=e*(1+2*theta*nu);
eu=e*(-theta*nu);
el=eu;
A=spdiags([el ed eu], -1:1, J-1,J-1);

ed=e*(1-2*(1-theta)*nu);
eu=e*((1-theta)*nu);
el=eu;
B=spdiags([el ed eu], -1:1, J-1,J-1);

tic;
for i=1:N
u=A\(B*u);
end

%Imponemos las condiciones de contorno.
u=[0; u; 0];
tiempo=toc;
end
El problema me pide para ejecutarlo con T=0.6 teta=0,1/2 e 1; nu=1/2, 1/2 e 5; mu=1/40, 1/4 respectivamente.
En realidad soy nuevo con MATLAB, no he aprendido en la carrera porque en mi Pais no se da y solamente lo encontre en el master y a veces tardo entender las cosas. Estare muy agradecido por vuestra ayuda.
Por ahora el error que me esta dando es:
??? Input argument "delta_t" is undefined.

Error in ==> theta_metodo at 2
N=1/delta_t;

Atentamente,
A. Onesimo

Para ser más claro con mi pregunta he decidido ponerte el problema que estoy intentando resolver, tal cual lo encontré:
Consideremos la ecuacion del calor
ut=uxx, (x,t)ϵ(0,1)x(0,T)
u(x,0) = u(1,t)=0, t>0
u(x,0)=u0(x), 0≤x≤1

Com condicion inicial
u0(x)=x(1 - x)

Vamos a obtener aproximaciones a la solucion utilizando el teta-método con teta=0 (método explicito), teta=1/2 (método de Crank-Nicolson) y teta=1 (método implícito). Escogemos tiempo final T=0.6 y consideremos mallados uniformes en la variable x com espaciado ∆x y en la variable t com espaciado ∆t. Llamemos J al número de nodos en la variable x y N al número de nodos en la variable t. Denotemos por
nu=∆t/(∆x)^2, mu=∆t/∆x
• Realice un gráfico en escala doblemente logarítmica representando J frente al error en norma infinito en tiempo T=0.6 com los 3 metodos para los valores:
Ө=0, ѵ=1/2, ѵ=1/6
Ө=1/2, ѵ=1/2, µ=1/40
Ө=1, ѵ=5, µ=1/4

• Realice un gráfico en escala doblemente logarítmica representando el tiempo de CPU necessário para obtener cada una de las aproximaciones frente al error en norma infinito en tiempo T=0.6 para los mismos valores del apartado anterior.

Bueno es programación nivel intermedio en Matlab. Deberías empezar por ejemplos mas sencillos para que veas como función archivos function, tu programa es tipo function en matlab

Si desafortunadamente es lo que estamos estudiando, si tener aprendido antes lo más basico. Mientras tanto estoy estudiando el basico desde los apuntes de tu blog "aprenda matlab 7.0". Por lo tanto si usted me da una aclaración sobre esta práctica sería muy agradecida.
Alcides

Hola Onesimo;

Me da la impresión que el error que obtienes es porque estas tratando de ejecutar directamente la función, en vez de usar o llamar a la función desde un script.

Espero que sea de alguna ayuda.

Saludos
Dave Correa
[email protected]
Servicios de Programación Matlab
http://fismatlab.blogspot.com/