Mein MATLAB Forum

//träger Hebelarm_euler_polygonzug
//Konstanten
theta=1;
l=0.5616;
m=13.1;
phim=0.34907;
fg=68.67;
dt=0.06;
fm=18000;
n=l/5;
y2(1)=0.01;
y1(1)=7;
ge=9.81;
h=0.1;
//Funktionen
function [zeit]=t(w)
    zeit=(w-1)*h;
endfunction

function [kf]=ef(i)
    kf=(-sin(y2(i))*l*l*m*cos(phim))/(theta+cos(y2(i))*l*l*m*cos(phim));
endfunction

function [gf]=eg(i,u)
    gf=(m*ge*cos(phim)*l-(fm/dt)*t(u)*cos(y2(i))*n+fg*0.5*cos(y2(i))*l)/(theta+cos(y2(i))*l*l*m*cos(phim));
endfunction


for i=2:10
    y1(i)=y1(i-1)+h*(eg(i-1,i)-ef(y2(i-1))*y1(i-1)*y1(i-1));
    y2(i)=y2(i-1)+h*y1(i-1);
end

-->for i=2:10
-->    y1(i)=y1(i-1)+h*(eg(i-1,i)-ef(y2(i-1))*y1(i-1)*y1(i-1));
-->    y2(i)=y2(i-1)+h*y1(i-1);
-->end
 !--error 21
ungültiger Index.

at line       2 of function ef called by :  
    y1(i)=y1(i-1)+h*(eg(i-1,i)-ef(y2(i-1))*y1(i-1)*y1(i-1));
at line      32 of exec file called by :    
NTO_SCILAB-5167202354040281151.sce', 1
while executing a callback