z[0]<-input$valz
k[0]<-trunc((input$vala*input$valz+input$valc)/input$valm)

for (i in 1:input$num-1) {

  k[i-1]<-trunc((input$vala*z[i-1]+input$valc)/input$valm)
  z[i]<-(input$vala*z[i-1]+input$valc)-input$valm*k[i-1]
  u[i]<-(z[i]/input$valm)
}