differential equations solution doesnt work

1 view (last 30 days)
i have for example
syms x y(x)
M4=(R1*x-p*g*pi*0.0810*(x-l1/2)-0.2626*p*g*pi*(x-K1-l2/2)-0.0630*p*g*pi*(x-K2-l3/2)-p*g*pi*((5*x - 19)*(3025*x^2 - 49225*x + 219217))/4704000)/E4/(pi*(b*x+c)^4/64);
%with all the variables known
and i type
ode = diff(y,x,2) == M4;
conds = [y(K4)== 1.7536e-04 Dy(K4)== -2.6143e-04];
ySol(x) = dsolve(ode,conds);
ySol = vpa(simplify(ySol))*10^(6)
yesterday i could solve the equation but today i try it and it continuous running without stopping
Also when i int M4 the result is different from the result yesterday this one does not have and log in it or fantastic number
how can i fix this?
  6 Comments
Dimitrios Adam
Dimitrios Adam on 23 Nov 2019
Edited: Walter Roberson on 23 Nov 2019
R1=78055; l1 = 1.6; g=9.81; p=7860; K1=l1; K2=3 ; l2=1.4 ; l3=0.8 ;b=-0.5/2.8; c= b*l1-1;
Yes i put it in the place of M4 defined before
Walter Roberson
Walter Roberson on 23 Nov 2019
The expression
I=vpa((R1*x-p*g*pi*0.0810*(x-l1/2)-0.2626*p*g*pi*(x-K1-l2/2)-0.0630*p*g*pi*(x-K2-l3/2)-p*g*pi*((5*x - 19)*(3025*x^2 - 49225*x + 219217));
is missing two ) somewhere. As we are not familiar with the underlying equations, we will have to get you to make the repair.

Sign in to comment.

Answers (0)

Community Treasure Hunt

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

Start Hunting!