Harjoitus 2

Uusia

• Tehtäviin 1 ja 5 liittyviä Matlab-komentoja (h2avt1.m)
• harj2av.mws

• Maple-ws html-muodossa Siinä näkyy diffyhtälöiden Maple-käsittelyä.
• Harj. 1 LV Tehtävän 1 Matlab-käsittelyä. Samalla malliesimerkki Matlab-työn html-dokusta.

Alkuviikko (AV)

Loppuviikko (LV)

Eulerin menetelmä funktiota meuler käyttäen

Tiedosto: f.m
----------------
>> type f

function xpr=f(t,x)
xpr=t-x.^2;
%----------------------------
%Vaihtoehtoisesti voidaan tehdä "inline":lla:

  f=inline('t-x.^2','t','x')
  type meuler  % Jos ei löydy, aseta polku!
  hold on;grid;m=2;
  a=0;b=1;ya=1;   % a:alkuhetki, b:loppuhetki, ya:alkuarvo
 [T,Y]=meuler(f,a,b,ya,m);plot(T,Y,T,Y,'o');[T',Y'],m=2*m
 [T,Y]=meuler(f,a,b,ya,m);plot(T,Y,T,Y,'o');[T',Y'],m=2*m
% Toista vaikka nn-iteroinnilla, se on hauskempaa ja antoisampaa kuin
% juuri  mikään muu.