function [node,elem]=mesh_generator

clear all
%initial shape_circle
theta=[0:pi/20:2*pi];
x_temp=12*cos(theta);
y_temp=12*sin(theta);

% remove the last
size_theta=length(theta);
n_node_temp=size_theta - 1;
n_element_temp=n_node_temp;
for i=1:n_node_temp
    node(i,1)=i;
    node(i,2)=x_temp(i);
    node(i,3)=y_temp(i);
end


for i=1:n_element_temp
    elem(i,1)=i;
    elem(i,2)=i;
    elem(i,3)=i+1;
end
elem(n_node_temp,3)=1;

