# _____________________________________________________________
# Kommenttirivi: Tama on octave ohjelmani Omalli13.m
# - Plottaaminen jatkuu ..
# _____________________________________________________________
graphics_toolkit('gnuplot')       # Vahentaa grafiikka ongelmia
clear ; clc              # Poistetaan ... Tyhjennetaan ...
x=pi*((0:50)/25.);                # x-vektori
y=sin(x);                         # y-vektori
z=cos(x);                         # z-vektori
e=0.*z+0.15;                      # virhe
# _____________________________________________________________
subplot(6,10,[2 19]);             # 1. plotin tila pienempi
set(gca,'fontsize',5);            # xy-akseli numerojen koko
axis([0.0,2.0*pi,-1.2,1.2]);      # xy-rajat (vahan tilaa)
hold on;                          # plot ei poista edellista
plot1=errorbar(x,y,e,'~');        # plottaa x,y,e 
hold off;                         # uusi plotti poistaisi
# errorbar komennossa loputtomasti vaihtoehtoja: '~'=default
set(plot1,'linestyle','none','marker','o','color','r');
title(['Funktio',' ','{sin x}'],'fontsize',20) # Paateksti
xlabel('x','fontsize',10);                      # x-teksti
ylabel('sin(x)','fontsize',10);                 # y-teksti
text(5.7, 0.7,'(a)','fontsize',15);# Kirjoitetaan jotain
# __________________________________________________________
subplot(6,10,[31 60]);             # 2. plotin tila: suurempi
set(gca,'fontsize',7);             # xy-akseli numerojen koko
axis([-0.5,2.0*pi+0.5,-1.8,1.8]);  # xy-rajat (paljon tilaa)
hold on;                           # plot ei poista edellista
plot2=errorbar(x,z,e);             # plottaa x,z,e 
hold off;                          # uusi plot poistaisi ...
set(plot2,'linestyle','none','marker','o','color','b');
title(['Funktio',' ','{cos x}'],'fontsize',25)  # Paateksti
xlabel('x','fontsize',15);                      # x-teksti
ylabel('sin(x)','fontsize',15);                 # y-teksti
txt=['Paljonko on ','{\Sigma_{i=1}^n a^i}','{, jos n=7 ja a=-1?}'];
text(0.0,1.5,txt,'fontsize',15);   # Kirjoitetaan jotain
# __________________________________________________________
print('Omalli13.pdf','-dpdf')      # Kuva tiedosto
waitforbuttonpress;                # Odota: Nappain poistaa
 
