Why the given code is not running

1 visualización (últimos 30 días)
MINATI
MINATI el 9 de Sept. de 2018
Comentada: MINATI el 28 de Dic. de 2018
function pdex4
m = 0;
x = [0 0.005 0.01 0.05 0.1 0.2 0.5 0.7 0.9 0.95 0.99 0.995 1];
t = [0 0.005 0.01 0.05 0.1 0.5 1 1.5 2];
sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
u1 = sol(:,:,1);
u2 = sol(:,:,2);
% --------------------------------------------------------------
function [c,f,s] = pdex4pde(x,t,u,DuDx)
c = [1; 1];
f = [0.024; 0.17] .* DuDx;
y = u(1) - u(2);
F = exp(5.73*y)-exp(-11.47*y);
s = [-F; F];
end
% --------------------------------------------------------------
function u0 = pdex4ic(x);
u0 = [1; 0];
end
% --------------------------------------------------------------
function [pl,ql,pr,qr] = pdex4bc(xl,ul,xr,ur,t)
pl = [0; ul(2)];
ql = [1; 0];
pr = [ur(1)-1; 0];
qr = [0; 1];
end
% figure
surf(x,t,u1)
title('u1(x,t)')
xlabel('Distance x')
ylabel('Time t')
% figure
% surf(x,t,u2)
% title('u2(x,t)')
% xlabel('Distance x')
% ylabel('Time t')
end

Respuesta aceptada

madhan ravi
madhan ravi el 9 de Sept. de 2018
Editada: madhan ravi el 9 de Sept. de 2018
It’s working for me . Make sure you have partial differential equations toolbox by typing ver in your command window.
  19 comentarios
Walter Roberson
Walter Roberson el 13 de Sept. de 2018
Which MATLAB version are you using? And what shows up for
which pdepe
MINATI
MINATI el 28 de Dic. de 2018
C:\Users\HP\pdepe.m

Iniciar sesión para comentar.

Más respuestas (0)

Community Treasure Hunt

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

Start Hunting!

Translated by