ode45 invalid call to script octave 6.1.0(GUI)
10 views (last 30 days)
Show older comments
I'm trying to run a program of a big system of odes with 43 equations and 165 constants,here is a preview of my code:
script1:I called it imene99
%dn/dt=npoint
function npoint=imene99(t,n,k1,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11,k12,k13,k14,k15,k16,k17,k18,k...)
%chemin réactionnels
npoint= [k1*n(2)*n(42)+k15*n(5)*n(9)+k16*n(5)*n(26)+k25*n(4)*n(26)+k29*n(4)*n(9)+k33*n(4)*n(11)+......)
....
script 2 :I called it call_imene99
I put my constants there
%costantes :
%temperature
T=17000
k1=2.36 *(1e-12)*(T/300)^(-0.29)*exp(17.60/T)
k2 =3.24 *(1e-12) *(T/300)^(-0.66)
k3=3.50 *(1e-12)*(T/300)^(-0.53)* exp(3.20/T)
k156=...
%initial conditions
[IC1,IC2,IC3,IC4,IC5,IC6,IC7,IC8,IC9,IC10...,IC43]=deal(10^19)
IC=[IC1 IC2 IC3 IC4 IC5 IC6 IC7 IC8 IC9 IC10 IC11 ..... IC43]
%time
t0 = 0;
tf = 1E-7;
t1 = 5E-12;
%resolving
[t,x]=ode45('imene99',[0 5E-12 1E-7],IC)
error: invalid call to script C:\Users\lenovo\imene99.m
error: called from
imene99
starting_stepsize at line 53 column 5
ode45 at line 196 column 25
call_imene99 at line 205 column 8
5 Comments
Answers (0)
See Also
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!