%% mlNl010R.m  Funktion nollakohdat
close all
format
f=@(x)12*(x-1).*sin(x./(x.^2+0.4*x+.1))
f=@(x)(x-1).*sin(12*x./(x.^2+0.4*x+.1))
x=linspace(-4,4);
plot(x,f(x)); grid on;shg
figure
x=linspace(-1,1,1000);
plot(x,f(x)); grid on;shg
%%
x=linspace(-.4,.1,1000);
plot(x,f(x)); grid on;shg

%% Alkupisteet kuvasta: -0.5, -0.2, 0, 1 (2 viimeista suoraan lausekkeesta) 
x1=fzero(f,-0.5)
%%
x2=fzero(f,-0.2)
%%
x3=fzero(f,0)
%%
x4=fzero(f,1)
%% Lisätään nollakohtarinkulat kuvaan:
hold on
plot([x1 x2 x3 x4],0,'or')
%% Onko tässä kaikki?
% 
hold off
x=linspace(-10,10);
plot(x,x./(x.^2+0.4*x+.1))
grid on
shg
%%
figure
x=linspace(-1,0);
plot(x,x./(x.^2+0.4*x+.1))
grid on
shg
%%
% Sini ei voi saada arvoa 0, kun x < -1 (tällöin sinin argumentti välillä
% (-3,0)), eikä arvoilla x > 0.
% Niinpä muita ei ole. 


