Recent content by Allamarein

You are right. I hope to achieve my goal quickly. Early I submit you a new problem! It's very similar..but it is not a system.

Plugging? Progress? Sincerely I think to be in the middle of the ocean and I know where is the coast. Would you like to see my M.files?

I tried to use your wise hints. I used taylor series for trigonometric functions \eqalign{ & \sin \left( \alpha \right) \approx 1 - {1 \over 2}\left( {\alpha - {1 \over 2}\pi } \right)^2 \cr & \cos \left( \alpha \right) \approx - \alpha + {1 \over 2}\pi \cr} I didn't discover a...

I suspected this behavior. In fact I tried to use different solver (ode45, ode113, ode15s,...). In order to learn ODE solving in Matlab I read a lot of tutorials, but they offer easy problems. I consider Matlab should be helpful when the system is complicated (like mine). I wouldn't be a...

Dr Transport seems very stong in this kind of problems. Please be patient with me: it' s the first time I tackle ODEs system and I should join the fray without a solid preparation. Could you suggest a reference to learn how solve this kind of problem? Anyway if I understand I should write a...

It seem a good hints. May you suggest a first guess of a simpler system? How could I linearize my system.jpg? I am not very experienced in this kind of thing. Anyway I am disapponited: I thought Matlab reports these kind of troubles when it solves ODE

d_v=dy(2); I used this expedient because I thought typing dy(2) directly, Matlab uses that dy(2) and not that calculated at precedent step. In a word it is a complicated story, but now I suppose defining d_v is useless. You'd be right, but these are my ICs. I supposed Matlab should give a...

ICs: y1=pi/2; y2=51.6; y3=1.96e-7 y4=300 I'd know my transposition from mathematical problem to Matlab is correct. My results tell me I am wrong.

Let's suppose I know every coefficients this script: and y1 alfa y2 v y3 A y4 T function dy=isaacsimply(s,y) dy = zeros(4,1) global ... dy(1)=(Cd_for*q_inf*h*(sin(y(1)))^2+... g*y(3)*(rho-rho_inf)*cos(y(1))... +E*U_inf*sin(y(3)))/... (-rho*y(3)*y(2)^2); %dalfa...