function my_pde
% rownanie na propagfacje potencjalu w przewodniku o pojemnosci C i
% opornosci R i promieniu r, na koncu 0 wstrzykjemy prad tak ze potencjal
% tam ma postac t.*exp(-a*t), na drugim koncu jest uziemienie
% Dv/Dt =r/(2RC) D/Dx Dv/Dx  po zadzialaniu funkcji, dane o poprzednim
% dtas\
% --

m = 0;
x = linspace(0,1000,100);
t = linspace(0,20000,100);
close all
sol = pdepe(m,@my_pde1pde,@my_pde1ic,@my_pde1bc,x,t);
u1 = sol(:,:,1); % tu zniekke po przestrzeni


figure;
surf(x,t,u1);
title('u1(x,t)');
xlabel('Distance x');
ylabel('Time t');
figure
for i=1:length(t)
plot(x,u1(i,:))
ylim([0 0.1])
%pause(0.1)
drawnow
end
% --------------------------------------------------------------------------

function [c,f,s] = my_pde1pde(x,t,u,DuDx)
% tu sa moje parametry
r=0.8;
R=1;
C=1;

%a tu zwracam postac funkcji c f i s
c = 1;                                 
f = r/(2*R*C) .* DuDx;                   
s = 0;                               

% --------------------------------------------------------------------------

function u0 = my_pde1ic(x)
u0 = 0;%exp(-((x-0.5).^2)/0.01);                                 

% --------------------------------------------------------------------------

function [pl,ql,pr,qr] = my_pde1bc(xl,ul,xr,ur,t)
a=0.05;
pl = ul-t.*exp(-a*t);                                 
ql = 0 ;%// leop-lepiej-1;                             
pr = ur;                            
qr = 0;                                   

