 34
 0
Main Question or Discussion Point
Hey, guys
Can you please help me to spot mistakes in numerical solution of following diffusivty equation:
∂P/∂t= 0.001127*k/(μ*ϕ*c_t )*((∂P/∂x)^2*c+(∂^2 P)/(∂x^2 )).
Matlab give the following command:
Undefined function or variable 'r'.
Error in function_handle2 (line 9)
for i=1:r
Error in odearguments (line 88)
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.
Error in ode15s (line 149)
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...
Error in Diffusivity_equation (line 41)
[t,p]=ode15s(@function_handle2,tout,p0,options);
I don't know how to interpret the last command. M files containing relevant commands are attached.
Can you please help me to spot mistakes in numerical solution of following diffusivty equation:
∂P/∂t= 0.001127*k/(μ*ϕ*c_t )*((∂P/∂x)^2*c+(∂^2 P)/(∂x^2 )).
Matlab give the following command:
Undefined function or variable 'r'.
Error in function_handle2 (line 9)
for i=1:r
Error in odearguments (line 88)
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.
Error in ode15s (line 149)
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...
Error in Diffusivity_equation (line 41)
[t,p]=ode15s(@function_handle2,tout,p0,options);
I don't know how to interpret the last command. M files containing relevant commands are attached.
Attachments

879 bytes Views: 306

479 bytes Views: 346