0

SciLab Compex ode

Dear all,

I have a problem with solving system equation. I cannot understand where is my mistake

clear clc function dRho=TwoLevelSystem(t, rho) dRho(1)=%i*Omega*(rho(3)-rho(4))+Gamma*rho(2); dRho(2)=-%i*Omega*(rho(3)-rho(4))-Gamma*rho(2); dRho(3)=-(gama+%i*Delta)*rho(3)-%i*Omega*(rho(2)-rho(1)); dRho(4)=-(gama-%i*Delta)*rho(4)+%i*Omega*(rho(2)-rho(1)); endfunction Omega=2; Delta=0; Gamma=1; gama=Gamma/2; t=0:0.1:10; y=ode('rk',[1;0;0;0],0,t,TwoLevelSystem)

thanks

Scilab

Login to add comment