//Dies ist ein Kommentar //Simulation eines Federpendels g=9.81; //Variablen, Startwerte und Konstanten deklarieren m=0.2; D=32.7; v=0; y=-0.13; dt=0.01; Y=[];V=[];T=[]; //Speicherplatz reservieren für Variablen; Achtung: Grossbuchst.! for t=0:dt:5, //von 0 bis 5 in dt-Schritten immer neu t, v, s berechnen a=-g-D/m*y; //denn a = F/m , F=-D*y Hooksches Gesetz v=v+a*dt; //d.h. v = v_alt + dv y=y+v*dt; //d.h. y = y_alt + ds V=[V,v]; //die aktuellen t-, y- und v-Werte werden gespeichert, Y=[Y,y]; //bevor sie im naechsten Durchlauf ueberschrieben werden T=[T,t]; end; //der Rest ist graphische Ausgabe clf(); //leeres Ausgabenfenster erstellen subplot(2,1,1) //beide Diagramme in ein Fenster, das erste oben links xgrid();xtitle('t-y-Diagramm', 't/s','y/m'); //mit Gitter und Beschriftungen plot2d(T,Y ,style=1); subplot(2,1,2) //das zweite unten rechts xgrid();xtitle('t-v-Diagramm', 't/s','v in m/s'); plot2d(T,V ,style=2);