- #1

- 34

- 0

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.