Info

This question is closed. Reopen it to edit or answer.

How do I solve this differential Equation?

1 view (last 30 days)
Callum Jackson
Callum Jackson on 6 Feb 2015
Closed: MATLAB Answer Bot on 20 Aug 2021
I have defined my variables a, b and c for a differential equation using the code below, where a is a 31x1 matrix, how do I now solve the differential equation: a*x + b*x^1/2 + c = 0 and then plot the solutions of the differential equation with respect to the thrustloading variable?
Code:
thrustloading = transpose(0:0.1:3);
delta = 1;
theta =1.0849;
CLmax=2.4;
kto=1.2;
friction=0.07;
Cdr=0.1;
Cd0=0.014;
K1=0.0375;
Mto=0.1;
tr=3;
beta=1;
gamma=1.4;
TR=1.087;
g=9.81;
sigma = delta/theta;
theta0=theta*(1+((gamma-1)/2)*Mto^2);
delta0=delta*(1+((gamma-1)/2)*Mto^2)^(gamma/(gamma-1));
density=1.225*sigma;
alpha=delta0*(1-0.49*(sqrt(Mto)));
CL=CLmax/kto^2;
Cd=(K1*CL^2)+Cd0;
XI=Cd-friction*CL;
a=-(beta/(sigma*density*g*XI))*log(1-XI./((((alpha/beta)*thrustloading)-friction)*CL))
b=tr*kto*sqrt(2*beta/(sigma*density*CLmax))
c = -3048

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by