|
|
 |
| found 21 articles for 0.302 sec |
  |
| Kartik,
There are two ways to accomplish that:
1) Passing additional parameters to PDEPE (and _all_ functions)
Modify appropriate lines in your code to look like these:
sol = pdepe(2,@mfhfun,@mfhic,@mfhbc,x,t,[],a); %% note [] placeholder for
options
function[c,f,s] = mfhfun(x,t,u,DuDx,a) %% note additional argument a in
_all_ functions
function u0 = mfhic(x,a)
function[pl,ql |
|
 |
|
 |
Group: de.soc.politik.misc · Group Profile · Search for Roht in de.soc.politik.misc
Author: STIMME DER WAHRHEIT
Date: May 16, 2008 13:06
I am facing a different problem while using pdepe. function mfh(P) x = [1e-4:1e-4:1e-2]; %% sphere radius from 100 microns to 10 mm%% t = [0:1:600]; %% time steps in sec%% a = P(1); sol = pdepe(2,@mfhfun,@mfhic,@mfhbc,x,t,a); u = sol(:,:,1); %% visualizing the results %% figure; plot(x,u(end,:)); title('Solution at t = 600 sec'); xlabel('Distance x'); ylabel('T'); function
|
| Show full article (1.17Kb) · Show article thread |
|
|
|
1 ·
2 ·
3 ·
next |
|
|
