Attempt to execute SCRIPT pdepe as a function:
Mostrar comentarios más antiguos
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 comentarios
madhan ravi
el 21 de Mzo. de 2020
which pdepe -all % what do you get ?
You have a script named "pdepe.m" probably rename or remove it.
Okhunjon Sayfidinov
el 21 de Mzo. de 2020
Respuesta aceptada
Más respuestas (0)
Categorías
Más información sobre Numerical Integration and Differential Equations en Centro de ayuda y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!

