plot differential equation need help
1 view (last 30 days)
Show older comments
Hi. I want to plot (q,y) for three values of x=1,2,3 with two boundary conditions for this differential equation?
how can I write a code
3 Comments
Bjorn Gustavsson
on 24 Jun 2019
Solve it analytically - ought to be simple enough to do by hand. That will give you y(x,q), plot for desired values of q a desired points x.
HTH
Accepted Answer
Pullak Barik
on 24 Jun 2019
I tried a way out, tell me if it helps-
clc
syms y(x);
syms q;
Dy = diff(y);
ode = diff(y,x,2) + q*y == 10*q; %The equation to solve
cond1 = y(0) == 0; %Initial condition 1
cond2 = Dy(2) == 0; %Initial condition 2
conds = [cond1 cond2];
ySol(x) = dsolve(ode, conds); %Solution in terms of x and q
ySol = simplify(ySol)
syms f1 f2 f3;
f1(q) = ySol(1) %substituting x = 1
f2(q) = ySol(2) %substituting x = 2
f3(q) = ySol(3) %substituting x = 3
%f1, f2, f3 now contain the required functions, plotting them below for brevity
figure;
hold on;
subplot(3,1,1)
title('x=1')
fplot(f1);
subplot(3,1,2)
fplot(f2);
subplot(3,1,3)
fplot(f3);
hold off;
You may refer to the following links for help in understanding the above code-
1) Usage of syms- syms
2) Solving differential equations in MATLAB- Solve Differential Equation
3) Plotting symbolic equations- fplot
4 Comments
Bjorn Gustavsson
on 24 Jun 2019
Yes, I do understand that you want the curve of y(x=2,q) for a range of values of q. Now, for each and every value of q, your ODE as written will be a second order differential equation with constant coefficients. That type of ODE have a very simple analytical solution, where you do not have to plug in the value of q to get the solution for y. This is, I think, something that you are expected to learn properly, my advice is that you still do this by pen and paper just to understand the concept of "characteristic polynomial" as well as the separation of particular and homogenous solutions and how you adjust the coefficients to satisfy the initial/boundary-conditions. But I can only suggest, it is your learning to do.
More Answers (0)
See Also
Categories
Find more on Calculus in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!