% Skript zur Musterlösung
%
% Numerische Mathematik, SS 2013
%
% Simulation eines Wasserkreislaufes
%
% Ingolf Busch, Juni 2013


% Aufgabe a)

p = 11;
c = [1,1,1,1,1];
S = [0,0,0,200];
t_end = 200;

[t,U,V,W,R] = my_wasser(t_end,p,c,S);

% Ausgabe als Plot

figure
h1=plot(t,U,'ro-'); hold on; 
h2=plot(t,V,'gd-');          
h3=plot(t,W,'bs-');          
h4=plot(t,R,'c^-'); hold off
legend([h1 h2 h3 h4],'U_n','V_n','W_n','R_n');
xlabel('Zeit')
ylabel('Wassermenge')
title('Aufgabe 1a)')

% Ausgabe Endwerte
fprintf(1,'Aufgabe 1a): \n\n');
fprintf(1,'U_end = %8.4f  \n',U(t_end));
fprintf(1,'V_end = %8.4f  \n',V(t_end));
fprintf(1,'W_end = %8.4f  \n',W(t_end));
fprintf(1,'R_end = %8.4f  \n',R(t_end));

% Aufgabe c)

p = 4;
c = [1,1,0.25,0.25,0.375];
S = [0,0,0,100];
t_end = 400;

[t,U,V,W,R] = my_wasser(t_end,p,c,S);

% Ausgabe als Plot

figure
h1=plot(t,U,'ro-'); hold on; 
h2=plot(t,V,'gd-');          
h3=plot(t,W,'bs-');          
h4=plot(t,R,'c^-'); hold off
legend([h1 h2 h3 h4],'U_n','V_n','W_n','R_n');
xlabel('Zeit')
ylabel('Wassermenge')
title('Aufgabe 1b)')

% Ausgabe Endwerte
fprintf(1,'\n Aufgabe 1b): \n\n');
fprintf(1,'p     = %8.4f \n',p);
fprintf(1,'c_1   = %8.4f \n',c(1));
fprintf(1,'c_2   = %8.4f \n',c(2));
fprintf(1,'c_3   = %8.4f \n',c(3));
fprintf(1,'c_4   = %8.4f \n',c(4));
fprintf(1,'c_5   = %8.4f \n\n',c(5));
fprintf(1,'U_end = %8.4f  \n',U(t_end));
fprintf(1,'V_end = %8.4f  \n',V(t_end));
fprintf(1,'W_end = %8.4f  \n',W(t_end));
fprintf(1,'R_end = %8.4f  \n',R(t_end));

% Aufgabe c)

p = 11;
c = [1,1,1,1,1];
S = [0,0,0,100];
t_end = 50;

[t,U,V,W,R] = my_wasser(t_end,p,c,S);

% Ausgabe als Plot

figure
h1=plot([0:t_end],U,'ro-'); hold on; 
h2=plot([0:t_end],V,'gd-');          
h3=plot([0:t_end],W,'bs-');          
h4=plot([0:t_end],R,'c^-'); hold off
legend([h1 h2 h3 h4],'U_n','V_n','W_n','R_n');
xlabel('Zeit')
ylabel('Wassermenge')
title('Aufgabe 1c)')

% Ausgabe Endwerte
fprintf(1,'\n Aufgabe 1c): \n\n');
fprintf(1,'U_end = %8.4f  \n',U(t_end));
fprintf(1,'V_end = %8.4f  \n',V(t_end));
fprintf(1,'W_end = %8.4f  \n',W(t_end));
fprintf(1,'R_end = %8.4f  \n',R(t_end));
