clear
clc
close all
IgaTop2D(30, 10, [1 1],[101 51], 4,0.3, 3, 2);