function [functionMesh,xMesh,yMesh] = meshGenerator(f,xMin,xMax,yMin,yMax,nuOfEle) % Define the x thingy xVec = linspace(xMin,xMax,nuOfEle); % Define the y thingy yVec = linspace(yMin,yMax,nuOfEle)'; % Define the meshes % Mesh for x xMesh = zeros(nuOfEle,nuOfEle); % Mesh for y yMesh = zeros(nuOfEle,nuOfEle); % Function mesh functionMesh = zeros(nuOfEle,nuOfEle); % Loop through stuff for index = 1:nuOfEle xMesh(index,:) = xVec; yMesh(:,index) = yVec; functionMesh(index,:) = f(xVec,(yVec(index)*ones(1,nuOfEle))); end end