Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d1_gauss_4.m @ 434:0bc37a25ed88

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bug fix sbp in time gauss default arg
author Martin Almquist <martin.almquist@it.uu.se>
date Fri, 17 Feb 2017 15:08:25 +0100
parents 42c4f0b545d6
children
comparison
equal deleted inserted replaced
421:b6a5dc423990 434:0bc37a25ed88
1 function [D1,H,x,h,e_l,e_r] = d1_gauss_4(L) 1 function [D1,H,x,h,e_l,e_r] = d1_gauss_4(L)
2 2
3 % L: Domain length 3 % L: Domain length
4 % N: Number of grid points 4 default_arg('L',1);
5 if(nargin < 2)
6 L = 1;
7 end
8 5
9 N = 4; 6 N = 4;
10 7
11 % Quadrature nodes on interval [-1, 1] 8 % Quadrature nodes on interval [-1, 1]
12 x = [ -0.8611363115940526; -0.3399810435848563; 0.3399810435848563; 0.8611363115940526]; 9 x = [ -0.8611363115940526; -0.3399810435848563; 0.3399810435848563; 0.8611363115940526];