Can't run my main Matlab file related to ODE solution

  • MATLAB
  • Thread starter Ein Krieger
  • Start date
  • #1
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.
 

Attachments

  • Diffusivity_equation.m
    879 bytes · Views: 457
  • function_handle2.m
    479 bytes · Views: 499

Answers and Replies

  • #2
DrClaude
Mentor
7,627
4,068
I get an error on line 10 of function_handle2. When I comment it out, I get an error because dss004 is undefined.
 
  • #3
I get an error on line 10 of function_handle2. When I comment it out, I get an error because dss004 is undefined.

Hey, mate

You can download dss004 m file from the following link:
https://www.physicsforums.com/showthread.php?p=4292439#post4292439.

I couldn't upload the dss004 file, so decided to redirect you to my previous post where I have already atttached dss004.
 
  • #4
Guys, I need to solve by Matlab the following PDE:
4e5c3bb7979f.png


I am desperate about finding solution; I have tried many approaches, and the files attached on topic start is my last approach.
 

Related Threads on Can't run my main Matlab file related to ODE solution

Replies
2
Views
3K
Replies
1
Views
2K
  • Last Post
Replies
1
Views
2K
Replies
7
Views
3K
  • Last Post
Replies
1
Views
2K
Replies
1
Views
6K
  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
3
Views
12K
Replies
5
Views
78K
Top