PureBoard

Code: Alles auswählen

Dim x.d(1000)
Dim y.d(1000)
Global Merker1.d
Global Merker2.d
Define test.d
Global  k.i


Procedure.f Funktion(x.f)
Define x.f
Define y.f

;hier die Formel eintragen zum Beispiel y=x*x-4.0

y=x*x-4
;y=x*x-3
;y=-1.5-Log10(x)
;Ende Formel
ProcedureReturn y
EndProcedure


;Startwerte
x(0)=0.0001
x(1)=7.0


For k = 1 To 100
y(k-1)=Funktion(x(k-1))
y(k)=Funktion(x(k))

If y(k-1)*y(k)>0
  Debug "Keine Nullstelle"
  Goto Fertig
EndIf
    
    
    x(k+1)=(x(k-1)*y(k)-x(k)*y(k-1))/(y(k)-y(k-1))
    ;Mittelwert Test
    ;x(k+1)=(x(k+1)+x(k))/2
    y(k+1)= Funktion(x(k+1))

  
  If y(k+1)=0
         Debug "Fertig"
;         Debug y(k+1)
         Debug x(k+1)
    Goto Fertig
  EndIf


  If y(k+1)*y(k-1)<0
      Debug "Ist kleiner"
      Debug x(k)
      Debug x(k+1) 
     x(k)=x(k-1)
     y(k)=y(k-1)
     
   Else
     Merker1=x(k+1)
     Merker2=y(k+1)
     x(k+1)=x(k)
     y(k+1)=y(k)
     x(k)=Merker1
     y(k)=Merker2
     Debug "Ist grösser"
     Debug x(k)
  EndIf

Next k


Fertig:
; Debug y(k-1)
; Debug y(k)
; Debug y(k+1)
; Debug x(k-2)
; Debug x(k-1)
Debug "k"
Debug k
Debug x(k+1)
Debug Funktion(x(k+1))
; Debug x(k+2)
; Debug x(k+3)
; Debug x(k+4)
; Debug x(k+5)
; Debug x(k+6)