Помощь в решении задачи. MATLAB

Нужна помощь в решении задачи, используя Matlab.

1) Есть сетка и в каждой точке надо найти значение C(i,j) и S(i,j).

Условия:

i=0,...,100

j=0,...,100

S(i,j)=S(i,j-1)+h(1-S(i,j-1))*C(i,j-1)

C(i,j)=C(i-1,j-1)-h(1-S(i,j))*C(i,j-1)

h=0,01

Значение С на оси j всегда равно 1, значение S находим по формуле.

Значение S на фронте равно 0, значение С находим по формуле.

Все значения за фронтом равны 0.

во всех остальных точках тоже находим S,C по формуле.

Первые несколько точек,для примера, посчитала руками.

2) После того, как найдены значения в этих точках, нужно построить график как изменяется С и S.

По горизонтальной оси аргумент x(наше i) или t(наше j), а по вертикальной оси - функция C или S. Интервал изменения обеих функций по вертикали от 0 до 1. По горизонтали для x: 0<=x<=1, для t: 0<=t<=100.

Результатом работы программы являются два двумерных массива значений С(x, t) и S(x,t). Каждому значению С и S, например, С=0,99643 соответствует пара значений х и t, которые должны быть получены в пункте 1.

Надо зафиксировать одну переменную (х или t) и строить график зависимости значений функции С или S от другой переменной (t или х).

1. Фиксируем х=1 и строим обычные плоские графики С(1, t) и S(1,t).

2. Фиксируем t=1 и строим обычные плоские графики С(х, 1) и S(х,1) при 0<=x<=1.

3. Фиксируем t=5, t=10, t=20, t=50, t=100 и строим график аналогично t=1.