Attempt to execute SCRIPT pdepe as a function:
3 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Okhunjon Sayfidinov
il 21 Mar 2020
Modificato: Cris LaPierre
il 21 Mar 2020
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 Commenti
madhan ravi
il 21 Mar 2020
which pdepe -all % what do you get ?
You have a script named "pdepe.m" probably rename or remove it.
Risposta accettata
Cris LaPierre
il 21 Mar 2020
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 Commenti
Cris LaPierre
il 21 Mar 2020
Modificato: Cris LaPierre
il 21 Mar 2020
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
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Eigenvalue Problems in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!