Uusia
VanhojaTiedosto: 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.