using Plots

L=0.24;
t1=4*60*60; #t1 in Sekunden
anz=30; #Zahl der Summanden 

D=0.79; #Ziegel
#D=0.197 #Gasbeton

K=5.27e-7; #Ziegel
#K=4.69e-7; #Gasbeton

C=8*L/(D*pi^2);
h=pi^2*K/(4*L^2);

function u(t,x)
    if t<t1 
        s=(L-x)/D;
        for n=0:anz
            s=s-C/(2*n+1)^2*exp(-h*(2*n+1)^2*t)*cos((2*n+1)/(2*L)*pi*x);
        end
    else 
        s=0;
        for n=0:anz
            f=exp((2*n+1)^2*h*(t1-t))-exp(-(2*n+1)^2*h*t)
            s=s+C/(2*n+1)^2*f*cos((2*n+1)/(2*L)*pi*x);
        end
    end
    return s;
end



t=range(0,4*t1, length=100);
y=map(t->u(t,0), t);
display(plot(t,y, title="Ziegelwand"))
M=5;
for k=1:M
  global x,y;
  x=k/M*L;
  y=map(t->u(t,x), t);
  display(plot!(t,y))
end

ylims!(0,0.5)

savefig("D:Dokumente/Lehre/Aufgabensammlung/Aufgaben/Bilder/ziegelwand.png")