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

 Akzeptierte Antwort

madhan ravi
madhan ravi am 9 Sep. 2018
Bearbeitet: madhan ravi am 9 Sep. 2018

0 Stimmen

It’s working for me . Make sure you have partial differential equations toolbox by typing ver in your command window.

19 Kommentare

MINATI
MINATI am 9 Sep. 2018
Thanks Madhan for your quick response By running the code the following Error shows Attempt to execute SCRIPT pdepe as a function: C:\Users\HP\pdepe.m
Error in pdex4 (line 11) sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
madhan ravi
madhan ravi am 9 Sep. 2018
Bearbeitet: madhan ravi am 9 Sep. 2018
change your function filename as pdex4.m
madhan ravi
madhan ravi am 9 Sep. 2018
If it’s according to your wish and working please accept the answer.
madhan ravi
madhan ravi am 9 Sep. 2018
Bearbeitet: madhan ravi am 9 Sep. 2018
Haven’t heard from you in a while @Pattnaik. If something is not clear let know.
MINATI
MINATI am 10 Sep. 2018
Dear Madhan Sorry for being late but file name is already given as yours. Can you do a favour by sending the modified code at minatipatra456@gmail.com
Thanks
madhan ravi
madhan ravi am 10 Sep. 2018
check your mail
MINATI
MINATI am 11 Sep. 2018
your code is running after download and putting RUN option but in my laptop in the following error occurs
Attempt to execute SCRIPT pdepe as a function: C:\Users\HP\pdepe.m
Error in pdex4 (line 7) sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
What to do
madhan ravi
madhan ravi am 11 Sep. 2018
Bearbeitet: madhan ravi am 11 Sep. 2018
Do you have partial differential equations toolbox? You should have one if you want to get the output. You can check it by typing ver in your command window.
MINATI
MINATI am 11 Sep. 2018
yes the version is present
madhan ravi
madhan ravi am 11 Sep. 2018
Bearbeitet: madhan ravi am 11 Sep. 2018
Just download the attached file,
1)Press the green button in the script file OR
2)Type pdex4 in the command window.
madhan ravi
madhan ravi am 13 Sep. 2018
Did you try what I said?
MINATI
MINATI am 13 Sep. 2018
Yes I hv tried but again the following error comes Attempt to execute SCRIPT pdepe as a function: C:\Users\HP\pdepe.m
Error in pdex4 (line 7) sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
I think something wrong is happened with my software but other codes are running smoothly
MINATI
MINATI am 13 Sep. 2018
Thanks Dear Madhan for your patience. Another thing are you OK with Matlab code to solve COUPLED PDE with TWO variables (y & t ) by using LAPLACE TRANSFORM. I need your help Please reply me on minatipatra456@gmail.com
madhan ravi
madhan ravi am 13 Sep. 2018
Bearbeitet: madhan ravi am 13 Sep. 2018
Please accept my answer if it solved your problem.
Walter Roberson
Walter Roberson am 13 Sep. 2018
You need to rename HP\pdepe.m as it is interfering with MATLAB routine of the same name.
madhan ravi
madhan ravi am 13 Sep. 2018
Bearbeitet: madhan ravi am 13 Sep. 2018
@PRADYUMNA I have no experience with solving PDE
MINATI
MINATI am 13 Sep. 2018
Hi Walter following error occurs Error using pdepe Too many input arguments.
Error in pdepe (line 6) sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
Which MATLAB version are you using? And what shows up for
which pdepe
MINATI
MINATI am 28 Dez. 2018
C:\Users\HP\pdepe.m

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Gefragt:

am 9 Sep. 2018

Kommentiert:

am 28 Dez. 2018

Community Treasure Hunt

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

Start Hunting!

Translated by