% numsym/05/H/harj3teht1.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Jatkuvan tehtävän ominaisfkt: 
% v(k,x)=sin(beta(k)*x), missä beta(k)=k*pi

clear all; close all
n=9;  % Muuttele tätä!

x=linspace(0,1); x=x';
v=zeros(length(x),n);
for k=1:n
  beta=k*pi;
  v(:,k)=sin(beta*x);
end;

plot(x,v)
hold on

%%%%%%%%%%%%%%%%%%%%%%%%%
% Diskreetti tehtävä:
%
X=linspace(0,1,n+2);
%X=X(2:end-1);   % Sisäpisteet (annetaan päätepisteidenkin olla)
K=1:n;
V=sin(pi*K'*X);
plot(X,V,'o')
title(['n= ',num2str(n)]) ;grid on; shg
%
X=X'; V=V'; % Jos halutaan operoida sarakkeittain, kuten edellä,
            % täytyy transpoida. Toinen vaihtoehto:
            % plot(x,v(:,1),'b',X,V(1,:),'o') 
            % ...      %%%         %%%%         
%%%%%%%%%%%%%%%%%%%
% Eri ikkunoihin:   Viritetään ensin tapaukseen n=4.
figure
subplot(221)
     plot(x,v(:,1),'b',X,V(:,1),'o')
     title('k=1')
     grid on
subplot(222)
     plot(x,v(:,2),'r',X,V(:,2),'o')
     title('k=2')
     grid on
subplot(223)
     plot(x,v(:,3),'g',X,V(:,3),'o')
     title('k=3')
     grid on
subplot(224)
     plot(x,v(:,4),'c',X,V(:,4),'o')
     title('k=4')
     grid on
shg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Yleisemmin:

clf
m=ceil(sqrt(n));  % Joudutaan piirtämään "liikaa", ellei
                  % n ole neliö.
for k=1:n
subplot(m,m,k)
  plot(x,v(:,k),X,V(:,k),'o')
  grid on
  title(['k= ',num2str(k)])
end;
shg
