Shooting Method (Matlab)

edited March 2011 in Science & Tech
Hi, how do I solve this equation (numerically) using Matlab? (can use shooting method)..

20*x^4*y''''[x] - x^4*y'''[x] + 3*x^2*y''[x] - 6*x*y'[x] +
6*y[x] == 0, y'[1] == 0, y'''[1] == 1, y''[5] == -50, y'''[5] == -20


By Mathematica it is just two lines:

sols = NDSolve[{20*x^4*y''''[x] - x^4*y'''[x] + 3*x^2*y''[x] -
6*x*y'[x] + 6*y[x] == 0, y'[1] == 0, y'''[1] == 1,
y''[5] == -50, y'''[5] == -20}, y, x];
Plot[Evaluate[y[x] /. sols], {x, 1, 5}]
Plot[Evaluate[y'[x] /. sols], {x, 1, 5}]

Thank you very much for your help!
Would need a Matlab expert for this..

Comments

Sign In or Register to comment.