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.
Problem Recent Solvers8
Suggested Problems
-
1777 Solvers
-
537 Solvers
-
Convert Roman to Arabic Numerals
84 Solvers
-
Magic is simple (for beginners)
5238 Solvers
-
Circular Primes (based on Project Euler, problem 35)
280 Solvers
More from this Author150
Problem Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!