C	DIE SUBROUTINE SQUARE WIRD ZUR BERECHNUNG DER FUNKTION 
C	FUER DEN JEWEILIGEN WIDERSTAND BENUTZT

	SUBROUTINE SQUARE(RTC1,RTC2,RTC3,RTM,RTB,RTH)

	RL1=RTC2-RTC1
	RL2=RTC3-RTC2

	IF (RL1) 30,10,20
10	IF (RL2) 20,15,20
15	RTM=0.
	RTB=0.
	RTH=0.
	RETURN

20	IF (RL1) 30,40,40
30	WRITE(1,35)
35	FORMAT('**FEHLER: C2-C1<0 ')
	STOP

40	IF (RL2) 50,60,60
50	RL2=0.
60	RTH=1./(RL2+RL1/2.)

	IF (RL1) 30,70,80
70	RTM=0.
	GOTO 90
80	RTM=RTH/RL1
90	RTB=-RTM*RTC1
	RETURN

	END