Attempt to execute SCRIPT pdepe as a function:

3 次查看(过去 30 天)
Hi dears
I have problem with executing this program. Can anyone explain me what is going on here? /
Attempt to execute SCRIPT pdepe as a function: /Applications/MATLAB_R2019b.app/toolbox/matlab/funfun/pdepe.m Error in kpz (line 7) sol = pdepe(m,@heat1D,@heatIC,@heatBC,x,t);
opts = odeset('RelTol',1e-8,'AbsTol',1e-8);
N=200; xi=0; xf=5; ti=0; tf=5; TIME=tf-ti;
x = linspace(xi,xf,N);
t = linspace(ti,tf,200);
m = 0;
tic
sol = pdepe(m,@heat1D,@heatIC,@heatBC,x,t);
toc
w = sol(:,:,1);
surf(x,t,w)
xlabel('x')
ylabel('t')
zlabel('w(x,t)')
view([40 55])
%plot(w(200,:))
function [c,f,s] = heat1D(x,t,w,dwdx)
D=4; %0.2;
c = 1;
f = D*dwdx;
s = 1.5*dwdx.^2; %% NO NOISE
%s = 1*dwdx.^2 - 0.5/t; %%White
%s = 0.5*dwdx.^2 - 1/(x*sqrt(t)); %%Pink
%s = -0.5*dwdx.^2 - 0.5/(x*x); %%Brown
%s = 10.5*dwdx.^2 - 0.1*exp(-x*x/t)/t; %%Gauss
end
function w0 = heatIC(x)
w0=2;
end
function [pl,ql,pr,qr] = heatBC(xl,wl,xr,wr,t)
pl = wl;
ql = 0;
pr = 0;
qr = 1;
end
  2 个评论
madhan ravi
madhan ravi 2020-3-21
which pdepe -all % what do you get ?
You have a script named "pdepe.m" probably rename or remove it.
Okhunjon Sayfidinov
I will try everything one more time. thanks for your support.

请先登录,再进行评论。

采纳的回答

Cris LaPierre
Cris LaPierre 2020-3-21
Your code runs as is for me. Perhaps you have overwritten the pdepe function? Try clearing it by typing
clear pdepe
Run the code again. If you still get the error, try typing
which pdepe
The result should be something like this
C:\Program Files\MATLAB\R2020a\toolbox\matlab\funfun\pdepe.m
  14 个评论
Okhunjon Sayfidinov
I got error on my Matlab still, then I tried this code to other laptop and it is working now. i guess there is a problem with my laptop.
What do you suggest to solve this problem on my laptop?
For now i will use other loptop to find solution to my equation.
Thanks dear a lot for your support.
Cris LaPierre
Cris LaPierre 2020-3-21
编辑:Cris LaPierre 2020-3-21
This error suggests that MATLAB is finding multiple files called pdepe.m. Let's see what files MATLAB is seeing. What is the result of
which pdepe -all

请先登录,再进行评论。

更多回答(0 个)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by