Contents

H2teht. 7

%{
Kirjoitetaan funktio 'omapolyder' omaan työhakemistoon tiedostoksi
 omaployder.m
Tätä varten
1) tarkistetaan komentoikkunan hakemisto ja siirrytään
tarvittaessa omaan (vaikkapa matlab-hakemistoon).
2) Avataan FILE-valikosta new m-file, jonka otsikkorivin pitää
olla tällainen. (Lisätään myös alkukommentit, kokeile help omapolyder)
function dp=omapolyder(p)
% dp=polyder(p) laskee polynomin p derivaattapolynomin(kerroinvektorin).
% Lisäksi piirretään ...
%}

Kommenttimerkkien '%{ %}' välissä oleva koodi tiedostoon omapolyder.m

%{
function dp=omapolyder(p)
% dp=polyder(p) laskee polynomin p derivaattapolynomin(kerroinvektorin).
% Lisäksi piirretään ...
% Esim:
% p=[2 0 1 -1 3], omapolyder(p)
n=length(p);
dp=(n-1:-1:1).*p(1:end-1);
% Tästä alkaa piirto, saat täydentää itse.

%}
% Testataan ensin helppiä
help omapolyder
% Sitten testataan funktiota helpin esimerkillä:
p=[2 0 1 -1 3], omapolyder(p)
  dp=polyder(p) laskee polynomin p derivaattapolynomin(kerroinvektorin).
  Lisäksi piirretään ...
  Esim:
  p=[2 0 1 -1 3], omapolyder(p)

p =
     2     0     1    -1     3
ans =
     8     0     2    -1