Problem 51783. Solve an ODE: equation B
Solution Stats
Problem Comments
-
2 Comments
Chris,
I solved the first test suite symbolically to make sure I was on the right path an got:
[1.0, 6.4839432567518947511762312927844, 9.0049224364695494006414430797622, 8.992494072596270056922480305016, 6.9752731606055591524878412503205, 3.6901871946858786387135688683041, 0].
I am not sure why it would be wrong.
syms y(x)
Dy = diff(y);
a=1;p=0;
X=linspace(1,4,7);
ode=diff(y,x,2)==-Dy/x+(p^2/x^2-a^2)*y;
conds=[y(1)==1, y(4)==0];
ySol(x)=simplify(dsolve(ode,conds));
Y=vpa(ySol(X));
Sorry David. I made a mistake with the signs in the last term of the equation. I "modified" the problem description (quotes intentional). This version should work, but please let me know if it doesn't.
Solution Comments
Show commentsProblem Recent Solvers8
Suggested Problems
-
Recurring Cycle Length (Inspired by Project Euler Problem 26)
159 Solvers
-
Pandigital Multiples of 11 (based on Project Euler 491)
52 Solvers
-
Highly divisible triangular number (inspired by Project Euler 12)
162 Solvers
-
Find terms in the Connell sequence
34 Solvers
-
Pick the die most likely to win
12 Solvers
More from this Author316
Problem Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!