function createArtisynthStiffness(centre,radius,c10stiff) outputfile = 'StiffnessParameters.csv'; nodeFile = 'nodePosition.csv'; elementfile = 'Tongue Elements.elem'; c10= 1037; c20= 486; c11= 0; c02= 0; c01= 0; % [node,distance,inside_nodes]=findnodedistances(centre,radius,nodeFile); % weights = findelements(inside_nodes,elementfile); % weights=[]; for i =1:740 [~,w] = findelementdistance(centre,radius,i); weights = [weights;w]; end z=zeros(740,1); elements = (1:740)'; T = [elements,c10 + weights*(c10stiff-c10),c20+z,c11+z,c02+z,c01+z]; writematrix(T,outputfile)