Hey guys, I am new to PF. I need to be able to model a stiff differential equation in MatLab. I haven't used MatLab before so I am not really sure how to set the function and boundary conditions for the equation:(adsbygoogle = window.adsbygoogle || []).push({});

y'' + (2/x)*Y' = (.7/x^2)*( (y^(-1/2)) - (.067)((1-y)^(-1/2) )

y(0)=0

y'(1)=1

The zero boundary is undefined and there is a lot of singularities in the solution itself.

So far I have:

function dy = f(x,y)

y = [1 1];

k = .7;

lambda = .67;

xspan = [0 10];

ic = [0 1];

dy=[y(2); (-2/x) * y(2) + (k/(x)^2)*( (y(1)^(-1/2)) - lambda*((1-y(1)))^(-1/2) )];

[x,y] = ode15s(@f,xspan,ic);

plot(x,y(:,1),'-o')

And I am getting a an error in the command console:

??? Input argument "x" is undefined.

Error in ==> f at 7

dy=[y(2); (-2/x) * y(2) + (k/(x)^2)*( (y(1)^(-1/2)) -

lambda*((1-y(1)))^(-1/2) )];

I do not know where to go from here.

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# MATLAB Second Order Differential Equations in MatLab

Can you offer guidance or do you also need help?

Draft saved
Draft deleted

**Physics Forums | Science Articles, Homework Help, Discussion**