MATLAB: Solving System ODE with ode45

MATLABodeode errorode45system ode

Hi All, I'm trying to solve a system ODE with ode45. I get this error for my code:

Error in ode45 (line 115) odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

Any help is appreciated. I have attached the m.files

Best Answer

With first look, in your function, there are many undefined variables, that are X0, y0, z0, and Ap0.

Related Solutions

MATLAB: Finding difficult to solve this coupled ode equations using ode45

Parameterize your function as below:

 P=10E12;
 b=0.268E-9;
 Q=285;
 M=3;
 K=1.38E-23;
 v=3.6E16;
 sigma_e=200;
 T=923;
 ec=0.04;
 c=7; %provide your value
 tspan=0:0.01:70;
 initial_conditions=[0.01,100]; %change the intial condition according to your values
 [t,y]=ode45(@(t,y)pair1(t,y,P,Q,b,c,v,sigma_e,M,K,T,c,ec),tspan,initial_conditions)   %function call
 plot(t,y(:,1))
 
 
function dy=pair1(t,y,P,Q,b,c,v,sigma_e,M,K,T,c,ec)           %function definition
dy=zeros(2,1);
dy(1)=(P*(b^2)*v*sigma_e*((y(2)^2*sigma_e).^(-1/3))/(M^2))*exp(-Q/K*T)*sinh((sigma_e)/((M*K*T)*(b^3))*(Q/((y(2)).^2*sigma_e)).^(1/3));
dy(2)=(1/(y(2)*ec))*((c^2)-(y(2))^2)*y(1);
end

Best Answer

Related Solutions

MATLAB: Finding difficult to solve this coupled ode equations using ode45

Related Question