Задача на Коши, Matlab

[elenapeteva]

Здравейте! Някой би ли могъл да ми помогне с решението на следната задача на Коши:
y''' + y'' + 5y' + 3y = 1
y(0) = 0
y'(0) = -3
y''(0) = 3
Начертайте графика на решението и определете колко корена има уравнението y(x) = 0 в интервала [0;12].

[peyo]

Здравейте! Някой би ли могъл да ми помогне с решението на следната задача на Коши:
y''' + y'' + 5y' + 3y = 1
y(0) = 0
y'(0) = -3
y''(0) = 3
Начертайте графика на решението и определете колко корена има уравнението y(x) = 0 в интервала [0;12].

Gnu Octave - подобно ...

pkg load symbolic

syms y(x)
Dy = diff(y,x);
D2y = diff(y,x,2);

ode = diff(y,x,3) + diff(y,x,2) + 5*diff(y,x,1) + 3*y == 1;
cond1 = y(0) == 0;
cond2 = Dy(0) == -3;
cond3 = D2y(0) == 3;
conds = [cond1 cond2 cond3];

uSol(x) = dsolve(ode,conds)

ezplot(uSol(x),[0,12]);

resdif332434.png (51.22 KiB) Прегледано 424 пъти

Така като гледам 6 корена.

[elenapeteva]

Здравейте! Някой би ли могъл да ми помогне с решението на следната задача на Коши:
y''' + y'' + 5y' + 3y = 1
y(0) = 0
y'(0) = -3
y''(0) = 3
Начертайте графика на решението и определете колко корена има уравнението y(x) = 0 в интервала [0;12].

Gnu Octave - подобно ...

pkg load symbolic

syms y(x)
Dy = diff(y,x);
D2y = diff(y,x,2);

ode = diff(y,x,3) + diff(y,x,2) + 5*diff(y,x,1) + 3*y == 1;
cond1 = y(0) == 0;
cond2 = Dy(0) == -3;
cond3 = D2y(0) == 3;
conds = [cond1 cond2 cond3];

uSol(x) = dsolve(ode,conds)

ezplot(uSol(x),[0,12]);

resdif332434.png

Така като гледам 6 корена.

Благодаря много!