Facebook
From Sexy Macaque, 4 Years ago, written in Plain Text.
Embed
Download Paste or View Raw
Hits: 219
  1. A=zeros(7,3);
  2. calka=c(b)-c(a);
  3. #for m=1:3
  4. for k=1:7
  5.   A(k,1)=asq(a,b,f,10^(-k));
  6.   x(1,k)=log(A(2)); #liczba przedziałów
  7.   x(2,k)=-log(abs(A(1)-calka)); #różnica dla metody asq
  8.   #c(b)-c(a) - wartosc całki policzona analitycznie, trzeba zdefiniować
  9.                        #razem z a,b,f itd
  10.   x(3,k)=abs( nsq(a,b,f,A(2))-calka);   #różnica dla metody NSQ                  
  11.  endfor
  12. #endfor
  13.   plot(A(1,:), A(2,:), ";ASQ", A(1,:), A(3,:), "'NSQ;")