Let's suppose I know every coefficients this script:(adsbygoogle = window.adsbygoogle || []).push({});

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

d_alfa=(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);

dy(2)=(q_inf*y(3)*sin(y(1))*cos(y(1))*d_alfa...

-g*y(3)*(rho-rho_inf)*sin(y(1))+...

E*U_inf*cos(y(1))-...

pi*h*rho*beta*(y(2)-U_inf*cos(y(1)))^2-...

y(2)*E)/(rho*y(2)*y(3)); %dv

d_v= (q_inf*y(3)*sin(y(1))*cos(y(1))*d_alfa...

-g*y(3)*(rho-rho_inf)*sin(y(1))+...

E*U_inf*cos(y(1))-...

pi*h*rho*beta*(y(2)-U_inf*cos(y(1)))^2-...

y(2)*E)/(rho*y(2)*y(3));

dy(3)=(E-y(3)*rho*d_v)/(rho*y(2)); %dA

d_A=(E-y(3)*rho*d_v)/(rho*y(2));

dy(4)=(E*Cp*y(4)+Nu/d*k*C*(T_inf-y(4))...

-rho*y(2)*Cp*y(4)*d_A-rho*y(3)*Cp*y(4)*d_v)...

/(rho*y(3)*Cp*y(2)); % dT

This script solves the system.jpg?

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Solving ODEs system at Matlab

Loading...

Similar Threads for Solving ODEs system |
---|

MATLAB A second order PDE |

Maple ODE System solution Maple |

Mathematica Solve Perturbation problem with mathematica |

MATLAB Solving equation with integration using MATLAB |

**Physics Forums | Science Articles, Homework Help, Discussion**