yexact=1

yexact =

     1

m=[1,2,4,8]

m =

     1     2     4     8

y=[17,9,5,3]

y =

    17     9     5     3

err=abs(yexact-y)

err =

    16     8     4     2

loglog(m,err)
loglog(m,err,'-*')
loglog(m,err,'-*',m,m.^(-1))
loglog(m,err,'-*',m,m.^(-2))
loglog(m,err,'-*',m,m.^(-1))
yrif=2

yrif =

     2

errrif=abs(y-yrif)

errrif =

    15     7     3     1

loglog(m,err,'-*',m,m.^(-1),m,errrif)
loglog(m,err,'-*',m,m.^(-1),m,errrif,'-o')
y=[y,2]

y =

    17     9     5     3     2

m=[m,16]

m =

     1     2     4     8    16

errrif=abs(y-yrif)

errrif =

    15     7     3     1     0

loglog(m,errrif,'-o')
errrif

errrif =

    15     7     3     1     0

LS
Warning: Could not find an exact (case-sensitive)
match for 'LS'.
/usr/local/matlab_R2008A/toolbox/matlab/general/ls.m
is a case-insensitive match and will be used instead.
You can improve the performance of your code by using
exact
name matches and we therefore recommend that you
update your
usage accordingly. Alternatively, you can disable
this warning using
warning('off','MATLAB:dispatcher:InexactCaseMatch').
This warning will become an error in future releases.
2010-03-11.txt	lejapoints.m   octave	      Video
aa0910		lejapoints.m~  oldlab.tar.gz
Documenti	Modelli        Pubblici
Immagini	Musica	       Scrivania

ls
2010-03-11.txt	lejapoints.m   octave	      Video
aa0910		lejapoints.m~  oldlab.tar.gz
Documenti	Modelli        Pubblici
Immagini	Musica	       Scrivania

cd aa0910/equazioni_differenziali/
ls
2010-03-04.txt	   pendolo.m
functionhandles.m  pendolo.m~

pendolo

ans =

    0.0779

uiopen('/home/accounts/personale/clrmrc90/aa0910/equazioni_differenziali/functionhandles.m', true);
functionhandles
??? Error using ==> mpower
Matrix must be square.

Error in ==> @(x,n)x^n


Error in ==> @(x)f(x,n)


Error in ==> <a href="matlab: opentoline('/usr/local/matlab_R2008A/toolbox/matlab/funfun/quad.m',77,0)">quad at 77</a>
y = f(x, varargin{:});

Error in ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa0910/equazioni_differenziali/functionhandles.m',4,0)">functionhandles at 4</a>
quad(@(x) f(x,n),0,1)

functionhandles

ans =

    0.5000


ans =

    0.3333

diary off