function [a,b,c]=mirkf4c
c = [0 1/2 1/2 1]';
a  = [zeros(1,4);1/2 zeros(1,3);0 1/2 zeros(1,2);0 0 1 0];
b = [1/6 2/6 2/6 1/6]';
