function cmatrix = ccoeffunction(region,state)
n1 = 9;
nr = numel(region.x);
cmatrix = zeros(n1,nr);
cmatrix(1,:) = ones(1,nr);
cmatrix(2,:) = 2*ones(1,nr);
cmatrix(3,:) = 8*ones(1,nr);
cmatrix(4,:) = 1+region.x.^2 + region.y.^2;
cmatrix(5,:) = state.u(2,:)./(1 + state.u(1,:).^2 + state.u(3,:).^2);
cmatrix(6,:) = cmatrix(4,:);
cmatrix(7,:) = 5*region.subdomain;
cmatrix(8,:) = -ones(1,nr);
cmatrix(9,:) = cmatrix(7,:);
I created a function for c coefficient in PDE toolbox using the above example given in MATLAB documentation. My problem is a system of parabolic equations. I want to make sure that parabolic non linear system of PDE can be solved using PDE toolbox in MATLAB. However when I write the state.u as mentioned above error shows up: Function handle specifying a coefficient must accept two input arguments and return one output argument.
However, when I solve a single non linear parabolic equation, I am getting results.
Kindly help me with a system of parabolic non linear PDE to be solved using PDE toolbox command line.
Thank you Luis for your reply. PDEPE solves PDE in one dimension. My problem is in 2 dimensions, i.e., x and y.
I am unable to solve nonlinear parabolic system of partial differential equations in 2D using solvepde. However, I am getting some results for a single nonlinear parabolic PDE (2 dimension) using solvepde. I want to make sure that I can solve such a system of equations using solvepde in MATLAB.
Also , I found the following file as a reference, where it is clearly mentioned that MATLAB’s PDE Toolbox does not have an option for solving nonlinear parabolic PDE, we can make use of its tools to develop short M-files that will solve such equations, (refer page 21).
http://www.math.tamu.edu/~phoward/m401/pdemat.pdf
This document dates back to 2010.
- Badges
- Science topic