As per my understanding, you are trying to find the set of solutions that satisfies the four set of equations specified by you.
From the graph attached by you, it is clearly seen that there are many solutions that satisfies the equations. For an example,
t = 1.2720733281966199317376699591669
Satisfies all the equations, but still it is not shown by the function “vpasolve”.
refer to the documentation for “vpasolve”:
It seems that the number of equations should be equal to the number of symbolic variables. In your case, you are using only one symbolic variable “t” and number of equations are equal to 4. As per my understanding, you should pass only one equation in the “vpasolve” function.
I found a workaround to solve the set of equations mentioned by you. I have found it to be working and attaching it at the end of this answer for your reference. (Note: the curve for the equation (dv=0) is plotted in "blue" and equation (d2v=0) is plotted in "red".
The explanation of the code is given below:
Divide the time interval into smaller intervals. For each interval, first solve the equation “eqn1” and find the point (values of “t”) that satisfies it. The solution should lie between the interval specified by “eqn3” and “eqn4”. As a next step, find the value of “d2v” variable for the evaluated solution. If the value of “d2v” is greater than 0, that means we have found one of the solution of equations. Therefore, store them in a vector. Keep on repeating this procedure for all the time intervals.
In the code, I have defined the time range as : 0 < t < 6, and also plotted the equations "dv=0" and "d2v=0" for your reference.
v = 1/2*(2*10^3)*(0.5*cos(t))^2 - 0.75*9.81*20*cos(t);
plot(time,subs(dv,t,time),'b.');
plot(time,subs(d2v,t,time),'r.');
while (starting_value<ending_value)
sol = vpasolve(eqn1,t,[starting_value,(starting_value+interval_length)]);
value_of_d2v = subs(d2v,t,sol);
solution_vector(end+1)=sol;
starting_value=starting_value+interval_length;
disp("The solution (values of t) which satisfies all the equations are given below:");
The solution (values of t) which satisfies all the equations are given below:
Additionally, please refer to the below documentation to understand the functions used in the code:
I hope it helps !
