####Aufghbe 41
h=rbind(c(810,796,8621,4938,2895),c(2107,1859,23304,10767,5074),c(1635,1648,22557,10020,5422))# Erzeugen der Tabelle als Matrix
# alternativ:  h= matrix(data = c(810,796,8621,4938,2895,2107,1859,23304,10767,5074,1635,1648,22557,10020,5422), nrow = 3, ncol = 5, byrow = TRUE)
            


print("          h:   ")
print(h)                     #Ausgabe von h


Z_altstadt_lehel=sum(h[1,])                # berechnen aller Zeilensummen
Z_ludwigvorstadt_ishrvorstadt=sum(h[2,])
Z_maxvorstadt=sum(h[3,])
S_kleiner6=sum(h[,1])                      # berechnen aller Spaltensummen
S_615=sum(h[,2])
S_1545=sum(h[,3])
S_4565=sum(h[,4])
S_groesser65=sum(h[,5])




Zeilensummen=rbind(sum(h[1,]),sum(h[2,]),sum(h[3,]))     # Erzeugen eines Spaltenvektors,der die Zeilensummen enthält
Spaltensummen=cbind(sum(h[,1]),sum(h[,2]),sum(h[,3]),sum(h[,4]),sum(h[,5])) # Erzeugen eines Zeilenvektors, der die Spaltensummen enthält

N=sum(h)       # Gesamtanzahl aller erfassten Bürger


Randverteilung=rbind(cbind(h,Zeilensummen) ,c(Spaltensummen,sum(Spaltensummen))) # Zusammensetzung der Gesamtmatrix
rownames(Randverteilung)=c("Altstadt-Lehel","Ludwigvorstadt-Isarvorstadt","Maxvorstadt","") # Vergabe von Zeilennamen
colnames(Randverteilung)=c("<6"," 6 - 15"," 15 - 45"," 45-65"," >= 65","") #Vergabe von Spaltennamen


print("      Randverteilung:")    
print(Randverteilung)              #Ausgabe der Randverteilung
print("Kontingenztafel der relativen Häufigkeiten:")
print(round(Randverteilung/N,digits=3))   #Ausgabe der Kontingenztafel der relativen Häufigkeiten: alle Matrixeinträge werden durch N dividiert