Mercurial > repos > public > sbplib

view convergence.m @ 407:e1db62d14835 feature/SBPInTimeGauss

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed bug with initial data and made gauss the default type
author Martin Almquist <martin.almquist@it.uu.se>
date Thu, 02 Feb 2017 19:30:16 +0100
parents ad1a41d135a1
children
line wrap: on
line source

function q = convergence(e,h)
    for i = 1:length(e)-1
        q(i) = log(e(i)/e(i+1))/log(h(i+1)/h(i));
    end
end