% scaletta per la lezione del 21 aprile 2005
%
% 

clear

% FUNZIONI PER VETTORI E MATRICI PARTICOLARI

v = 1:10

u = 10:-1:1

empty = 10:1

x = 0 : 0.1 : 1

x = linspace(0, 1, 11)

a = 1;
b = 5;
n = 5;

x = linspace(a, b, n);

x = linspace(a, b);

x =  linspace(a, b, n)';



% USO DELLA NOTAZIONE ":" PER SCRIVERE ISTRUZIONI IN MODO COMPATTO

A = [ 1,2,3; 4,5,6; 7,8,9];

A(1, :) = 2:2:6;


% CANCELLARE ELEMENTI

v = 1:10

v(2:5) = []

A(:,1) = []



rig = 3;
col = 5;

xlog = logspace(0.1,1000,5);

ZEROMN = zeros(rig,col);

UNI    = ones(rig,col);

IDENT  = eye(rig);

IDENT2 = eye(col);

DIAGO  = diag(v);

QUADR  = magic(3);




who

whos



% FUNZIONI

lookfor cos

n = 5;
x = linspace(0, pi, n);

coseno = cos(x);
seno   = sin(x);
 
[ x' coseno' seno' ]


% approssimazione della funzione coseno

y = (15120 - 6900*x.^2 + 313*x.^4) ./ (15120 + 660*x.^2 + 13*x.^2);

[ x' y' ]

% creazione di una funzione 

%%%%%%%%%%%%%%%%%%
%function y = apprcos(x)
%y = (15120 - 6900*x.^2 + 313*x.^4) ./ (15120 + 660*x.^2 + 13*x.^2);

% uso delle funzioni

% i file usati sono scaricabili dalla rete:
% http://www.mcgraw-hill.it/naldi

% scompattare i file con un programma apposito:
% per esempio: ark in un sistema linux
%              winzip in un sistema windows

% Cartella filem/cap1


% nota: correggere il file mediav.m fornito
% la parola function all'inizio del file e' incompleta
help mediav

m1 = mediav(v);



help mediam

m2 = mediam(A);


help mediam2

m2bis = mediam(A)


% calcolare le medesime quantita' usando la funzione SUM


% un po' di successioni ricorsive

help successione

a1 = 9;
as = 1e4;

% nota: il file successione.m fornito e' incompleto. 
% Aggiungere un'ulteriore riga:
% x=a;

seq = successione(a1, as)

% operatori relazionali

a1 < as               % minore di

a1 <= 9               % minore di o uguale a

a1 > as               % maggiore

a1 >= as              % maggiore di o uguale a

a1 == as              % uguale a

a1 ~= as              % diverso da



% operatori logici (booleani)

a2 = 4;

a1 < as & a1 ~= as 

a1 < as | a1 ~= as

~(a1 < as & a1 ~= as)

xor(a1 ~= as, a2 > as)


help successione2

seq2 = successione2(a1)

seq2 = successione2(100);


% GRAFICA ELEMENTARE

plot(seq2, ':o');



x = linspace(0,pi,101);

y = cos(x);

y_app = apprcos(x);


figure

plot(x,y, 'r', x,y_app,'k')

