%% Harj.1, teht. 4, H1T4R.m , mlG11R.tex

%% a)
%
x=linspace(0,5);
y=5*cos(3*pi*x);
plot(x,y); shg
%% b)
%
y=1./(1+x.^2);
plot(x,y)
shg
%% c)
%
close all
x=linspace(-2,2,200);
oso=sin(7*x)-sin(5*x);
nimi=cos(7*x)+cos(5*x);
plot(x,oso./nimi)
grid on
title('Y-rajaus on tarpeen')
shg
%%
disp('No niin, nimittäjällähän on 0-kohtia.')
% Katsotaan, missa
figure   % Uusi grafiikkaikkuna
plot(x,nimi,'r')
grid on
shg
%%
figure(1) % Palataan ekaan ikkunaan.
plot(x,oso./nimi)
axis([-2 2 -10 10])   % y-rajaus jättää pois 'räjähdyspisteiden'
%                      ympäristöt
grid on
shg
%% Mutta miksi ei räjähdä muiden nimittäjän 0-kohtien lähellä?
% Piirretään nimittäjäikkunaan myös osoittaja eri värillä.
figure(2)
hold off
plot(x,oso,'b',x,nimi,'r')
legend('Osoittaja','Nimittaja')
grid on
shg
Muissa nimittäjän nollakohdissa osoittaja näyttää olevan myös 0.
%% Taulukoidaan osoittaja ja nimittäjä jälmimmäisen pienuuskohdissa.
ind=find(abs(nimi)<0.1)
[oso(ind)' nimi(ind)']
    
%%


