ode45 invalid call to script octave 6.1.0(GUI)

29 views (last 30 days)
Imene Yed
Imene Yed on 22 May 2021
Commented: Imene Yed on 22 May 2021
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

Sign in to comment.

Answers (0)

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by