schlot=zeros(30,30) ;
for i=1:1:30
schlot(12:30,13:17) = 0.3 ;
schlot(10,14:16) =0.6  % gelber fitzel im krater


%schlot(6:10,14.1:15.4) = 0.3  erste zahl: y-achse, zweite zahl: x-achse
%obere hälfte des vulkanes

schlot(6:15,1:13) = 0.6
schlot(6:15,17:30) = 0.6
schlot(6:15,17:30) = 0.6
schlot(6:9,14:16) = 0.6
schlot(1:5,1:30) = 0.6

%schwarzer teil des kraters
schlot(13:15,5:14) = 0
schlot(13:15,16:25) = 0
schlot(11:13,7:14) = 0
schlot(11:13,16:23) = 0
schlot(8:11,9:10) = 0
schlot(8:11,20:21) = 0


%dunkelroter teil des kraters

schlot(11:12,10:14) = 0.1    
schlot(9:10,10:13) = 0.1
schlot(8:9,10:11) = 0.1
schlot(11:12,16:20) = 0.1
schlot(9:10,17:20) = 0.1
schlot(8:9,19:20) = 0.1




%vulkan restliche helle rote form unter dem boden

schlot(28:30,5:13) = 0.3
schlot(25:27,9:13) = 0.3

schlot(28:30,16:25) = 0.3
schlot(25:27,16:21) = 0.3


colormap(hot)
colorbar

%image(64*schlot)
imagesc(schlot)
hold on
%drawnow
%hold off
end