singular Jacobian encountered bvp4c

1 次查看（过去 30 天）

pxg882 2012-12-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/55751-singular-jacobian-encountered-bvp4c

I'm following the bvp4c example 5 found here:

http://www.mathworks.co.uk/matlabcentral/fileexchange/3819-tutorial-on-solving-bvps-with-bvp4c/content/BVP_tutorial/BVP_examples/ex5bvp.m

I have tried to modify the code for my problem, however I receive 'a singular Jacobian encountered' message when I try to run my code.

This is what I have:

    function nNbvp
    etainf = 55;
    options = bvpset('RelTol',1e-6);
    solinit = bvpinit(linspace(0,etainf),[0 0.506869 1 -0.624297 0]);
    sol = bvp4c(@nN,@nNbc,solinit,options);
    eta = sol.x;
    f = sol.y;
    fprintf('\n');
    fprintf('Value of H(etainf) is = %7.5f.\n',f(5,end))
    clf reset
    hold on
    plot(eta,f(1,:));
    plot(eta,f(3,:));
    plot(eta,f(5,:));
    axis([0 55 0 2]);
    xlabel('\eta')
    hold off
    shg
% --------------------------------------------------------------------------
function Yprime = nN(eta,Y)  
n = 0.9;
X(1) = (Y(2)^2+Y(4)^2)^((1-n)/2);
X(2) = Y(5)+((1-n)/(n+1))*eta*Y(1);
X(3) = (n-1)/(Y(2)^(2)+Y(4)^(2));
X(4) = Y(2)*Y(4)*X(3);
Yprime = [ Y(1)
           X(1)/n*((Y(1)^2-Y(3)^2+X(2)*Y(2))*(1+X(3)*Y(4)^2)-X(4)*(2*Y(1)*Y(3)+X(2)*Y(4)))
           Y(3)
           X(1)/n*((2*Y(1)*Y(3)+X(2)*Y(4))*(1+X(3)*Y(2)^2)-X(4)*(Y(1)^2-Y(3)^2+X(2)*Y(2)))
           -2*Y(1)-((1-n)/(n+1))*eta*Y(2)];
% --------------------------------------------------------------------------
function res = nNbc(Y0,Yinf)
n = 0.9;
res = [Y0(1)
       Y0(3)-1
       Y0(5)
       Yinf(2) - (n/(inf*(n-1)))*Yinf(1)
       Yinf(4) - (n/(inf*(n-1)))*Yinf(3)];